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INTRODUCTION 


Background:  Traumatic  brain  injury  (TBI)  has  been  labeled  as  a  “signature  wound”  in  the  anti-terrorism 
wars  in  Iraq  and  Afghanistan.  15%  -  25%  of  surveyed  returning  service  members  have  been  reported  to 
have  possible  long-term  mild  traumatic  brain  injury  (mTBI)  or  concussion.  In  the  civilian  sector,  the 
prolonged  neuro-cognitive  and  functional  symptoms  following  mTBI  affects  over  1.2  million  Americans 
annually.  However,  our  understanding  of  the  neuropathology  of  mTBI  and  its  recovery  process  is  still 
very  limited  due  to  the  lack  of  sensitive  clinical  and  diagnostic  tools.  First,  current  TBI  classification 
scales  are  mainly  based  on  results  from  computed  tomography  (CT)  scans  of  moderate  to  severe  TBI 
patients,  and  thus  can  hardly  be  applied  to  mTBI  cases.  Secondly,  current  laboratory  biological  markers, 
clinical  CT  and  conventional  magnetic  resonance  imaging  (MRI)  are  either  insensitive  or  not  specific  to 
the  subtle  abnormalities  in  mTBI  and  poorly  predict  patient’s  long-term  outcome.  It  is  in  urgent  need  of 
developing  a  set  of  advanced  MR  imaging  biomarkers  that  can:  i)  be  sensitive  enough  to  differentiate 
mTBI  patient  population  from  normal  healthy  ones,  ii)  have  the  potential  for  outcome  prediction;  and  iii) 
assist  in  management  of  mTBI  patients  in  acute  settings. 

Objective/Hypothesis:  The  overall  goal  of  the  current  research  is  to  develop  imaging-based  biomarkers 
to  improve  diagnosis  and  outcome  prediction  of  mTBI.  Advanced  MRI  techniques  reveal  many  more 
details  than  conventional  imaging.  For  example,  susceptibility  weighted  imaging  (SWI)  and  susceptibility 
mapping  (SWIM)  can  detect  and  quantify  temporal  reduction  of  hemorrhagic  lesions  associated  with 
mTBI  patients’  functional  outcome  and  diffusion  tensor  imaging  (DTI)  can  detect  axonal  injury  at  the 
acute  stage  and  possible  changes  over  time.  As  the  next  generation  of  SWI  developed  in  our  lab,  SWIM 
will  be  used  to  quantify  iron  in  microbleeds  and  oxygen  saturation  in  major  veins  throughout  the  brain. 
Our  central  hypothesis  is  that  axonal  injury  (measured  by  DTI)  and  vascular  damage  (detected  as 
hemorrhagic  lesions  by  SWI/SWIM)  are  important  pathologies  in  mTBI  that  are  associated  with  patients’ 
neurocognitive  and  clinical  symptoms  in  their  recovery. 

Study  Design:  Adult  mTBI  patients  will  be  screened  and  enrolled  at  the  acute  stage  (within  24  hours  after 
injury)  from  the  emergency  department  of  our  local  Level-One  Trauma  Center.  Both  up-to-60-minute 
MRI  scan  and  neuropsychological/clinical  assessment  will  be  conducted  in  mTBI  patients.  In  the  subacute 
(1  month)  and  chronic  (6  months)  stages  after  injury,  the  patients  will  be  brought  back  to  repeat  both  the 
MRI  scans  and  neurocognitive  evaluations.  Age/gender/education-matched  healthy  controls  will  be 
recruited  and  followed  up  using  the  same  imaging  and  neuropsychological  protocol  of  evaluation. 

Specific  Aims:  Specific  Aim  1  is  to  assess  whether  the  advanced  MRI  data  (SWI  and  DTI)  acquired 
acutely  (within  24  hours  after  injury)  are  more  sensitive  than  conventional  imaging  (CT  and  cMRI)  in 
detecting  mTBI.  Specific  Aim  2  is  to  determine  at  what  time  point  after  the  injury  (i.e.,  acute,  subacute 
and  chronic)  these  advanced  MR  techniques  can  differentiate  mTBI  patients  from  healthy  controls. 
Specific  Aim  3  is  to  use  susceptibility  weighted  imaging  and  mapping  (SWIM),  our  next  generation  of 
SWI  technique,  to  quantify  the  amount  of  iron  in  microbleeds,  to  monitor  any  changes  (evolution)  of  the 
microbleed  over  time  and  to  monitor  oxygen  saturation  compared  to  normal  volunteers.  And  Specific 
Aim  4  is  to  determine  which  MRI  indices  have  statistically  significant  associations  with  neurocognitive 
outcome  in  mTBI  patients  over  all  time  points. 
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BODY 


Summary  of  the  Statement  of  Work  for  Phase  1 : 

Phase  I  (Years  1-2):  recruit  a  cohort  of  50  mTBI  patients  from  the  emergency  setting  of  our  local  Level- 
One  trauma  center,  and  50  non-TBI  volunteers  to  establish  and  validate  an  image-based  diagnostic 
approach  to  brain  injury. 

a)  Screen  and  enroll  mTBI  patients;  perform  the  MRI  protocol  (conventional  imaging,  DTI  and  SWI); 
conduct  NP  tests  and  clinical  outcome  measures. 

Due  to  the  late  approval  of  IRB,  we  started  enrollment  fairly  late  (March,  2012).  To  date,  we  have 
successfully  recruited  two  mTBI  patients  and  scanned  them  at  the  acute  stage  (within  24  hours  after 
injury). 

b)  Screen  and  enroll  healthy  volunteers  and  perform  the  same  procedures  as  mentioned  above. 

We  have  successfully  enrolled  7  control  subjects  in  the  past  month. 

c)  Follow  up  mTBI  patients  and  controls  at  1 -month  and  6-months  after  injury  using  the  same  imaging 
protocol,  NP  tests  and  outcome  measures. 

We  have  successfully  performed  1 -month  follow  up  in  our  two  patients  and  two  controls. 

d)  Analyze  the  data  and  establish  a  baseline  imaging  database  for  mTBI  patients  and  normals. 

We  have  established  the  quantification  schemes  necessary  to  analyze  this  data  for  DTI  and  SWIM.  We 
are  in  the  process  of  creating  an  internal  database  from  which  we  can  gather  all  relevant  quantitative 
information  such  as  total  lesion  load,  FA  histograms  and  mean  FA. 

e)  Report  preliminary  data  by  presenting  it  to  conferences  or  publish  it  in  TBI  and  imaging  related 
journals. 

We  have  published  2  papers  and  2  book  chapters.  We  have  also  given  9  presentations  related  to 
this  work  during  the  last  year  by  Prof.  Haacke  and  Prof.  Kou.  The  details  are  given  in  the 
reportable  results  section. 

f)  Develop  and  use  SWI  and  SWIM  for  use  in  evaluating  abnormal  vasculature  and  microbleeds  in  mTBI 
patients. 

We  are  delighted  to  report  that  we  have  made  major  progress  in  the  evaluation  of  cerebral  microbleeds 
using  a  new  iterative  SWIM  algorithm.  The  reader  is  referred  to  the  paper  by  J.  Tang  et  al.  “Improving 
susceptibility  mapping  using  a  threshold-based  K-space/image  domain  iterative  reconstruction 
approach.”,  MRM,  June  26,  2012,  epub,  ahead  of  print  for  a  detailed  description  of  this  new  concept. 
However,  the  basic  concept  is  given  below. 

While  we  were  pending  for  the  IRB  approval  of  our  human  data  collection,  we  spent  considerable  time  on 
the  methodology  itself  and  have  made  some  dramatic  progress  in  this  direction.  We  have  just  published  a 
paper  on  an  advanced  form  of  susceptibility  mapping  (SM)  which  we  refer  to  as  “Iterative  SWIM”.  The 
concept  is  to  take  the  first  guess  of  the  SM  and  keep  only  the  highest  components  and  then  do  a  forward 


6 


model.  Now  we  use  the  ill-defined  portions  of  k-space  from  this  forward  data  and  embed  it  into  the 
original  k-space  data  and  reconstruct.  This  has  a  dramatic  effect  in  terms  of  removing  artifacts.  An 
example  data  simulation  is  shown  below  (Figure  1).  Note  that  even  after  the  first  iteration  there  is  a 
dramatic  improvement  in  the  image. 


*i(r)  map 


a  =  0.03 

*i(r)  map 


a  =  0.1 

*i(r)  map 


a  =  0.2 

*i(r)  map 


a  =  0.3 
^Ti(r)  map 


5=  0 
rad 


Figure  1:  Simulations  showing  the  comparison  of  the  calculated  susceptibility  distributions  for  a  cylinder 
perpendicular  to  B0,  indicated  by  a  white  long  arrow,  at  different  threshold  values  (a)  applied  to  g(k)  as 
well  as  the  initial  *i(r)  map.  The  diameter  of  the  cylinder  is  32  pixels  with  a  susceptibility  Ay  of  0.45 
ppm  inside  the  cylinder,  a)  The  top  row  of  images  shows  calculations  with  no  added  phase  noise  and  the 
second  and  the  third  row  shows  calculations  with  added  phase  noises  5=  0.025  and  0.05  radians, 
respectively.  The  first  column  of  images  are  initial  Afi(r)  maps  for  reference,  b)  The  variation  of  the  mean 
calculated  susceptibility  inside  the  cylinder  with  different  threshold  value,  a.  The  mean  susceptibility 
value  is  independent  of  the  noise  level;  therefore,  only  mean  value  from  6=  0  is  provided.  The  error  bar 
shows  the  standard  deviation  (SD)  of  the  calculated  susceptibility  values,  c)  The  variation  of  the  root 
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mean  squared  error  (RMSE)  of  the  susceptibility  values  outside  the  cylinder  as  a  function  of  the  threshold 
value,  a,  and  the  noise  level. 

g)  Simultaneous  acquisition  of  both  MR  angiography  and  venography  data. 

We  have  developed  a  dual  echo  SWI  technique  that  will  provide  the  ability  to  not  only  perform 
SWI  and  SWIM  but  also  to  calculate  T2*.  This  is  not  the  ordinary  double  echo  sequence  but  a  fully  flow 
compensated  double  echo  sequence  so  no  interference  from  arteries  takes  place.  It  can  image  both  the 
arteries  and  veins  that  the  same  time,  which  could  potentially  assess  a  relative  complete  picture  of  vessel 
damage  after  brain  injury  (see  Figure  2).  Specifically,  we  have  designed  a  flow  rephase/dephase 
interleaved  double-echo  GRE  sequence  which  offers  fully  flow  rephased  SWI  data  for  venography,  and 
TOF  MRA  plus  flow  dephased  data  for  enhancing  angiography  using  a  non-linear  subtraction  (NLS) 
method.  To  our  knowledge,  currently  no  MR  method  other  than  the  interleaved  double-echo  sequence  we 
proposed  can  provide  all  the  data  needed  for  both  MRAV  imaging  and  selective  MRA  enhancement  in  a 
single  scan.  Three  different  types  of  image  contrast,  i.e.,  TOF  MRA,  SWI  and  flow  dephased  dark  blood 
images,  could  be  obtained  simultaneously  with  perfect  alignment  to  each  other  so  that  no  realignment 
process  is  necessary.  In  conjunction  with  NLS,  one  can  obtain  high  resolution  SWI  and  enhanced  TOF 
MRA  without  the  need  for  contrast  agent  (Figure  2  a-b),  and  the  process  can  be  easily  automated  without 
much  additional  computational  power.  NLS  offers  a  clean  angiography  free  of  venous  contamination  as  if 
acquired  with  maximal  TOF  effects.  Most  noticeably,  NLS  seems  to  even  outperform  the  low  dose  CE 
TOF  data  by  its  higher  artery-tissue  contrast  (Figure  2  c-e),  extra  revelation  of  some  small  arteries  (e.g.  at 
the  frontal  lobe),  uniform  background  signal  and  the  lack  of  veins  (except  for  the  easily  distinguishable 
superior  sagittal  sinus).  For  future  works,  the  acquisition  time  may  be  further  reduced  to  provide  whole 
brain  coverage  using  a  segmented  echo  planer  imaging  (seg-EPI)  readout  approachor  compress  sensing 
(CS)  since  the  subtracted  images  have  highly  sparse  information. 


Figure  2:  MIP  images  of  (a)  normal  TOF  MRA;  (b)  SWI;  (c)NLS  and  (d)  low  dose  CE  MRA  with  5cc 
gadopentetatedimeglumine.  (a-c)  were  collected  with  a  single  scan  of  the  interleaved  double-echo 
sequence  before  the  use  of  contrast  agent,  (e)  Normalized  profiles  of  one  representative  artery,  extracted 
from  the  correspondent  single  slice  images. 


h)  Investigate  the  early  response  of  the  neurovascular  system  to  brain  concussion 

In  our  preliminary  study,  we  investigated  the  neural  vascular  response  to  mild  TBI  at  the  acute  stage. 
Specifically,  we  adopted  a  combined  use  of  both  SWIM  and  ASL  to  answer  this  question. 

Subject  Enrollment:  Twelve  mTBI  patients  were  recruited  in  emergency  setting  in  our  Level- 1  trauma 
center.  Eighteen  aged-  and  gender-matched  healthy  controls  were  also  recruited  for  comparison.  All 
patients  met  the  definition  of  mTBI  by  the  American  Congress  of  Rehabilitation  Medicine  (ACRM)  with 
Glasgow  Coma  Scale  (GCS)  score  of  13-15  at  emergency  entry.  All  patients  completed  a  short 
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neurocognitive  test,  standard  assessment  of  concussion  (SAC),  at  the  acute  stage.  SAC  is  a  short 
instrument  to  measure  patients  orientation,  delayed  recall,  attention,  and  executive  function. 

SWIM  Processing:  All  SWI  images  were  processed  to  produce  SWIM  data  by  using  our  previously 
developed  method  in  our  in-house  developed  SPIN  software  [1],  Specifically,  the  following  steps  were 
used:  a)  High  pass  filtering  of  original  phase  image  to  get  SWI  phase  image,  b)  Inverse  transform  to 
generate  susceptibility  maps  of  the  veins,  c)  Regularization  and  interpolation  of  the  data  in  k-space  of  the 
phase  image  to  reduce  reconstruction  artifacts,  and  d)  finally,  minimum  intensity  projection  (mIP)  of 
slices  over  8  mm  thickness. 

After  SWIM  processing,  eight  major  veins  and  small  transmeduallar  veins  were  selected  by  using 
semi-automated  region  of  interest  (ROI)  analysis  to  quantify  blood  oxygen  saturation  of  each  one  (see 
Figure  3).  A  student  t-test  was  used  to  compare  the  control  and  patient  group  data. 

ASL  Data  Processing:  First,  the  skull-stripped  untagged  T2  weighted  image  was  used  to  create  a  mask  to 
skull-strip  the  relative  CBF  (rCBF)  images.  The  skull-stripped  untagged  T2  weighted  images  were  then 
normalized  to  the  T1  weighted  ICBM  template  by  using  SPM8  to  generate  a  transformation  matrix.  Next, 
the  transformation  matrix  was  applied  to  the  rCBF  images.  Then,  the  Wake  Forest  University  (WFU) 
Pickatlas,  which  is  in  the  same  space  as  ICBM  template,  was  used  to  extract  the  rCBF  values  from 
specific  predefined  structures  and  lobes.  rCBF  values  from  frontal,  parietal,  occipital  and  temporal  lobes 
were  extracted.  In  addition,  values  from  thalamus  and  basal  ganglia,  including  caudate  nucleus,  putamen 
and  globus  pallidus,  were  also  extracted. 


Figure  3:  Measurement  of  susceptibility  of  major  veins.  As  numbered  on  the  figure,  the  major  veins  are 
1-Septal  Vein  (right),  2-Septal  Vein  (left),  3-  Central  Septal  Vein,  4-Thalamic  Striatal  Vein  (right),  5- 
Thalamic  Striatal  Vein  (left),  6-Internal  Cerebral  Vein,  7-Basilar  Vein  of  Rosenthal  (right),  8-Basilar  Vein 
of  Rosenthal  (left),  9-Transmedullar  Vein  (right),  and  10-Transmedullar  Vein  (left). 
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Results: 


Our  data  demonstrate  that:  a)  patient  group  has  significant  lower  susceptibility  signal  than  control  group 
in  left  thalamus  striatal  veins  (See  Figure  4),  which  implies  that  patient  group  has  lower  deoxygenation  in 
thalamus  region  than  controls;  b)  Patient  group  has  significant  higher  CBF  levels  than  control  group  in 
left  frontal  lobe  and  left  striatal  region,  including  caudate,  putamen  and  globuspallidus  (See  Table  1);  c) 
patient  group  demonstrated  significantly  reduced  SAC  scores  than  published  normative  data,  especially  in 
delayed  recall  (See  Figure  5);  and  d)  the  imaging  data  is  not  correlated  with  SAC  data. 


Table  1.  rCBF  of  Major  Brain  Structures 


Structure 

Thal-L 

Thal-R 

Stri-L 

Stri-R 

Caud-L 

Caud-R 

Put-L 

Put-R 

GP-L 

GP-R 

Control 

2184.25 

2217.42 

2127.25 

2145.00 

2114.42 

2123.83 

2136.83 

2154.42 

2151.58 

2125.83 

Patients 

2236.86 

2272.43 

2190.50 

2201.00 

2172.67 

2182.00 

2188.50 

2216.00 

2190.33 

2131.67 

T-Test 

0.23 

0.28 

0.02 

0.14 

0.02 

0.09 

0.03 

0.14 

0.12 

0.96 

Susceptibility  of  Left  Thalamic  Striatal  Vein 


150 


100 


50 


0 


■Controls 
■  Acute  Patients 
■Patient  Follow  up 


Subject  Groups 

Figure  4:  Comparison  of  left  Thalamic  Striatal  vein  at  different  time  points.  Star  sign  indicates  statistical 
significance. 
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Figure  5:  Testing  scores  of  Delayed  Recall  in  the  Standard  Assessment  of  Concussion  (SAC).  As 
indicated  by  the  bar  chart,  patients  at  acute  stage  have  significantly  lower  scores  of  delayed  recall,  which 
suggests  of  memory  problem. 

i)  Investigating  medullary  vein  damage 

A  new  dimension  in  imaging  technology  is  the  discovery  of  vascular  damage  in  the  form  of  extravasated 
blood,  thrombus  and  reduced  levels  of  oxygenation  that  may  better  illuminate  our  understanding  of  brain 
injury,  particularly  related  to  trauma.  Using  the  SWI  and  SWIM  technology  discussed  above,  we  have 
recognized  that  there  are  major  changes  in  iron  content  and  local  oxygen  saturation,  important  for 
evaluating  potential  changes  in  local  perfusion  or  tissue  function  (similar  to  what  is  seen  in  stroke  using 
SWI).  SWIM  can  also  be  used  to  monitor  changes  in  iron  content  over  time  which  can  be  used  to  see  if 
previous  iron  deposition  is  being  resorbed  or  if  bleeding  continues,  both  important  diagnostic  pieces  of 
information  for  the  clinician. 

In  this  current  project,  we  have  demonstrated  that  there  is  a  lower  impact  load,  either  inertia  or  direct 
impact  forces,  that  may  damage  only  veins,  and  we  have  shown  medullary  vein  damage  that  hasn’t  been 
visualized  with  other  techniques  (Figure  6).  The  medullary  veins  drain  the  frontal  white  matter  of  the 
brain,  so  reduced  blood  flow  here  could  possibly  impair  the  higher  level  frontal  neurocognitive  function. 
Treatments  that  improve  blood  flow  to  the  brain,  then,  might  be  a  promising  direction  to  pursue. 
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Figure  6:  An  example  mild  traumatic  brain  injury  case  in  the  presence  of  both  blast  and  subsequent 
traumatic  brain  injury.  The  left  septal  vein  and  medullary  veins  draining  into  it  are  both  show  abnormal 
signal  representing  either  iron  deposition  or  increased  levels  of  deoxyhemoglobin  or  both.  Left  image: 

SWI  data.  Middle  image:  SWI  filtered  phase.  Right  image:  SWIM  data  showing  quantified  iron  content. 

j)  Investigating  medullary  vein  structure  in  three  dimensions 

Overall  results  from  modeling  of  major  vessels  suggests  that  blood  vessels  increase  the  stiffness  of  the 
brain  and  help  to  resist  both  rotational  and  linear  impulses  and  may  also  increase  the  threshold  of  neural 
injury.  However,  there  is  no  data  that  describe  the  mechanical  interaction  between  the  vessels  and  brain. 
Medullary  vein  damage  may  be  due  to  the  presence  of  internal  shear  forces  (intracerebral  stresses).  How 
this  occurs  likely  depends  on  the  stiffness  and  natural  frequency  of  the  structure  impacted.  For  impacts 
with  structures  of  relatively  soft  stiffness  the  brain  tends  to  move  with  the  skull  because  of  the 
predominantly  low  frequencies  involved  in  these  impacts.  These  movements  result  in  intracerebral  strains 
(ICS)  that  can  produce  lesions  in  the  central  parts  of  the  brain.  Knowing  the  tissue  properties  and  the  local 
structures  of  the  vessels  and  how  they  are  embedded  in  the  surrounding  tissue  is  necessary  for  proper 
modeling  of  the  problem.  The  SWI  and  MRA  data  that  we  have  collected  makes  it  possible  to  create  a 
detailed  model  of  the  medullary  veins  down  to  250  microns  in  size  (see  Figure  7).  Armed  with  this  model 
and  estimates  for  the  tissue  properties  it  may  be  possible  in  the  future  to  predict  local  venous  damage.  The 
determination  of  vessel  vulnerability  could  help  engineer  safety  devices  such  as  helmets  to  protect  the 
brain  from  forces  that  cause  venous  damage. 
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Figure  7:  Modeling  the  medullary  veins.  An  SWI  projection  image  is  used  to  find  the  medullary  veins 
and  then  they  are  extracted  in  3D.  a)  An  example  image  from  SWI  data  showing  the  septal  vein  and  the 
confluence  of  several  medullary  veins,  b)  A  projection  of  a  3D  model  of  the  medullary  veins  branching 
off  the  septal  vein.  The  resolution  in  the  model  is  62.5  microns  and  the  vessels  range  in  size  from  1mm 
(the  septal  vein)  to  the  first  branchings  (0.5mm)  to  the  second  branchings  (0.25mm).  These  veins  can  be 
damaged  in  frontal  or  side  impact. 
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KEY  RESEARCH  ACCOMPLISHMENTS 


The  following  accomplishments  have  been  made  in  the  past  year: 

•  Final  approval  of  study  protocol  by  three  Institutional  Review  Boards  (IRBs). 

•  Recruitment  and  follow  up  of  patients  and  controls  is  established  and  in  progress. 

•  Analysis  of  preliminary  data.  Our  data  demonstrated  an  increased  blood  flow  as  early  response  of 
mild  TBI  patients. 

•  A  new  finding  of  medullary  vein  damage  has  been  discovered  that  is  present  not  only  in  moderate 
and  severe  TBI  but  also  in  mild  TBI. 

•  Technical  development  of  susceptibility  weighted  imaging  and  mapping  (SWIM),  including: 

>  A  method  of  simultaneous  acquisition  of  both  MR  angiography  and  venography  has  been 
developed. 

>  An  advanced  iterative  approach  has  been  developed  to  improve  the  post  processing  of 
SWIM  data. 
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REPORTABLE  OUTCOMES 


Papers: 

a)  A  review  on  SWI  and  fMRI  was  published:  E.  M.  Haacke  and  Y.  Ye.  "The  role  of  susceptibility 
weighted  imaging  in  functional  MRI".  Neuroimage:  2012  (epub,  ahead  of  print). 

b)  DTI  paper  on  mild  TBI.  We  are  currently  working  on  a  paper  related  with  DTI  study  of  mild  TBI 
by  using  our  DTI  analysis  approach  as  mentioned  above. 

c)  A  book  chapter  entitled  “The  presence  of  venous  damage  and  microbleeds  in  traumatic  brain 
injury  and  the  potential  future  role  of  angiographic  and  perfusion  magnetic  resonance  imaging”  is 
in  press. 

d)  A  book  chapter  is  published:  Zhifeng  Kou,  Randall  R  Benson,  E  Mark  Haacke.  "Advanced 
Neuroimaging  of  Mild  Traumatic  Brain  Injury”,  in  “Biomarkers  for  Traumatic  Brain  Injury,” 
edited  by  Svetlana  Dambinova,  Ronald  L  Hayes,  Kevin  K.W.  Wang  and  David  E  Thurston. 
Royal  Society  of  Chemistry,  May,  2012. 

e)  A  paper  on  the  advancement  of  SWIM  has  been  accepted  in  MRM.  J.  Tang  et  al.  “Improving 
susceptibility  mapping  using  a  threshold-based  K-space/image  domain  iterative  reconstruction 
approach.”,  MRM,  Juen  26,  2012,  epub,  ahead  of  print. 

f)  A  manuscript  comparing  ASL  with  SWIM  data  is  in  progress. 

Presentations: 

Our  group  has  given  a  number  of  presentations  in  national  and  international  meetings. 

2011 

a)  The  project  PI,  Dr.  E  Mark  Haacke,  was  invited  by  the  Human  Brain  Mapping  society  to  give  a 
presentation  on  SWIM  on  June  29th,  2011,  Quebec  City,  Canada. 

b)  Dr.  Haacke  also  met  with  members  of  Senator  Harry  Reid’s  office  in  Las  Vegas  regarding  his 
research  on  August  30th,  201 1 

c)  Dr.  Zhifeng  Kou  gave  a  plenary  presentation  on  brain  trauma  in  2011  NABIS  meeting  on  September 
15th,  2011,  New  Orleans,  Louisiana. 

d)  Dr.  Haacke  organized  and  held  a  TBI  Workshop  on  November  18,  2011  from  8:00am  -  4:00pm  at 
Wayne  State  University.  There  were  four  topic  sessions  with  five  speakers  in  each  section  for  around 
twenty  short  talks.  Dr.  Haacke  spoke  about  the  new  venous  concepts  and  Dr.  Kou  about  progress  in  mild 
TBI. 

e)  Wayne  State  University  -  WSU  Program  in  Traumatic  Brain  Injury  Research  (PTBIR).  December  8th, 
2011. 

2012 

f)  Henry  Ford  Hospital.  The  role  of  the  vasculature  in  neurodegenerative  disease  with  a  focus  on  TBI  and 
CCSVI.  January  5th,  2012. 

g)  The  Society  for  Cardiovascular  Magnetic  Resonance.  CCSVI:  Abnormal  Venous  Flow  and 
Neurodegenerative  Disease.  February  5th,  2012. 

h)  SWI  fMRI.  American  Society  of  Functional  Neuroradiology  (ASFNR)  6th  Annual  Scientific  Meeting. 
March  7th,  2012.  Orlando,  FL,  USA. 

i)  SWI:  Advanced  Applications,  Oxygenation,  Phase  Quantitation.  American  Society  of  Functional 
Neuroradiology  (ASFNR)  6th  Annual  Scientific  Meeting.  March  8th,  2012.  Orlando,  FL,  USA. 

j)  Moderator  at  the  American  Society  of  Functional  Neuroradiology  (ASFNR)  6th  Annual  Scientific 
Meeting.  Case  Based  SWI  and  Advanced  Imaging.  March  8th,  2012.  Orlando,  FL,  USA. 
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k)  Susceptibility  weighted  imaging  at  the  19th  Annual  Audiology  Conference  at  the  VA  Hospital  in 
Johnson  City,  Tennessee,  on  June  21st,  2012. 

Grants,  Honors  and  Awards 

•  Dr.  Zhifeng  Kou  obtained  a  Seed  Grant  from  the  International  Society  for  Magnetic  Resonance  in 
Medicine  (ISMRM)  in  the  annual  conference  in  Melbourne,  Australia  in  May  2012. 

•  Dr.  Zhifeng  Kou  won  the  prestigious  Charles  Hayes  award  from  the  North  American  Brain  Injury 
Society  (NABIS)  in  2011. 
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CONCLUSION 


Our  work  during  this  first  year  has  already  led  to  some  major  breakthroughs  in  understanding  traumatic 
brain  injury.  First,  we  have  observed  medullary  vein  damage  even  in  mild  TBI.  The  implications  for  this 
are  two-fold:  we  may  be  able  to  understand  why  certain  individuals  are  suffering  from  frontal  white 
matter  neurocognitive  problems  and  we  may  be  able  to  propose  a  treatment  for  these  individuals.  Second, 
we  have  developed  a  power  quantitative  susceptibility  mapping  advance  called  iterative  SWIM  which 
opens  the  door  to  mapping  bleeding  in  TBI.  This  has  the  importance  of  monitoring  current  damage  and 
patient  recovery  in  the  sense  that  one  can  follow  tissue  recovery  and  clearing  of  iron  from  the  tissue. 
Third,  we  have  married  SWI  and  MRA  to  allow  for  a  complete  vascular  study  of  TBI.  The  impact  here  is 
that  we  should  be  able  to  better  understand  the  brain's  hemodynamics  post  trauma.  In  summary,  we  are  on 
track  with  our  technical  development  and  are  now  on  track  to  collecting  the  human  data.  Although  this 
may  take  some  time,  we  anticipate  being  able  to  complete  this  data  collection  by  the  end  of  year  two.  Our 
future  work  will  continue  along  the  lines  proposed  in  the  original  statement  of  work  where  we  will 
compare  the  imaging  at  different  time  points  and  look  for  biomarkers  of  patient  outcomes. 
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FULL  PAPER 


Magnetic  Resonance  in  Medicine  000:000-000  (2012) 


Improving  Susceptibility  Mapping  Using  a 
Threshold-Based  K-Space/hnage  Domain  Iterative 
Reconstruction  Approach 

J.  Tang,1  S.  Liu,1  J.  Neelavalli,2  Y.  C.  N.  Cheng,2  S.  Buch,1  and  E.  M.  Haacke1-3* 


To  improve  susceptibility  quantification,  a  threshold-based  k- 
space/image  domain  iterative  approach  that  uses  geometric 
information  from  the  susceptibility  map  itself  as  a  constraint 
to  overcome  the  ill-posed  nature  of  the  inverse  filter  is  intro¬ 
duced.  Simulations  were  used  to  study  the  accuracy  of  the 
method  and  its  robustness  in  the  presence  of  noise.  In  vivo 
data  were  processed  and  analyzed  using  this  method.  Both 
simulations  and  in  vivo  results  show  that  most  streaking  arti¬ 
facts  inside  the  susceptibility  map  caused  by  the  ill-defined 
inverse  filter  were  suppressed  by  the  iterative  approach.  In 
simulated  data,  the  bias  toward  lower  mean  susceptibility  val¬ 
ues  inside  vessels  has  been  shown  to  decrease  from  around 
10%  to  2%  when  choosing  an  appropriate  threshold  value  for 
the  proposed  iterative  method.  Typically,  three  iterations  are 
sufficient  for  this  approach  to  converge  and  this  process 
takes  less  than  30  s  to  process  a  512  x  512  x  256  dataset. 
This  iterative  method  improves  quantification  of  susceptibility 
inside  vessels  and  reduces  streaking  artifacts  throughout  the 
brain  for  data  collected  from  a  single-orientation  acquisition. 
This  approach  has  been  applied  to  vessels  alone  as  well  as  to 
vessels  and  other  structures  with  lower  susceptibility  to  gen¬ 
erate  whole  brain  susceptibility  maps  with  significantly 
reduced  streaking  artifacts.  Magn  Reson  Med  000:000-000, 
2012.  ©2012  Wiley  Periodicals,  Inc. 

Key  words:  oxygen  saturation;  susceptibility  mapping; 
susceptibility  weighted  imaging 

Susceptibility  weighted  imaging  (SWI)  using  phase  infor¬ 
mation  has  become  an  important  clinical  tool  (1-3). 
However,  the  use  of  phase  information  itself  has  stimu¬ 
lated  great  interest  both  as  a  source  of  contrast  (4-6)  and 
a  source  for  producing  susceptibility  maps  (SM)  (7-24). 
The  impetus  for  solving  the  inverse  problem  from  mag¬ 
netic  field  perturbation  came  from  the  work  described  by 
Deville  et  al.  (25).  This  was  noted  by  Marques  and  Bow- 
tell  in  2005  (26).  Salomir  et  al.  (27)  were  the  first  group 
to  utilize  this  concept  in  MRI.  Unfortunately,  this 
inverse  process  is  ill-posed  and  requires  a  regularization 
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procedure  to  estimate  the  SM.  There  are  a  variety  of 
approaches  to  tackle  this  problem  (7-24).  One  unique 
method  uses  a  multiple  orientation  data  acquisition  to 
remove  the  singularities  (17).  Constrained  regularizations 
(14,20,22,23)  have  shown  good  overall  results,  but  they 
require  longer  reconstruction  times  and  assumptions 
about  the  contrast  in  or  near  the  object  to  be  detected. 
Threshold-based,  single-orientation  regularization  meth¬ 
ods  (TBSO)  (11,15,18,24)  provide  the  least  acquisition 
time  and  the  shortest  computational  time  to  calculate 
SM.  However,  their  calculated  SMs  lead  to  underesti¬ 
mated  susceptibility  values  (Ay)  and  display  severe 
streaking  artifacts  especially  around  structures  with  large 
Ax,  such  as  veins  or  parts  of  the  basal  ganglia. 

Based  on  TBSO  approaches,  we  propose  an  iterative 
method  to  overcome  the  singularities  in  the  inverse  filter 
and  produce  improved  accuracy  for  susceptibility  map¬ 
ping.  In  this  approach,  we  iteratively  replace  k-space  val¬ 
ues  associated  with  the  SM,  x(k),  near  the  singularities 
to  obtain  an  almost  artifact  free  SM,  xM-  The  k-space 
values  used  for  substitutions  are  estimated  using  struc¬ 
tural  information  from  the  masked  version  of  xM-  Simu¬ 
lations  using  2D  cylinders  and  full  brain  3D  models 
were  performed  to  examine  the  efficacy  of  this  iterative 
approach.  High  resolution  human  data  are  also 
evaluated. 


MATERIALS  AND  METHODS 

Briefly,  the  expression  for  the  susceptibility  distribution 
(26,27)  derived  from  the  phase  data  can  be  written  as 
(for  a  right  handed  system  (28)): 


X(r)  =  FT- 


, ,  x  FT 

d,(r) 

[g(k) 

■ yB0iE\ 
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g(k) 


1  kz2 

3  kx2  +  k2  +  kz2 
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and  4>(r)  is  the  phase  distribution,  TE  is  the  echo  time,  7 
is  the  gyromagnetic  ratio  for  hydrogen  protons,  B0  is  the 
main  field  strength,  kx,  ky,  and  kz  are  coordinates  in  k- 
space,  and  g[k)  is  the  Green’s  function  or  filter.  Clearly, 
the  analytic  inverse  filter  g_1(k)  =  l/g(k),  is  ill-posed 
when  g[k)  is  equal  or  close  to  zero,  i.e.,  points  on  or  near 
two  conical  surfaces  in  k-space  at  the  magic  angles  of 
54.7°  and  125.3°  from  the  B0  axis.  This  ill-posedness 
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FIG.  1 .  Illustration  of  the  iterative  reconstruction  algorithm  to  obtain  artifact  free  yfr)  maps. 


leads  to  severe  artifacts  (including  severe  streaking)  in 
x(k)  and  noise  amplification  (29).  Thus,  for  a  proper  pixel- 
by-pixel  reconstruction  of  xW,  recovering  the  correct  val¬ 
ues  of  x(k)  near  the  region  of  singularities  is  critical. 

K-space  Iterative  Approach 

If  the  shapes  of  the  structures  of  interest  are  known,  then 
one  can  use  this  information  in  the  SM  to  create  a  more 
accurate  k-space  of  said  SM  in  the  conical  region.  The 
structure  of  the  vessels  is  obtained  directly  from  the  first 
pass  SM  Xi=oW-  The  detailed  steps  of  the  iterative 
method  are  discussed  below  and  shown  in  Fig.  1. 

•  Step-1:  An  initial  estimate  of  the  SM,  Xj=oW,  is 
obtained  by  applying  a  regularized  version  of  the 
threshold-based  inverse  filter,  gfgg{k)  (18),  in  Eq.  1 
using  the  suggested  threshold  value,  thr  =  0.1.  The 
subscript  “i”  denotes  the  SM  after  the  ith  iteration 
(“i  =  0”  denotes  the  initial  step  before  doing  the 
iterative  method  and  i  =  1  for  the  first  iteration). 

•  Step-2:  The  geometry  of  the  structures  of  interest  is 
extracted  from  Xj=oW  using  a  binary  vessel  mask, 
i.e.  outside  the  veins,  the  signal  in  the  mask  is  set  to 
zero,  and  inside  it  is  set  to  unity.  Since  streaking 
artifacts  associated  with  veins  in  the  SM  are  usually 
outside  the  vessels,  after  multiplying  the  x;W  map 
by  the  mask,  little  streaking  remains  in  the  SM.  This 
leads  to  Xvm.iW  as  shown  in  Fig.  lb.  Vessel  mask 
generation  will  be  addressed  in  the  next  section. 

•  Step-3:  Xvm.iW  is  obtained  by  Fourier  transforma¬ 
tion  Of  Xvm.iW  (Fig.  lc). 

•  Step-4:  The  predefined  ill-posed  region  of  k-space 
in  Xvm.iW  is  extracted  (Fig.  Id).  These  extracted  k- 
space  data  are  denoted  by  Xvm.cone.iW-  The  size  of 
Xvm.cone.iW  is  decided  by  a  threshold  value,  a, 
which  is  assigned  to  g(k).  For  the  matrix  size  512  x 
512  x  512,  the  percentages  of  the  cone  region  in 


k-space  for  a  given  a,  are  2.4%  (a  =  0.01),  24.1%  (a 
=  0.1),  47.1%  (a  =  0.2),  and  70.6%  (a  =  0.3),  respec¬ 
tively.  When  a  increases,  the  size  of  Xvm.cone.iW 
increases  too.  If  a  is  increased  too  much  then  most 
of  the  original  information  will  be  lost. 

•  Step-5:  Data  from  Xvm.cone.iW  and  Xi=o W  (Fig.  le) 
are  merged.  This  means  part  of  Xi=oW  has  been 
replaced  by  Xvm.cone.iW-  The  merged  data  are 
denoted  by  Xmerged.iW  (Fig.  If). 

•  Step-6:  Inverse  Fourier  transformation  of  XmergediW 
gives  the  improved  SM,  Xi+iW  (Fig.  lg). 

•  Step-7:  XiW  in  step-1  is  replaced  by  Xi+iW  from 
step-6  and  the  algorithm  is  repeated  until 

[(x;-W-xmW)2]/iv  <  a  [3] 

where  N  is  the  number  of  pixels  in  x;W  and  £  is  the 
tolerance  value  (chosen  here  to  be  0.004  ppm). 

Binary  Vessel  Mask  Generation 

The  binary  vessel  mask  was  generated  using  threshold¬ 
ing  from  the  xW  map  itself.  The  detailed  steps  are  dis¬ 
cussed  below  and  shown  in  Fig.  2. 

•  Step-1:  A  threshold,  tha,  is  applied  to  Xr=oW  to  cre¬ 
ate  an  initial  binary  vessel  mask,  M0.  The  pixels 
whose  susceptibility  values  are  lower  than  thi  will 
be  set  to  zero  while  those  greater  than  or  equal  to 
tha  will  be  set  to  unity.  In  this  study,  a  relatively 
low  susceptibility  of  0.07  ppm  is  used  for  tht  to  cap¬ 
ture  most  vessels.  However,  this  choice  of  threshold 
inevitably  includes  other  brain  structures  in  M0  that 
have  high  susceptibility. 

•  Step-2:  A  morphological  closing  operation  is  per¬ 
formed  to  fill  in  holes  in  M0  to  generate  an  updated 
mask  Mi. 
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FIG.  2.  Illustration  of  the  binary  vessel  mask  generation  process. 


•  Step-3:  A  median  filter  is  applied  to  remove  noise 
in  Ma  and  create  M2. 

•  Step-4:  False  positive  data  points  from  M2  are 
removed  as  follows:  First,  the  xM  map  is  Mipped 
over  five  slices  centered  about  the  slice  of  interest  to 
better  obtain  contiguous  vessel  information,  as  seen 
in  XmipW-  Second,  another  threshold,  th2  =  0.25 
ppm,  is  performed  on  XmipM  to  create  a  new 
Xmip_viiiM  and  binary  mask  MP,  which  only  contains 
predominantly  vessels.  Here,  0.25  ppm  was  chosen 
to  isolate  the  major  vessels  in  the  MIP  image.  Third, 
each  slice  from  M2  is  compared  with  MP  on  a  pixel- 
by-pixel  basis  to  create  M3.  If  a  data  point  from  M2 
does  not  appear  on  MP,  this  data  point  will  be 
treated  as  a  false  positive  and  removed  from  M2,  oth¬ 
erwise  this  point  is  retained.  This  process  can  be 
equally  well  applied  to  extract  other  tissues  by 
choosing  appropriate  values  for  tha  and  th2. 

2D  Cylinder  Simulations 

Simulation  of  a  two  dimensional  cylinder  and  its 
induced  phase  was  first  performed  using  a  8192  x  8192 
matrix.  A  lower  resolution  complex  image  was  then 
obtained  by  taking  the  Fourier  transform  of  this  matrix 
and  applying  an  inverse  Fourier  transform  of  the  central 
512  x  512  matrix  in  k-space.  This  procedure  is  to  simu¬ 
late  Gibbs  ringing  effects  caused  by  finite  sampling 
which  we  usually  see  in  conventionally  required  MR 
data  sets.  Gibbs  ringing  comes  from  discontinuities  in 
both  the  magnitude  and  phase  images.  To  avoid  Gibbs 
ringing  from  magnitude  discontinuities,  we  used  a  mag¬ 
nitude  image  with  a  uniform  signal  of  unity.  Cylinders 
with  diameters  32,  64,  128,  256,  512,  and  1024  were 
simulated  on  8192  x  8192  matrices  and  their  effective 


diameters  were  2,  4,  8,  16,  32,  and  64  on  512  x  512  mat¬ 
rices.  All  phase  simulations  were  performed  using  a  for¬ 
ward  method  (8,26,27,30)  with  B0  =  3  T,  Ay  =  0.45  ppm 
in  SI  units,  TE  =  5  ms,  and  the  cylinder  perpendicular  to 
the  main  magnetic  field.  The  susceptibility  value  of  0.45 
ppm  represents  venous  blood  when  the  hematocrit  (Hct) 
=  0.44,  Axdo  =  4tt  0.27  ppm  (31)  and  the  oxygen  satura¬ 
tion  level  =  70%,  where  Axdo  is  the  susceptibility  differ¬ 
ence  between  fully  deoxygenated  and  fully  oxygenated 
blood  (32).  A  relatively  short  echo  time  was  chosen  to 
avoid  phase  aliasing  that  can  affect  the  estimated  suscep¬ 
tibility  values. 

Selection  of  a  TBSO  Method  to  Generate  the  Xi=oM  Map 

TBSO  methods  (11,15,18,24)  use  a  truncated  g[k)  to  solve 
the  singularity  problem  in  the  inverse  filter  g_1(k)  when 
g(k)  is  less  than  a  predetermined  threshold  value,  thr. 
When  g{k)  <  thr,  g_1(k)  is  either  set  to  zero  (11,24);  or  to 
1/thr  (15);  or  set  to  g_1(k)  =  1/thr  first  and  then  g^1(k)  is 
brought  smoothly  to  zero  as  k  approaches  kzo.  This 
smoothing  is  accomplished  by  multiplying  g_1(k)  by 
a2(kz)  with  a(kz)  =  (kz  -  kzo)/  I  kzthr  -  kzo  I  where  kz  is 
the  z  component  of  that  particular  point  in  k-space,  kzo 
is  the  point  at  which  the  function  g_1(k)  becomes  unde¬ 
fined,  and  kzthr  is  the  kz  coordinate  value  where  I  g(k)  I 
=  thr  (18). 

SMs  using  the  methods  in  Refs.  11,15,18  were  calcu¬ 
lated  based  on  Eq.  1  using  the  2D  cylindrical  model. 
Equation  1  can  be  used  to  calculate  the  SM  for  the  simu¬ 
lated  2D  cylinder  model  perpendicular  to  the  main 
field  since  the  2D  perpendicular  model  is  a  special  case 
of  the  3D  model  with  1  slice  [10].  Streaking  artifacts  are 
obvious  in  all  three  SMs  (figures  are  not  shown).  The 
calculated  mean  susceptibility  values  inside  the  cylinder 
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FIG.  3.  a:  A  transverse  view  of  the  3D  brain  model,  b:  The  simulated  phase  map  from  the  model  using  parameters:  B0  =  3  T  and  TE  = 
18  ms  which  are  consistent  with  imaging  parameters  in  the  real  data  (c).  Images  (b)  and  (c)  have  the  same  window  level  setting. 


are  around  0.40  ±  O.Olppm  for  all  SMs.  The  background 
noise  levels,  (i.e.,  standard  deviation  of  the  susceptibility 
values)  measured  from  a  region  outside  the  streaking  ar¬ 
tifact  in  SM  using  Ref.  18  are  around  1/2  to  2/3  of  the 
background  noise  levels  in  SMs  using  Refs.  11,15  with 
thr  =  0.06,  0.07,  and  0.1,  which  are  the  optimal  thresh¬ 
old  values  suggested  in  Refs.  11,12,18.  Given  this  result, 
the  method  in  (18)  was  chosen  to  generate  a  Xi=oM  map. 

Finding  an  Optimal  Threshold  Value 

To  find  the  optimal  threshold,  a  series  of  xM  maps  were 
reconstructed  by  the  iterative  method  using  threshold 
values  a  of  0.01,  0.03,  0.07,  0.1,  0.15,  0.2,  0.25,  and  0.3. 
The  larger  this  threshold,  the  closer  the  final  estimate  for 
x(r)  will  be  to  XvmW-  The  optimal  threshold  value  was 
found  by  comparing  the  accuracy  of  the  estimated  sus¬ 
ceptibility  values  as  well  as  the  effects  on  reducing 
streaking  artifacts  in  the  reconstructed  xM  maps.  To 
study  the  effect  of  noise  in  xM  maps  due  to  the  noise  in 
phase  images,  complex  datasets  for  cylinders  of  diameter 
2,  4,  8,  32  voxels,  respectively,  were  simulated  with 
Gaussian  noise  added  to  both  real  and  imaginary  chan¬ 
nels.  Noise  was  added  in  the  complex  images  to  simulate 
a  SNRmagnitude  of  40:1,  20:1,  10:1,  and  5:1  in  the  magni¬ 
tude  images.  Since  o-p^ase  =  l/SNRmagnitude,  this  corre¬ 
sponds  to  crphase  =  0.025,  0.05,  0.1,  and  0.2  radian. 

To  estimate  the  improvement  in  the  SM  by  the  itera¬ 
tive  method,  we  used  a  root  mean  squared  error  (RMSE) 
to  measure  streaking  artifacts  outside  the  cylinder.  Back¬ 
ground  noise  in  the  SM  is  measured  in  a  region  away 
from  all  major  sources  of  streaking  artifacts  to  compare 
the  noise  measured  in  the  phase  image  (i.e.,  so  we  can 
correlate  noise  in  the  phase  with  the  expected  noise 
enhancement  from  the  inversion  process). 

Effect  of  High-Pass  Filter 

The  effect  of  high-pass  (HP)  filtering  the  phase  data  on 
the  xM  map  generated  by  the  iterative  method  was  also 
studied.  Phase  images  of  a  cylindrical  geometry  with 
diameters  of  2,  4,  8,  16,  32,  and  64  voxels  were  simu¬ 
lated.  Homodyne  HP  filters  (33)  with  a  2D  hanning  filter 
(full  width  at  half-maximum,  FWHM  =  4,  8,  16,  and  32 
pixels)  were  applied  on  these  phase  images  in  both  in¬ 


plane  directions.  SM  reconstructions  were  stopped  based 
on  the  criteria  in  step  7  of  the  iterative  process. 

Three  Dimensional  Brain  Model  Simulations 

To  address  the  potential  of  the  iterative  technique  to 
improve  the  SM  of  general  structures  such  as  the  basal 
ganglia,  a  3D  model  of  the  brain  was  created  including 
the:  red  nucleus  (RN),  substantia  nigra  (SN),  crus  cerebri 
(CC),  thalamus  (TH),  caudate  nucleus  (CN),  putamen 
(PUT),  globus  pallidus  (GP),  gray  matter  (GM),  white 
matter  (WM),  cerebrospinal  fluid  (CSF),  and  the  major 
vessels  (34).  The  structures  in  the  3D  brain  model  were 
extracted  from  two  human  3D  TV  weighted  and  T2 
weighted  data  sets.  Basal  ganglia  and  vessels  are  from 
one  person;  GM  and  WM  are  from  the  other  person’s 
data  set.  Since  all  structures  are  from  in  vivo  human 
data  sets,  this  brain  model  represents  realistic  shapes 
and  positions  of  the  structures  in  the  brain.  Susceptibil¬ 
ity  values  in  parts  per  million  (ppm)  for  the  structures 
SN,  RN,  PUT,  and  GP,  were  taken  from  Ref.  12  and 
others  were  from  measuring  the  mean  susceptibility 
value  in  a  particular  region  from  SMs  using  Ref.  18  from 
in  vivo  human  data:  RN  =  0.13,  SN  =  0.16,  CC  =  —0.03, 
TH  =  0.01,  CN  =  0.06,  PUT  =  0.09,  GP  =  0.18,  vessels  = 
0.45,  GM  =  0.02,  CSF  =  -0.014,  and  WM=0.  All  struc¬ 
tures  were  set  inside  a  512  x  512  x  256  matrix  of  zeros. 
The  phase  of  the  3D  brain  model  was  created  by  apply¬ 
ing  the  forward  method  (8,26,27,30)  to  the  3D  brain 
model  with  different  susceptibility  distributions  using 
the  imaging  parameters:  TE  =  5  ms  and  B0  =  3  T.  A  com¬ 
parison  between  the  phase  maps  from  this  brain  model 
and  a  real  data  set  is  shown  in  Fig.  3.  To  match  the 
imaging  parameters  of  the  real  data  set,  B0  =  3  T  and  TE 
=  18  ms  were  applied  for  the  results  presented  in  Fig.  3. 
Except  for  Fig.  3,  all  other  figures  in  the  paper  associated 
with  the  3D  brain  were  simulated  by  using  TE  =  5  ms. 

In  Vivo  MR  Data  Collection  and  Processing 

A  standard  high-resolution  3D  gradient  echo  SWI 
sequence  was  used  for  data  acquisition.  A  transverse  0.5 
mm  isotropic  resolution  brain  dataset  was  collected  at  3 
T  from  a  23-year-old  healthy  volunteer.  The  sequence  pa¬ 
rameters  were:  TR  =  26  ms,  flip  angle  =  15°,  read 
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FIG.  4.  Simulations  showing  the  comparison  of  the  calculated  susceptibility  distributions  for  a  cylinder  perpendicular  to  B0  at  different 
threshold  values  (a)  applied  to  g(k)  as  well  as  the  initial  x/  =  oW  map.  The  direction  of  B0  is  indicated  by  a  black  long  arrow.  The  suscep¬ 
tibility,  Ax.  inside  the  cylinder  is  0.45  ppm.  a:  The  comparison  of  the  converged  xi-bir)  map  with  the  x/=oM  map  for  a  diameter  of  32- 
pixel  cylinder,  where  b  is  the  iterative  step  required  to  reach  convergence.  In  this  data,  b  =  2  when  a  =  0.03,  b  =  3  when  a  =  0.1  and  b 
=  4  when  a  =  0.2  when  <rPhase  =  0.  The  top  row  of  images  shows  simulations  with  no  phase  noise.  The  second  and  the  third  row  show 
simulations  with  added  phase  noises  <jphase  =  0.025  and  0.05  radian,  respectively.  The  first  column  of  images  show  initial  xj-o(r)  maps 
for  reference,  b:  The  variation  of  the  mean  calculated  susceptibility  inside  the  cylinder  with  different  threshold  value,  a,  for  diameter  (d) 
=  2,  4,  8,  and  32  pixels  cylinders.  The  mean  susceptibility  value  is  independent  of  the  noise  level;  therefore,  only  mean  values  from 
'’'phase  =  0  were  provided,  c:  The  variation  of  the  RMSE  of  the  susceptibility  values  outside  the  cylinder  as  a  function  of  the  threshold 
value,  a,  and  the  noise  level.  The  d  =  32  pixels  cylinder  was  used  to  generate  (c).  The  range  of  the  gray-scale  bars  is  chosen  to  high¬ 
light  the  artifacts  in  the  images.  It  does  not  reflect  the  quantified  higher  susceptibility  values  inside  cylinders. 


bandwidth  =  121  Hz/pixel,  TE  =  14.3  ms,  192  slices,  and  a 
matrix  size  of  512  x  368.  To  reconstruct  X;=oM  with  mini¬ 
mal  artifacts,  the  following  steps  were  carried  out: 

1.  The  unwanted  background  phase  variations  were 
removed  using  either:  (a)  a  homodyne  HP  filter 
(FWHM  =  16  pixels)  (33)  or  (b)  Prelude  in  FMRIB 
Software  Library  (FSL)  (35)  to  unwrap  the  phase,  fol¬ 
lowed  by  the  process  of  Sophisticated  Harmonic  Arti¬ 
fact  Reduction  for  Phase  data  (SHARP)  (36)  with  a 
filter  radius  of  6  pixels.  To  reduce  artifacts  in  the  cal¬ 
culated  SMs,  regions  with  the  highest  phase  devia¬ 
tions  due  to  air/tissue  interfaces  were  removed 
manually  from  the  HP  filtered  phase  images  and  the 
phase  in  those  regions  were  set  to  zero. 

2.  A  complex  threshold  approach  (37)  was  used  to 
separate  the  brain  from  the  skull. 

3.  The  phase  image  with  an  original  matrix  size  of  512 
x  368  x  192  was  zero  filled  to  512  x  512  x  256  to 


increase  the  field-of-view  and  to  avoid  streaking 
artifacts  caused  by  the  edge  of  brain  to  alias  back  to 
the  reconstructed  SM. 

4.  The  regularized  inverse  filter,  g^glk)  (18)  was 
applied  to  obtain  X;=oM,  followed  by  the  iterative 
process  using  a  =  0.1.  For  in  vivo  data,  the  iterative 
program  was  terminated  at  the  third  iterative  step. 

RESULTS 

Selection  of  Threshold  Level  Based  on  Simulations 

To  find  the  optimal  threshold  value,  SMs  were  recon¬ 
structed  using  a  =  0.01,  0.03,  0.07,  0.1,  0.15,  0.2,  0.25, 
and  0.3,  respectively,  with  different  noise  levels  (Fig.  4). 
The  streaking  artifacts  shown  in  X;=oM  (the  first  column 
in  Fig.  4a)  have  been  significantly  reduced  by  the  itera¬ 
tive  method  and  fall  below  the  noise  level  when  a  >= 
0.1.  Also,  when  a  >=  0.1,  the  mean  susceptibility  value 
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inside  the  cylinder  was  found  to  increase  to  0.44  ppm 
when  the  diameter  of  the  cylinder  was  larger  than  8  pix¬ 
els  (Fig.  4b)  and  this  trend  is  independent  of  the  object 
size  and  the  noise  in  the  phase  image.  The  optimal  result 
in  terms  of  obtaining  the  true  susceptibility  value  was 
with  a  threshold  of  0.1.  Figure  4c  shows  a  plot  of  RMSE 
of  the  susceptibility  values  from  the  whole  region  out¬ 
side  the  32-pixel  cylinder  using  different  a.  The  RMSE 
of  the  susceptibility  values  decreases  as  a  increases. 
Therefore,  for  vessels,  a  value  of  a  =  0.3  would  be  the 
optimal  value.  Flowever,  a  large  threshold  value  means 
replacing  more  original  k-space  with  the  k-space  only 
consisting  of  vessel  information  which  will  reduce  the 
signals  from  other  brain  structures  and  blur  these  struc¬ 
tures.  Since  the  SM  using  a  =  0.1  already  reveals  the 
optimal  susceptibility  value  for  the  vessels  and  an  ac¬ 
ceptable  RMSE,  it  is  appropriate  to  choose  0.1  for  more 
general  applications  to  study  the  entire  brain. 

Figure  4a  compares  the  converged  X/=z>(r)  map  with  the 
Xi=oM  map,  where  b  is  the  iterative  step  required  to 
reach  convergence.  In  this  data,  b  =  2  when  a  =  0.01 
and  0.03,  b  =  3  when  a  =  0.07,  0.1,  and  0.15  and  b  =  4 
when  a  =  0.2,  0.25,  and  0.3  for  Up^ase  =  0-  When  Uphase 
increases,  more  iterative  steps  were  required  to  reach 
convergence.  For  instance,  the  maximum  iterative  step 
number  is  9  when  uphase  =  0-2  radians.  Using  a  noise 
level  of  0.025  radian  in  the  phase  image  as  an  example, 
greg(k)  (18)  leads  to  a  susceptibility  noise  of  roughly 
0.025  ppm  in  the  Xj=oM  map.  The  iterative  approach 
leads  to  a  slight  decrease  in  background  noise,  0.021 
ppm,  in  X;=3(r)  map  when  a  =  0.1.  The  background  noise 
was  measured  in  a  region  outside  the  streaking  artifact 
indicated  by  the  black  circle  in  Fig.  4a.  The  overall 
decrease  in  RMSE  in  the  background  (Fig.  4c)  is  a  conse¬ 
quence  of  both  a  decrease  in  streaking  artifacts  and  a 
reduction  in  thermal  noise  contribution. 

Selection  of  the  Optimal  Iterative  Step 

The  inverse  process  (18)  was  applied  to  the  dipole  field 
in  Fig.  5a  to  give  the  X;=oM  map  shown  in  Fig.  5b; 
prominent  streaking  artifacts  are  evident  in  this  image. 
Streaking  artifacts  are  significantly  reduced  at  each  step 
of  the  iterative  method  quickly  reaching  convergence 
(Fig.  5c-e).  The  largest  improvement  is  seen  in  the  first 
iterative  step,  which  is  verified  by  Fig.  5f,  showing  the 
difference  between  Fig.  5c  (x;=iM  map)  and  Fig.  5b 
(Xj=oM  map).  After  the  second  iteration,  we  can  see 
some  minor  streaking  reductions  (Fig.  5g,  the  difference 
between  the  Xi=iM  map  and  Xi-'zM  map).  The  mean  sus¬ 
ceptibility  value  approaches  0.44  ppm  in  a  single  step. 
Similar  results  (not  shown)  are  also  obtained  when  the 
iterative  method  is  run  with  different  aspect  ratios 
between  the  in-plane  resolution  and  the  through  plane 
resolution  (such  as  1:2  and  1:4).  The  iterative  results 
always  lead  to  higher  final  susceptibility  values  com¬ 
pared  to  the  initial  value  in  Xi=oM-  Finally,  even  when 
an  HP  filter  is  applied,  up  to  a  10%  increase  in  the  sus¬ 
ceptibility  is  realized  (Fig.  5i).  The  SMs  of  large  vessels 
benefit  from  a  low  order  HP  filter  (FWMH  =  4  pixels) 
and  small  vessels  up  to  8  pixels  benefit  from  a  HP  filter 
(FWMH  =  16  pixels). 


Effect  of  the  Iterative  Approach  on  Surrounding  Brain 
Tissues  in  the  3D  Brain  Model 

SM  Reconstruction  Using  a  Vessel  Mask  Only 

Figure  6a, d  represents  X;=oM>  without  noise  and  with 
0.025  radians  of  noise  in  phase  images.  Figure  6f  is  the 
vessel  map.  Streaking  artifacts  (delineated  by  the  black 
arrows)  are  obvious  in  Fig.  6a, d  and  significantly 
reduced  in  the  Xr=3W  maps  (Fig.  6b, e)  using  a  =  0.1.  Fig¬ 
ure  6c  is  the  Xr=3(r)  map  using  a  =  0.2.  As  can  be  seen, 
when  a  increases,  the  iterative  method  still  works  for 
vessels,  but  brain  tissues  become  more  blurred.  Figure  7a 
plots  the  mean  susceptibility  values  inside  the  vessel 
(vein  of  Galen),  GP,  SN,  RN,  PUT,  and  CN  from  X;=3(r) 
maps  generated  by  using  a  =  0.1,  0.15,  0.2,  0.25,  and  0.3, 
respectively.  The  susceptibility  value  in  the  brain  model 
and  Xi=oM  map  are  also  provided  in  the  plot  as  referen¬ 
ces.  Generally,  the  susceptibility  values  of  brain  tissues 
except  vessels  decrease  as  a  increases  while,  for  vessels, 
the  susceptibility  value  is  0.41  ppm  in  the  Xj=oM  map 
and  is  increased  to  0.45  ppm  in  the  Xi=3(r)  maps. 

SM  Reconstruction  Using  a  Mask  Including  Vessels  and 
Brain  Structures 

The  iterative  method  is  not  limited  to  improving  SM 
from  just  vessels;  it  can  also  be  applied  to  the  entire 
brain.  Figure  6g  shows  a  coronal  view  of  the  Xi=oM  map 
for  the  brain  model.  The  Xi=3W  map  using  a  mask  keep¬ 
ing  all  major  structures  (GP,  SN,  RN,  PUT,  CN)  and  ves¬ 
sels  is  shown  in  Fig.  6h.  In  practice,  this  is  equivalent  to 
setting  thresholds  in  the  Xj=oM  map  to  be  greater  than 
0.09  ppm  to  extract  all  these  high  susceptibility  struc¬ 
tures  from  the  Xj=oM  map  to  create  the  mask.  Figure  6h 
reveals  that  streaking  artifacts  associated  with  veins  as 
well  as  all  major  structures  have  been  reduced.  Figure  6i 
shows  the  difference  between  Fig.  6g,  h.  In  addition, 
streaking  artifacts  sometimes  cause  the  appearance  of 
“false”  structures.  For  instance,  there  is  no  internal  cap¬ 
sule  (IC)  included  in  the  model  (Fig.  61),  yet  we  see  an 
IC  like  structure  in  the  X;=oM  map  (Fig.  6j)  (indicated  by 
a  dashed  white  arrow  in  Fig.  6j).  The  iterative  method 
removes  the  streaking  artifacts  and  the  “false”  IC  (Fig. 
6k).  Figure  7b  shows  susceptibility  values  in  each  struc¬ 
ture  in  the  brain  model  for  Xi=oM  and  Xi=3(r)  when  the 
mask  includes  vessels  and  all  major  structures.  The 
underestimated  susceptibility  values  of  all  major  struc¬ 
tures  and  vessels  in  the  X;=oM  map  have  been  recovered 
by  the  iterative  method  in  the  X;=3W  map. 

Effect  of  Errors  in  the  Vessel  Map 

Accurately  extracting  vessels  from  Xi=oM  is  critical  for 
the  iterative  method.  Figure  8b-d  and  the  corresponding 
enlarged  views  (Fig.  8f-h)  from  the  rectangular  region 
indicated  in  Fig.  8a  show  the  Xr=3(r)  maps  using  an  accu¬ 
rate  (Fig.  8j),  a  dilated  (Fig.  8k),  and  an  eroded  (Fig.  81) 
vessel  map  to  show  the  effect  of  errors  in  the  vessel 
mask  on  the  Xi=  sM  map.  The  dilated  and  eroded  vessel 
maps  were  generated  using  Matlab  functions  based  on  a 
3-by-3  square  structuring  element  object.  The  susceptibil¬ 
ity  values  measured  from  a  vein  indicated  by  the 
white  arrow  in  Fig.  8e  are  0.40  ±  0.03  ppm  (Fig.  8e), 


Susceptibility  Mapping  Using  Iterative  Approach 


7 


FIG.  5.  a:  Phase  images  from  a  cylinder  with  a  diameter  of  32  pixels  are  simulated  with:  Ax  =  0.45  ppm,  S0  =  3  T,  and  TE  =  5  ms.  The 
cylinder  is  perpendicular  to  the  main  field.  No  thermal  noise  was  added  in  these  images,  b:  The  initial  xi-o(r)  map.  c:  The  SM  from  the 
first  iteration,  x/=iM  map,  (d)  x/=zW  map  and  (e)  x/=3(r)  maP  using  threshold  value  a  =  0.1.  The  SM  has  converged  at  x/= a(f)  map.  The 
streaking  artifacts  are  reduced  as  the  number  of  iterative  steps  increases,  f:  The  difference  image  of  (c)  subtracted  from  (b)  illustrates 
that  the  streaking  artifacts  were  reduced  by  the  iterative  procedure  and  the  largest  improvement  happens  in  this  first  iterative  step,  g: 
The  difference  image  of  the  x/=iM  map  subtracted  from  the  xt=z(r)  map  indicates  that  the  streaking  artifacts  were  further  reduced  by  the 
second  iterative  step,  h:  The  difference  image  of  x,= 2W  map  subtracted  from  x/= alf)  map  shows  much  less  improvement  at  the  third  iter¬ 
ative  step.  Thus  it  indicates  a  convergence  of  the  iterative  procedure.  All  images  were  set  to  the  same  window  level  setting  for  direct 
comparisons  and  for  enhancing  the  presence  of  the  streaking  and  the  remnant  error,  i:  The  effect  of  the  iterative  approach  on  the 
changes  in  susceptibility  values  from  HP  filtered  phase  images.  Differences  between  the  values  in  iterative  and  non-iterative  susceptibil¬ 
ity  map  reconstruction  (i.e.,  XconvergedW  -  x;-oM)  from  HP  filtered  phase  images  are  plotted  for  different  filter  sizes.  Results  for  four  filter 
sizes  (FWHM  =  4,  8,  16,  and  32  pixels)  are  shown  here.  Applying  an  HP  filter  leads  to  an  underestimation  of  Ax  (18).  The  iterative 
approach  helps  to  improve  the  accuracy  of  the  estimated  susceptibility  values.  The  range  of  the  gray-scale  bars  is  chosen  to  highlight 
the  artifacts  in  the  images.  It  does  not  reflect  the  quantified  higher  susceptibility  values  inside  cylinders. 


0.45  ±  0.03  ppm  (Fig.  8f),  0.45  ±  0.03  ppm  (Fig.  8g),  and 
0.40  ±  0.07  ppm  (Fig.  8h),  respectively.  The  iterative 
method  still  works  if  the  vessel  is  slightly  enlarged  but 
does  little  to  change  the  original  Xr=oM  map  if  the  ves¬ 
sels  are  too  small  or  absent  in  the  mask.  As  we  just  dis¬ 
cussed,  streaking  artifacts  produced  “false”  vessels  indi¬ 
cated  by  the  dashed  black  arrow  in  Fig.  8e  since  these 
vessels  are  not  in  the  model  (Fig.  8i).  These  false  vessels 
disappeared  in  Fig.  8f. 

Results  from  the  In  Vivo  Dataset 

In  the  in  vivo  example,  we  compare  the  differences 
between  SHARP  (Fig.  9a-d)  and  a  homodyne  HP  filter 
(FWHM  =  16  pixels)  (Fig.  9e-h).  Compared  to  the  trans¬ 
verse  view,  streaking  artifacts  are  more  obvious  in  the 
sagittal  or  coronal  view.  Figure  9a  shows  the  Xi=oM  map 
with  severe  streaking  artifacts.  The  streaking  artifacts 


were  significantly  reduced  in  the  Xi=3W  map  (Fig-  9b) 
using  a  =  0.1.  The  streaking  artifacts  associated  with  the 
superior  sagittal  sinus  vein  (indicated  by  two  black 
arrows  in  Fig.  9a)  were  significantly  decreased  in  Fig. 
9b, d.  The  subtracted  image  (Fig.  9c),  Fig.  9a  minus  Fig. 
9b,  reveals  the  removed  streaking  artifacts.  These  streak¬ 
ing  artifacts  are  one  of  the  reasons  why  the  Xi=oM  maps 
appear  noisy.  In  the  Xj=3M  map,  the  reduction  in  streak¬ 
ing  artifacts  from  individual  veins  leads  to  a  decrease  of 
noise  therefore  an  increased  SNR  of  veins.  If  veins  are 
the  only  interest,  even  a  threshold  of  0.2  can  work  rea¬ 
sonably  well  (Fig.  9d).  Two  relatively  big  veins,  Vl  and 
V2,  indicated  by  the  white  dashed  and  white  solid 
arrows,  respectively,  in  Fig.  9b,  were  chosen  to  measure 
the  susceptibility  values.  Results  are  provided  in  Table 
1.  The  susceptibility  values  of  these  two  veins  have  been 
improved  by  roughly  16%  by  the  iterative  method.  The 
standard  deviation  of  the  susceptibility  values  measured 
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FIG.  6.  Results  before  and  after  the  iterative  method  using  a  region  of  interest  map  which  consists  of  either  only  vessels  or  specific 
brain  structures  (in  this  case  the  basal  ganglia)  plus  vessels,  a:  The  initial  xi-o(r)  map  without  noise  added  in  the  original  simulated 
images,  b:  x/= 3(1)  map  of  (a)  using  threshold  value  a  =  0.1.  c:  Similar  to  (b),  a  =  0.2.  d:  The  initial  x/=oM  map  with  noise  added  in  original 
images,  resulting  a  standard  deviation  of  0.025  radian  in  phase  images,  e:  x;= 3W  map  of  (d)  using  a  =  0.1.  f:  The  associated  vessel 
map.  g:  The  x/=oM  map  in  the  coronal  plane  as  a  reference.  The  streaking  artifacts  are  clearly  shown  in  every  structure,  h:  The  \i= 3W 
maps  created  by  using  a  region  of  interest  map  which  consists  of  GR  SN,  RN,  PUT,  CN,  and  vessels,  i:  The  difference  image  of  (g)  and 
(h).  j:  The  initial  x/=o(r)  maP  in  the  transverse  plane  has  ‘false’  internal  capsule  (1C)  (pointed  by  an  arrow)  around  GP;  (k)  The  x/= 3(1)  map 
shows  no  “1C.”  This  matches  the  originally  simulated  model  (I).  No  noises  were  added  to  images  from  (g)  to  (I). 


from  a  uniform  region  inside  the  WM  decreased  from 
0.042  ppm  in  X;=oM  map  to  0.035  ppm  and  0.023  ppm 
in  the  Xr=3W  map  with  a  =  0.1  and  0.2,  respectively. 
The  baseline  susceptibilities  of  the  major  structures  are 
higher  with  SffARP  than  with  the  HP  filter.  The  iterative 
method  works  for  brain  structures  too  when  the  structure 
is  included  in  the  mask.  For  instance,  the  mean  suscepti¬ 
bility  values  of  the  GP  and  SN  have  been  increased  from 
0.155  ±  0.058  ppm  and  0.162  ±  0.067  ppm  in  the  Xr=oW 


map  to  0.163  ±  0.070  ppm  and  0.186  ±  0.083  ppm  in 
the  Xi— 3M  map,  from  the  dataset  processed  using 
SHARP.  The  result  after  HP  filtering  (Fig.  9e)  shows 
more  edge  artifacts  indicated  by  the  left  arrow  in  Fig.  9e. 
Much  of  this  error  was  reduced  by  the  iterative  method 
(Fig.  9f).  It  seems  that  the  iterative  method  compensated 
for  the  worse  first  guess  (Fig.  9e)  and  ended  up  with 
almost  the  same  result  (Fig.  9f,h)  as  having  started  with 
SHARP  (Fig.  9b, d)  from  the  image  perspective.  Since  a 
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FIG.  7.  The  plots  of  mean  susceptibility  values  inside  the  vessel 
(vein  of  Galen),  GP,  SN,  RN,  PUT,  and  CN  from  x/= 3W  maps.  The 
first  two  data  points  of  each  curve  is  the  value  inside  each  struc¬ 
ture  from  the  brain  model  and  the  X/=oM  map,  respectively,  a: 
X,= 3(r)  maps  generated  by  applying  a  region  of  interest  map  which 
consists  only  vessels  using  a  =  0.1,  0.15,  0.2,  0.25,  and  0.3, 
respectively,  b:  Xis(r)  maps  generated  by  applying  a  region  of 
interest  map  which  consists  of  the  GP,  SN,  RN,  PUT,  CN,  and 
vessels  using  a  =  0.1 . 


small  sized  HP  filter  cannot  remove  rapid  phase  wrap¬ 
ping  at  air-tissue  interfaces;  we  had  to  cut  out  the  region 
near  the  sinuses  in  the  phase  images. 


DISCUSSION 

In  this  article,  a  threshold-based  k-space/image  domain 
iterative  approach  has  been  presented.  Simulations  and 
in  vivo  results  show  that  the  ill-posed  problems  of 
streaking  artifacts  and  biases  in  the  estimates  of  suscepti¬ 
bilities  can  be  significantly  reduced.  The  replacement  of 
the  x(k)  values  near  the  singularities  by  XvmM-  which  is 
obtained  from  the  geometric  information  from  the  xM 
map  itself,  obviates  many  of  the  current  problems  seen 
in  the  TBSO  methods.  Since  XvmW  contains  little  streak¬ 
ing  artifacts  itself,  the  values  used  inside  the  thresholded 
regions  in  xffl  now  contain  no  artifact  either.  In  this 
sense,  we  obtain  an  almost  perfect  k-space  without  bad 
data  points  in  the  region  of  singularities.  This  explains 
why  this  method  converges  quickly  and  the  major 
improvement  is  in  the  first  iterative  step  (Fig.  5). 

The  proposed  iterative  approach  is  different  from  the 
other  threshold-based  methods  (11,15,18,19,24]  which 
fill  a  predefined  conical  region  using  a  constant,  zero  or 
1/thr  threshold  (11,15,24)  or  the  first-order  derivative  of 
g_1(k)  (19).  The  iterative  method  uses  full  geometry  in¬ 
formation  from  the  SM  (vessels  or  predefined  structures 
and  not  edge  information)  to  iteratively  change  k-space 
values  in  the  conical  region  using  the  forward  model. 
This  is  also  quite  different  than  other  currently  proposed 
solutions  (9,12,20,22).  Even  though  spatial  priors  such  as 
gradients  of  the  magnitude  are  used  (9,12,20,22),  in  those 
methods,  the  meaningful  values  of  the  singularity 
regions  in  k-space  are  obtained  through  solving  the  com¬ 
plex  cost  function  problem.  However,  the  iterative 
method  uses  priors  not  from  the  magnitude  image  but 
from  the  SM.  The  missing  data  in  the  singularity  regions 
are  obtained  through  iterating  back  and  forth  between 


0.02ppm 


-0.03ppm 
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FIG.  8.  Comparison  of  the  reconstructed  x/=3(r)  rnaps  using  (j)  accurate,  (k)  dilated,  and  (I)  eroded  vessel  maps.  Their  corresponding 
vessel  maps  and  the  enlarged  views  from  the  rectangular  regions  are  provided  in  (b)-(d)  and  (f)— (h).  (a)  and  (e)  The  initial  xi-o(r)  maps 
and  (i)  the  original  brain  model  as  references.  The  circle  in  the  midbrain  in  the  x(r)  maps  represents  the  RN  and  is  indicated  by  a  black 
arrow  in  (i).  Other  hyper-intense  regions  in  SMs  are  vessels. 
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FIG.  9.  Comparisons  of  SMs  using  SHARP  or  a  HP  filter  (FWHM  =  16  pixels)  to  remove  the  background  field.  The  iterative  method 
with  a  =  0.1  and  0.2  is  applied  after  the  background  is  removed,  (a)-(d)  and  (e)-(h)  are  results  after  the  application  of  SHARP  and  the 
HP  filter,  respectively,  (a)  and  (e)  the  initial  x/=oM  rnaps.  (b)  and  (f)  the  xi= 3W  maps  generated  from  the  iterative  method  with  a  =  0.1.  (c) 
and  (g)  the  differences  of  images  between  (a)  and  (b),  and  between  (e)  and  (f),  respectively.  These  two  images  show  the  successful 
reduction  of  the  streaking  artifacts,  (d)  and  (h)  the  x/= 3(1)  maps  generated  from  the  iterative  method  with  a  =  0.2.  The  range  of  the  gray¬ 
scale  bars  is  chosen  to  highlight  the  artifacts  in  the  images.  It  does  not  reflect  the  quantified  susceptibility  values  inside  veins. 


the  SMs  and  their  k-space.  The  advantage  of  cost  func¬ 
tion  approaches  is  that  they  do  not  need  to  predefine  the 
singularity  region  in  k-space  which  is  solved  by  the  opti¬ 
mization  process  automatically  (although  the  optimiza¬ 
tion  process  itself  is  usually  quite  time-consuming].  On 
the  other  hand,  the  iterative  method  is  the  most  time-ef¬ 
ficient.  It  is  fast  enough  to  reconstruct  SMs  for  a  512  x 
512  x  256  data  set  using  an  Intel  Core  i7  CPU  3.4  GHz 
processor  in  less  than  30  s,  since  in  practice  usually 
three  iterations  are  good  enough  to  generate  decent 
results. 

The  threshold  value  also  plays  a  key  role.  A  threshold 
value  of  0.1  is  a  reasonable  choice  since  a  lower  thresh¬ 
old  value  leads  to  an  increase  in  noise  and  a  higher 


threshold  value  leads  to  a  blurring  of  the  object  (Figs. 
6c,9d,h). 

It  is  known  that  the  ill-posedness  of  the  inverse  filter 
will  increase  the  noise  level  from  the  phase  to  the  SM. 
Based  on  both  simulations  and  real  data,  we  find  that 
there  is  a  factor  of  4  increase  in  noise  in  the  SM  relative 
to  the  original  phase  data.  This  result  and  the  fact  at  B0 
=  3  T,  Te  =  5  ms,  and  aXJ=0(r)  =  0.025  ppm  make  it  pos¬ 
sible  to  write  the  total  noise  in  the  background  region  in 
Xi=oM  as  0.025-4  (3/Bo)  (5/TE)  in  ppm.  The  noise  in 
Xj=3M  wiH  be  less  than  this  value  since  the  iterative 
method  will  reduce  streaking  artifact  in  SM. 

The  iterative  method  can  be  used  to  remove  streaking 
artifacts  associated  with  not  only  vessels  but  also  other 
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Table  1 

Ax  Measured  In  Vivo  in  Two  Veins  in  xi=o  if)  Maps  and  x/= 3  M 
Maps  With  Different  Threshold  Values 


50= o  ( r ) 
map 

50=3  ( r ) 

map /a  =  0.1 

50=3  (r) 
map /a  =  0.2 

VI  (SHARP) 

0.32  ±  0.07 

0.37  ±  0.08 

0.38  ±  0.09 

VI  (HP) 

0.24  ±  0.05 

0.28  ±  0.06 

0.28  ±  0.06 

V2  (SHARP) 

0.35  ±  0.04 

0.40  ±  0.05 

0.41  ±  0.05 

V2  (HP) 

0.25  ±  0.05 

0.31  ±  0.06 

0.30  ±  0.06 

Mean  and  standard  deviation  for  the  susceptibility  values  (in  ppm) 
of  two  veins  processed  using  SHARP  and  a  HP  filter  (FWHM  =  16 
pixels),  respectively,  were  chosen  from  the  0.5  mm  isotropic  reso¬ 
lution  data.  VI  and  V2  are  shown  in  Fig.  9.  The  susceptibility  val¬ 
ues  of  these  two  veins  have  been  increased  by  the  iterative 
method.  There  is  not  much  variation  of  the  susceptibility  value 
with  different  threshold  values. 


brain  structures  as  well.  Figure  6h  shows  a  reduction  in 
artifacts  associated  specifically  with  iron-rich  regions 
such  as  the  GP  and  CN. 

Accurately  extracting  vessels  from  the  Xr=o(r)  map  is 
critical  for  the  iterative  method  (Fig.  8).  In  this  study, 
vessels  were  segmented  directly  from  the  SM  (Fig.  2).  It 
may  also  possible  to  segment  veins  from  original  magni¬ 
tude  images  (9,12,20,22),  phase  images,  and/or  SWI 
images.  Extraction  of  accurate  anatomic  information  from 
phase  data  sometimes  is  difficult  since  phase  is  orienta¬ 
tion  dependent  and  phase  changes  are  generally  nonlo¬ 
cal.  SWI  images  work  better  for  an  anisotropic  dataset 
rather  than  an  isotropic  dataset  since  phase  cancellation 
is  needed  to  highlight  vessel  information.  Therefore,  we 
may  consider  combining  SMs  with  magnitude  images, 
phase  images,  and/or  SWI  images  together  to  segment 
the  veins,  since  different  types  of  images  can  compensate 
for  missing  information. 

The  iterative  method  appears  to  help  even  in  the  pres¬ 
ence  of  non-isotropic  resolution  with  partial  volume 
effects  and  to  a  minor  degree  when  an  HP  filter  is 
applied.  A  smaller  sized  HP  filter  would  be  better,  since 
a  larger  HP  filter  will  significantly  underestimate  the 
susceptibility  value  (Table  1).  SHARP  gave  us  better 
results  compared  with  the  HP  filter  (FWHM  =  16  pixels) 
(Fig.  9),  but  SHARP  requires  phase  unwrapping  which 
can  be  time  consuming  and  is  noise  dependent  (19). 
From  this  perspective,  an  HP  filter  has  the  advantage 
since  it  does  not  need  unwrapped  phase.  If  the  forward 
modeling  approach  of  Neelavalli  et  al.  (38)  can  be  used 
to  reduce  air/tissue  interface  fields,  then  it  may  be  possi¬ 
ble  to  use  a  small  size  HP  filter  (FWHM  =  8  pixels) 
which  may  provide  similar  results  to  SHARP. 

Severe  streaking  artifacts  associated  with  structures 
having  high  susceptibility  values  such  as  veins  can  lead 
to  major  changes  in  the  appearance  of  the  brain  struc¬ 
tures  with  low  susceptibility.  Practically,  the  susceptibil¬ 
ity  of  the  veins  is  a  factor  of  2.5-20  times  higher  than 
other  structures  in  the  brain.  Therefore,  even  a  10% 
streaking  artifact  can  overwhelm  the  information  in  the 
rest  of  the  brain  and  create  false  appearing  structures  as 
in  (Fig.  6j)  and  in  (Fig.  8e).  The  reduction  of  these  arti¬ 
facts  makes  a  dramatic  difference  in  the  ability  to  prop¬ 
erly  extract  the  susceptibility  of  other  tissues. 


In  conclusion,  both  simulations  and  human  studies 
have  demonstrated  that  the  proposed  iterative  approach 
can  dramatically  reduce  streaking  artifacts  and  improve 
the  accuracy  of  susceptibility  quantification  inside  the 
structures  of  interest  such  as  veins  or  other  brain  tissues. 
Given  its  relatively  fast  processing  time,  it  should  be 
possible  to  expand  its  use  into  more  daily  clinical  prac¬ 
tice.  With  the  improved  accuracy  of  the  susceptibility 
values  inside  veins,  this  method  could  be  used  poten¬ 
tially  to  improve  quantification  of  venous  oxygen  satura¬ 
tion  (18). 
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The  development  of  functional  brain  magnetic  resonance  imaging  (fMRI)  has  been  a  boon  for  neuroscientists 
and  radiologists  alike.  It  provides  for  fundamental  information  on  brain  function  and  better  diagnostic  tools 
to  study  disease.  In  this  paper,  we  will  review  some  of  the  early  concepts  in  high  resolution  gradient  echo  im¬ 
aging  with  a  particular  emphasis  on  susceptibility  weighted  imaging  (SWI)  and  MR  angiography  (MRA).  We 
begin  with  the  history  of  our  own  experience  in  this  area,  followed  by  a  discussion  of  the  role  of  high  resolu¬ 
tion  in  studying  the  vasculature  of  the  brain  and  how  this  relates  to  the  BOLD  (blood  oxygenation  level  de¬ 
pendent)  signal.  We  introduce  the  role  of  SWI  and  susceptibility  mapping  (SWIM)  in  fMRI  and  close  with 
recommendations  for  future  high  resolution  experiments. 

©  2012  Elsevier  Inc.  All  rights  reserved. 
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Introduction 

In  the  early  1990s,  our  research  interests  were  in  developing  and 
understanding  3D  gradient  echo  imaging  methods  to  study  neurolog¬ 
ical  diseases.  It  would  be  no  surprise  to  the  reader  that  we  were  de¬ 
lighted  to  see  the  early  work  from  Seiji  Ogawa  at  Bell  Labs  (Ogawa 
et  al„  1993)  and  Jack  Belliveau  and  colleagues  at  MGH  (Belliveau  et 
al„  1991)  regarding  the  use  of  magnetic  resonance  imaging  (MRI) 
for  functional  MRI  (fMRI)  of  the  brain.  The  news  came  to  our  atten¬ 
tion  quickly,  as  we  were  then  focusing  on  the  effects  of  changes  in 
local  fields  in  gradient  echo  imaging  at  Case  Western  Reserve  Univer¬ 
sity  (CWRU)  more  as  artifacts  than  as  a  useful  tool.  At  that  time,  we 
had  interests  in  017  imaging  (Hopkins  et  al„  1988;  Kwong  et  al„  1991) 
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using  gradient  echo  imaging  for  brain  function  and  Amos  Hopkins,  the 
major  instigator  in  this  direction,  became  interested  immediately  in  col¬ 
laborating  with  people  at  MGH  on  the  concept  of  fMRI  with  O17  imag¬ 
ing.  Some  of  the  experiments  done  in  this  area  showed  that  we  could 
see  major  changes  in  the  pial  veins  between  activation  and  resting 
state  (Frahm  et  al„  1993;  Haacke  et  al„  1994;  Lai  et  al.,  1993),  and  that 
the  early  explanation  of  diffusion  as  the  source  of  signal  change  in 
fMRI  was  likely  not  correct.  This  began  what  was  then  known  as  the 
brain/vein  debate. 

It  took  some  years  later  and  a  number  of  key  papers  to  show  that 
the  intravascular  and  extravascular  effects,  at  least  at  1 .5  T,  were  the 
major  sources  of  signal  change  in  fMRI.  To  demonstrate  this,  re¬ 
searchers  applied  strong  dephasing  gradients  and/or  saturation 
pulses  to  basically  remove  the  arterial  signal  from  the  images  and 
then  evaluated  the  fMRI  response  (Duong  et  al.,  2003).  After  suppres¬ 
sing  blood  flow,  there  was  little  fMRI  response  remaining,  implicating 
a  major  intravascular  source.  An  example  of  the  effect  of  saturation 
pulses  applied  outside  the  slices  of  interest  for  a  3  T  fMRI  experiment 
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Fig.  1.  GE-  (a,  b)  and  SE-EP1  (c,  d)  t-maps  obtained  without  (a,  c)  and  with  (b,  d)  blood 
saturation,  showing  a  reduced  BOLD  response  when  blood  signal  is  saturated.  TE  of  GE- 
and  SE-EPI  was  24  ms  and  96  ms  respectively.  The  blood  saturation  was  done  by  plac¬ 
ing  two  saturation  bands  above  and  below  the  imaging  slab,  and  the  visual  stimulation 
was  8  Hz  flashing  checkerboard  with  30  s/30  s  on/off  block  design. 


in  the  visual  cortex  is  shown  in  Fig.  1 .  Clearly,  there  is  a  decrease  in 
the  signal  even  for  the  gradient  echo  EPI  sequence.  It  should  be 
noted  that  for  the  long  TR  experiments  used  for  whole  brain  multi¬ 
slice  EPI,  there  is  a  strong  inflow  of  fresh  blood,  which  is  referred  to 


as  “inflow  BOLD  enhancement”  (Duong  et  al„  2003;  Duyn  et  al., 
1994).  By  suppressing  the  signal  from  blood,  the  inflow  BOLD  en¬ 
hancement  effects  are  also  suppressed,  leaving  only  the  extravascular 
effects.  During  the  last  20  years,  many  papers  have  been  written 
studying  the  role  of  intravascular  and  extravascular  BOLD  effects, 
even  up  to  today,  especially  since  higher  and  higher  field  strengths 
have  become  available  and  more  and  more  signal  is  generated  from 
extravascular  effects  (Baudendistel  et  al„  1998;  Boxerman  et  al., 
1995;  Casciaro  et  al.,  2008;  Song  et  al.,  2007;  Ugurbil  et  al.,  2000). 

The  role  of  high  resolution 

One  of  our  goals  has  been  to  look  for  methods  to  increase  resolu¬ 
tion  and  yet  keep  imaging  time  reasonable.  One  early  example  of  this 
was  partial  Lourier  imaging  (Haacke  et  al.,  1999;  Xu  and  Haacke, 
2001).  This  approach  allowed  us  (and  still  does)  to  take  an  asymmet¬ 
ric  echo,  keep  flow  compensation  for  short  echoes,  and  collect  many 
points  after  the  echo  to  maintain  high  resolution  in  the  read  direction. 
We  also  focused  on  segmented  EPI  to  improve  speed  for  high  resolu¬ 
tion  imaging  (Haacke  et  al.,  1986;  Xu  and  Haacke,  2008).  We  believed 
that  new  information  awaited  us  structurally  and  functionally  if  we 
could  push  resolution  down  to  100  pm.  Today,  high-field  and  parallel 
imaging  are  re-opening  the  door  to  that  goal.  The  group  in  Minnesota 
pioneered  this  direction  in  the  visual  cortex  by  studying  ocular  dom¬ 
inance  columns  with  high  resolution  fMRI  (Menon  et  al.,  1997; 
Yacoub  et  al.,  2008).  Of  course,  the  most  interesting  approach  would 
be  to  create  high  resolution  functional  imaging  with  anatomical  qual¬ 
ities,  whether  it  be  vascular,  metabolic,  or  brain  function.  No  doubt 
researchers  in  this  field  will  continue  to  fine  tune  the  MR  microscope 
to  push  these  limits  further  and  further  in  the  future  thanks  to  the 
continued  technical  developments  that  lead  to  higher  signal-to- 
noise,  faster  imaging  and  better  resolution. 

Why  high  resolution?  Evidently  the  workhorse  for  fMRI  today  re¬ 
mains  a  low  resolution  scan  with  roughly  3  mm  resolution  in  each  di¬ 
rection.  The  argument  for  this  is  that  the  signal-to-noise  (SNR)  has  to 
be  high  enough  that  one  can  detect  changes  in  signal  as  small  as  1%. 
However,  it  depends  on  whether  or  not  there  is  a  scale  invariant 


Fig.  2.  Activations  in  the  primary  motor  cortex  using  a  block  design  finger  tapping  paradigm,  (a)  Filtered  3D  GRE  images  of  single  subtraction  between  activation  and  resting  state, 
with  voxel  size  of  1.75  x  1.75  x  2  mm3;  (b)  phase  image  from  the  same  slice  showing  the  motor  cortex;  (c)  t-value  map  calculated  from  the  GE-EPI  data  acquired  using  3.5  mm  iso¬ 
tropic  voxels  and  a  TE  of  30  ms;  (d)  same  as  in  (c)  but  with  modified  window  level  settings;  (e)  sagittal  reformat  of  the  3D  GE  subtraction;  and  (f)  sagittal  MP-RAGE  image.  Imaging 
parameters  for  the  3D  GRE  sequence  were:  TR/TE  =  30/20  ms,  flip  angle  =  1 5°,  BW  =  80  Hz/pixel,  2  x  GRAPPA,  resolution  =  0.875  x  0.875  x  2  mm3,  Nz  =  64,  acquisition  time  for  each 
block  was  96  s.  The  average  signal  change  in  GE-EPI  activation  (d)  is  about  0.8%  with  maximum  of  3.5%,  and  the  signal  change  in  the  3D  GRE  data  (e)  is  about  14%. 
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Fig.  3.  (a)  MaxlP  images  of  a  normal  TOF;  (b)  TOF  with  a  fixed  arterial  saturation  pulse;  and  (c)  their  subtraction  image.  The  imaging  parameters  were:  0.5  x  0.5  x  0.5  mm3  isotropic 
voxels,  TR/TE  =  20/1 1  ms,  flip  angle  =  15°,  and  fully  flow  compensated.  All  MaxlP  was  done  over  1 6  slices.  Although  most  major  vessels  are  suppressed,  one  can  see  some  remnant 
arterial/venous  signal  even  with  saturation.  The  reduced  signal  from  the  veins  in  this  case  actually  leads  to  a  poorer  signal  cancellation  in  the  usual  SW1  data  (not  shown  here) 
which  relies  on  having  a  strong  intravascular  component. 


response.  For  capillary  structures,  resolution  is  probably  not  impor¬ 
tant.  But  what  are  we  really  seeing  in  current  low  resolution  EP1 
fMRl  results?  In  1999,  Hoogenraad  et  al.  (1999)  showed  that  a  very 
high  resolution  gradient  echo  scan  could  be  used  to  collect  fMRI 
data  (albeit  much  more  slowly  than  EP1)  and  then  collapsed  back  to 
the  same  low  resolution  as  EPI  and  that  the  results  were  basically 
the  same.  This  suggests  that  the  signal  from  the  major  veins  still 
plays  a  key  role  in  the  fMRI  response  even  for  low  resolution  imaging 
(see  Fig.  2). 

As  mentioned  earlier,  the  only  way  to  really  remove  the  intravas¬ 
cular  effects  is  by  saturation  or  dephasing  or  both,  and  in  this  way  one 
can  be  sure  that  the  response  being  imaged  is  predominantly  extra- 
vascular.  But  saturation  alone  may  not  explain  the  fact  that  there  re¬ 
mains  some  signal  from  what  appears  to  be  the  same  location  as 
before  in  Fig.  1  but  just  reduced  in  amplitude.  Consider  the  high  res¬ 
olution  angiogram  in  Fig.  3.  The  role  of  saturation  can  be  appreciated 
in  this  figure  and  its  effect  is  demonstrated  by  subtraction  of  the  sat¬ 
urated  data  from  the  unsaturated  data  (Fig.  3c).  The  second  image 
(Fig.  3b)  shows  that  the  signal  from  the  arteries  and  veins  is  almost 
isointense  to  the  surrounding  tissue  but  certainly  not  completely 
nulled  either  due  to  Ti  recovery.  To  remove  any  remnant  effects 
from  the  venous  blood,  one  could  then  apply  dephasing  gradients. 
This  should  null  any  remaining  intravascular  contribution  to  the  sig¬ 
nal  (see  also  (Song  et  al.,  1996)  and  Fig.  4).  This  was  indeed  the  argu¬ 
ment  in  the  early  brain/vein  debate. 


The  results  in  Fig.  2  however  are  very  suggestive.  They  are  identi¬ 
cal  in  nature  to  the  earlier  results  from  Frahm  et  al.  (1993),  Haacke  et 
al.  (1994),  Lai  et  al.  (1993),  only  acquired  faster  at  3  T.  The  unfiltered 
EPI  acquired  with  lower  resolution  show  pretty  much  the  same  areas 
of  activation.  The  filtered  EPI  activation  t-maps  have  lost  the  ability 
for  good  spatial  localization  and  show  activation  badly  blurred  over 
that  portion  of  the  brain.  Perhaps  one  of  the  best  clinical  applications 
for  these  high  resolution  fMRI  methods  would  be  in  the  study  of  elo¬ 
quent  cortex  when  the  neurosurgeon  must  operate  on  a  tumor 
(Baudelet  et  al.,  2006;  Hou  et  al.,  2006;  Sherman  et  al.,  2011).  Given 
the  results  seen  in  the  3D  data  with  complete  visualization  of  the 
brain,  one  wonders  why  this  method  has  not  been  used  more  in  the 
last  20  years.  Can  we  use  high  resolution  gradient  echo  methods  prac¬ 
tically?  To  understand  when  this  is  possible,  we  need  to  first  address 
the  role  of  vessel  size. 

The  work  by  Cheng  and  Haacke  (2001)  tried  to  address  this  ques¬ 
tion  of  scale  variance  more  carefully,  showing  that  the  expected  fMRI 
response  is  scale  invariant  when  the  voxels  are  larger  than  the  largest 
veins,  hence  making  EPI  the  method  of  choice  after  that  point 
(2x2x2  mm3  is  probably  large  enough  to  meet  that  limit  and 
would  provide  far  better  spatial  delineation  than  what  people  are 
using  today  which  is  often  just  a  broadly  spread  out  inflow  enhanced 
vascular  response).  However,  Haacke  et  al.  (1994)  and  Cheng  and 
Haacke  (2001)  also  showed  that  when  the  scale  becomes  on  the 
order  of  the  vessel  size  and  smaller,  there  will  be  a  dramatic  increase 


Fig.  4.  (a)  MRA  MaxlP  images  acquired  without  dephasing  gradients,  and  (b)  MinlP  images  acquired  with  dephasing  gradients,  with  a  VENC  value  of  0.94  cm/s.  Both  arteries  and 
even  medullary  veins  are  well  suppressed,  (c)  is  the  result  of  subtracting  (b)  from  (a),  removing  most  of  the  background  tissues  while  accentuating  the  vessels.  All  MaxlP  was  done 
over  16  slices  of  0.5  mm  thickness.  In  this  subtraction  both  arteries  and  veins  are  evident  although  the  venous  structures  are  not  as  bright  as  the  arteries  structures  in  part  because 
the  signal  from  the  veins  is  already  suppressed  because  of  their  low  T2*  value. 
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in  signal  response.  In  fact,  in  Lai  et  al.  (1993),  the  authors  showed  an 
image  from  Duvernoy  et  al.  (1981 )  showing  the  pial  vein,  venules  and 
capillaries,  and  commented  that  as  resolution  improves  so  will  the 
ability  to  localize  the  fMRI  signal  once  the  resolution  is  on  the  order 
of  the  size  of  the  vessel  (vein)  of  interest.  The  general  argument 
against  this  has  been  that  there  is  not  enough  SNR  or  temporal  reso¬ 
lution  to  collect  the  data  at  such  high  resolution  (and  Cheng  and 
Haacke's  paper  even  suggests  this  to  a  point).  However,  this  has  all 
changed  with  the  advent  of  high  fields  such  as  3  T  and  even  more 
so  with  7  T,  where  very  high  SNR  and  very  high  resolution  fMRI  has 
become  viable.  It  is  a  testament  to  the  progress  made  in  MRI  that 
the  picture  from  Duvernoy  showing  the  pial  veins,  venules  and  capil¬ 
laries  has  basically  been  replicated  now  with  SWI  filtered  phase 
(Haacke  et  al.,  1995;  Haacke  and  Reichenbach,  2011)  and  gradient 
echo  images  at  7  T  (Ge  et  al.,  2008;  Lee  et  al.,  2010)  and  is  beginning 
to  be  seen  at  3  T  (Haacke  and  Reichenbach,  2011 ). 

In  the  early  days,  much  work  went  into  suppressing  the  signal 
from  the  veins  (Duong  et  al.,  2003;  Glover  et  al.,  1996).  But  now  it 
is  possible  to  do  what  was  alluded  to  in  the  work  of  Haacke  et  al. 
(1994),  image  the  small  draining  veins  adjacent  to  the  activated  tis¬ 
sue,  especially  if  the  goal  is  to  do  the  neuroscience  to  understand 
the  sources  of  the  local  neural  activity  involved  in  activation.  With 
these  high  SNR,  high  resolution,  rapid  scanning  methods,  it  should 
be  possible  to  begin  to  see  local  changes  in  venous  oxygen  saturation 
(which  is  the  underlying  mechanism  of  BOLD  effect)  at  the  level  of 
medullary  veins  (several  hundred  microns)  to  perhaps  even  local 
changes  in  venules  (50  pm).  It  turns  out  that  the  optimal  resolution 
according  to  Cheng  and  Haacke  is  about  four  times  larger  than  the 
vessel  volume.  So  a  vein  that  is  250  pm  could  be  seen  and  imaged 
with  a  resolution  of  500  pm  (as  it  is  now  with  SWI  for  example 
(Belliveau  et  al.,  1992;  Haacke  et  al.,  1997,  2004,  2009;  Mittal  et  al., 
2009)  or  a  vessel  that  is  50  pm  could  be  seen  with  a  resolution  of 
100  pm.  Practically,  the  results  are  better  than  Cheng  and  Haacke  had 
anticipated.  The  7  T  data  of  today  can  already  reveal  venules  with  a  res¬ 
olution  of  0.2  x  0.2  x  1  mm3  (Geetal.,2008;  Leeetal.,  2010).  We  can  ex¬ 
pect  to  see  many  more  publications  in  the  future  on  high  resolution 
functional  MRI  using  3D  gradient  echo  imaging  techniques. 

There  are  some  practical  issues  that  need  to  be  dealt  with  to  make  im¬ 
aging  with  high  resolution  more  attractive.  One  of  these  is  scanning  time 
and  the  other  is  SNR.  To  cover  the  entire  brain  at  3  T  with  64  2  mm  thick 
slices  and  an  in-plane  resolution  of  0.5  mm  x  0.5  mm  will  take  about 

5  min.  Keeping  a  single  activation  task  running  for  5  min  is  not  easy  and 
then  of  course  another  5  min  to  obtain  the  resting  state  images.  However, 
with  the  marriage  of  parallel  imaging  and  segmented  echo-planar  imag¬ 
ing,  we  anticipate  this  time  being  reduced  to  one  and  a  half  minutes.  Re¬ 
ducing  the  resolution  slightly  to  1  mm  x  1  mm  x  2  mm  would  further 
reduce  the  time  by  about  30  s.  Now  this  allows  for  a  more  conventional 
block  design  of  30  s  activation  and  30  s  resting  for  five  cycles,  allowing 
the  entire  3D  scan  to  take  place  in  just  5  min.  This  is  similar  to  current  im¬ 
aging  times  in  most  fMRI  experiments. 

High  resolution  MR  angiography:  the  role  of  saturation  and 
dephasing 

Before  we  address  the  practical  issue  of  high  resolution  fMRI  with 
3D  gradient  echo  imaging  and  SWI,  it  is  important  to  understand 
what  is  possible  on  the  MR  angiography  side  from  the  perspective 
of  collecting  high  resolution  data,  enhancing  contrast,  and  nulling  sig¬ 
nal  from  blood.  In  fMRI  with  EPI,  one  gets  excellent  intravascular  sig¬ 
nal  contributions  thanks  to  the  inflow  enhancement  from  the  arterial 
side  and  hence  venous  signal.  But  with  gradient  echo  imaging,  blood 
signal  will  be  suppressed  with  respect  to  gray  matter  due  to  blood’s 
longer  T,.  One  can  do  MRA  both  with  and  without  contrast  agent  to 
best  enhance  the  vascular  signal  and  more  importantly,  the  contrast 
between  vessels  and  surrounding  tissue.  Recall  that  if  the  vascular 
signal  is  too  low,  the  enhancement  from  the  intravascular  BOLD  effect 


will  be  minimized.  Fig.  3  and  particularly  Fig.  4  show  the  capability  of 
high  resolution  imaging  in  mapping  out  not  just  the  major  arteries 
but  also,  most  importantly,  smaller  branches  that  are  the  basic  source 
of  blood  supply  to  the  tissue. 

These  images  can  be  further  improved  by  using  a  magnetization 
transfer  contrast  (MTC)  pulse,  using  a  saturation  band,  dephasing 
the  signal  from  the  moving  spins,  or  other  means  to  suppress  the 
background.  Consider  the  use  of  MTC  as  a  first  approach.  The  problem 
with  MTC  is  that  it  lengthens  the  scan  time  and  increases  the  specific 
absorption  rate  (SAR)  making  it  less  desirable  to  use  at  high  magnetic 
fields.  A  second  approach  is  to  perform  two  scans:  one  with  high  ves¬ 
sel  signal  and  the  other  with  low  signal  and  subtract  them.  One  can 
use  saturation  bands  to  suppress  the  arterial  signal  or  venous  signal 
by  saturating  the  inflowing  blood,  or  use  a  strong  set  of  bipolar  gradi¬ 
ents  to  dephase  moving  spins  in  both  arteries  and  veins.  The  advan¬ 
tage  of  the  bipolar  gradients  over  saturation  bands  is  that  it 
suppresses  fast  flowing  blood,  requires  no  SAR  increases  and  can  be 
acquired  in  an  interleaved  TR  manner  (i.e.,  one  TR  with  dephasing 
and  the  next  TR  without)  without  interfering  with  equilibrium  of 
the  magnetization  and  without  suffering  from  misregistration  arti¬ 
facts.  The  VENC  concept  in  flow  quantification  methods  (Haacke  et 
al.,  1999)  can  be  used  to  describe  the  dephasing  effects,  that  flowing 
blood  with  a  velocity  higher  than  the  VENC  value  will  be  dephased. 
Setting  the  VENC  values  to  about  1  cm/s  or  less  should  do  a  reason¬ 
ably  good  job  dephasing  laminar  flow  even  for  small  vessels  (250  to 
500  pm  sized  arteries  and  veins).  The  disadvantage  is  that  it  still  re¬ 
quires  twice  the  time  (except  in  some  special  double  or  triple  echo  se¬ 
quence  designs  (Kimura  et  al.,  2009)),  and  both  arteries  and  veins 
will  be  shown  in  the  subtracted  images.  Fig.  4a  shows  an  example 
of  the  usual  arterial  inflow  effects,  Fig.  4b  a  MinIP  image  showing 
that  most  vessels  are  dephased,  and  Fig.  4c  the  subtraction  of  Fig.  4b 
from  4a.  As  a  side  note,  those  in  the  fMRI  field  often  quote  a  b-value 
rather  than  the  VENC  value.  However,  the  bipolar  gradients  are  not 
being  used  as  diffusion  gradients,  nor  is  the  diffusion  effect  the  under¬ 
lying  mechanism  for  signal  suppression  here.  Therefore,  we  have 
used  the  more  conventional  nomenclature  from  MR  angiography  of 
VENC  value  above  which  aliasing  occurs.  With  a  pair  of  bipolar  rect¬ 
angular  gradients,  for  simplicity,  one  can  relate  the  VENC  value  to 
the  b-value  via  VENC  =  n/(6bT)1/2  (since  VENC  =  tt/(27Gt2)  and 
b  =  2-y2G2T3/3),  where  t  is  the  duration  of  one  gradient  lobe.  So,  for 
example,  if  the  applied  gradient  is  25  mT/m  with  5  ms  duration, 
then  b  =  3.73  s/mm2  and  VENC  =  0.94  cm/s.  Since  the  peak  flows  in 
pial  veins  is  on  the  order  of  1-4  cm/s,  one  can  expect  their  signal  to 
be  fully  suppressed  with  a  VENC  of  0.94  cm/s. 

As  a  fourth  approach,  a  contrast  agent  can  be  used  to  enhance  ves¬ 
sel  signal  throughout  the  brain  without  being  concerned  about  satu¬ 
ration  effects  even  for  slower  flowing  blood  in  small  vessels.  The 
challenge  in  imaging  the  vessels  in  vivo  is  to  strive  to  mimic  the  ex¬ 
ample  of  the  arteries  from  Salamon  (1971)  as  shown  in  Fig.  5a.  In 
this  sagittal  view,  one  can  see  the  radial  branching  of  the  arteries 
into  the  brain.  Many  of  the  major  arteries  that  are  a  few  millimeters 
in  size  bifurcate  into  arteries  that  are  perhaps  only  100  to  200  pm  in 
size.  These  in  turn  branch  into  vessels  only  a  few  dozen  microns  in 
size  and  these  finally  to  the  capillaries.  Using  the  high  resolution 
MRA  approach  with  contrast  agent  injection  yields  the  comparative 
results  shown  in  Fig.  5b.  Although  no  dephasing  example  has  been 
used  to  subtract  away  the  remaining  background  tissue,  there  is  still 
clear  evidence  of  small  arteries  in  these  images  that  are  approaching 
100  to  200  pm  in  size.  A  key  component  in  understanding  the  role  of 
the  vasculature  in  BOLD  imaging  is  to  be  able  to  map  out  the  arteries 
and  veins  separately.  Further,  since  a  major  component  of  the  BOLD 
signal  comes  from  intravascular  effects,  the  use  of  a  Ti  reducing  con¬ 
trast  agent  to  increase  the  venous  signal  should  in  principle  lead  to  an 
improvement  in  the  BOLD  response. 

The  thrust  for  high  resolution  MR  angiographic  approaches  is  not 
just  to  better  understand  the  vascular  structure  and  source  of  vascular 
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Fig.  5.  (a)  High  resolution  radiographic  image  of  the  arteries  obtained  by  injecting 
radio-opaque  dye  in  a  cadaver  brain  within  a  few  hours  after  death  (Salamon,  1971). 
(b)  A  representative  sagittal  MaxIP  image  from  a  high  resolution  MRA  dataset  acquired 
in  a  transverse  plane  post-contrast  with  a  voxel  size  of  0.25  x  0.5  x  0.5  mm3  and  then 
reformatted  to  match  this  saggital  view.  Note  the  radial  pattern  of  the  arteries  into 
the  thalamus  and  surrounding  regions.  Many  very  small  vessels  can  be  seen  on  the 
order  of  only  a  few  hundred  microns.  No  subtraction  was  performed  here. 

signal,  but  also  to  consider  the  novel  approach  of  using  the  arterial 
signal  as  a  means  to  study  inflow  changes  induced  by  the  activation 
state.  This  has  been  considered  in  the  past  for  arterial  spin  labeling 
applications  in  fMRI  (Detre  and  Wang,  2002;  Obata  et  al.,  2004)  and 
more  recently  for  task  induced  arterial  inflow  changes  (Kang  et  al., 
2010).  Being  able  to  analyze  the  inflow  as  well  as  cerebral  blood  re¬ 
serve  (such  as  CBV  and  CBF)  may  sente  as  a  bridge  to  connect  arterial 
input  and  venous  BOLD  effects,  as  well  as  to  understand  the  underly¬ 
ing  vasculature  basis  for  the  significant  inter-individual  variation 
commonly  seen  in  BOLD  fMRI. 

Susceptibility  weighted  imaging  (SWI) 

Throughout  this  early  period  our  interest  was  to  understand  BOLD 
imaging  by  studying  the  veins.  This  focus  eventually  led  to  pushing 
the  resolution  limits  until,  in  1997  (Reichenbach  et  al.,  1997),  we 
came  up  with  the  concept  of  susceptibility  weighted  imaging  (SWI). 
One  of  the  key  participants  in  this  work  was  Jurgen  Reichenbach,  then 
a  post-doctoral  fellow  in  St.  Louis,  Missouri  at  the  Mallinckrodt  Institute 
of  Radiology  (he  is  now  a  full  Professor  running  his  own  imaging  labo¬ 
ratory  at  the  Friedrich-Schiller  Universitat  in  Jena,  Germany).  Other  re¬ 
searchers  who  were  involved  at  this  time  included  Song  Lai  (Lai  et  al., 
1993)  and  Frank  Hoogenraad  (Hoogenraad  et  al.,  1999).  Although  the 
precursor  to  this  had  been  a  paper  in  1995  studying  the  susceptibility 
effects  that  could  be  visualized  with  phase  (Haacke  et  al„  1995),  the  po¬ 
tential  clinical  applications  started  first  with  SWI.  During  the  next  ten 
years,  we  saw  SWI  show  its  value  in  a  wide-ranging  set  of  conditions 
such  as  aging,  multiple  sclerosis,  stroke,  trauma  and  tumors  (Mittal  et 
al.,  2009).  In  1997,  we  used  phase  for  both  flow  and  susceptibility  to 
study  changes  in  oxygen  saturation  in  fMRI  experiments  using  a 
method  we  referred  to  as  SAIF  for  susceptibility  and  interleaved 
flow  imaging  (Haacke  et  al.,  1997). 

In  the  late  1990s,  we  proposed  that  multi-echo  fMRI  (Barth  et  al., 
1999)  would  be  of  interest  and  would  serve  as  a  new  tool  to  investi¬ 
gate  the  echo  time  dependence  of  the  BOLD  response  and  how  it  led 
to  different  spatial  response  patterns.  In  that  work,  we  also  evaluated 
phase  data  showing  that  phase  could  serve  as  a  good  marker  for  sus¬ 
ceptibility  (oxygen  saturation)  changes  in  fMRI.  In  1995  (Haacke  et 
al„  1995),  1997  (Haacke  et  al.,  1997;  Reichenbach  et  al.,  1997),  and 
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in  2004  (Haacke  et  al.,  2004),  we  published  again  with  an  effort  to 
push  the  SWI  filtered  phase  as  a  useful  tool  in  itself.  We  have  also 
used  the  phase  as  a  means  to  map  the  primary  motor  cortex  as  the 
iron  content  in  the  gray  matter  is  quite  high  in  this  region  and  the 
phase  images  give  excellent  gray  matter/white  matter  contrast 
there  (Ogg  et  al.,  1999).  More  recently,  there  has  been  increased  in¬ 
terest  and  attention  paid  to  phase  in  fMRI  (Rowe  and  Haacke,  2009; 
Rowe  and  Logan,  2004).  The  final  thrust  in  this  direction  is  an  inter¬ 
esting  return  to  quantifying  oxygen  saturation  using  susceptibility 
mapping  (de  Rochefort  et  al.,  2008;  Haacke  et  al.,  2010)  or  it  has 
also  been  called  susceptometry  (Fernandez-Seara  et  al.,  2006).  Per¬ 
haps  one  day  direct  measurements  of  changes  in  oxygen  saturation 
at  the  level  of  the  smallest  veins  will  be  a  viable  means  to  study 
local  neural  responses. 

Practically,  SWI  has  been  used  to  create  exquisite  images  of  the  veins 
without  contrast  agent.  The  standard  clinical  approaches  use  a  resolu¬ 
tion  of  0.5  x  1 .0  x  2  mm3  or  when  there  is  time  0.5  x  0.5  x  2  mm3.  The 
usual  TE,  TR  and  flip  angle  values  at  3  T  are  20  ms,  30  ms  and  1 5°.  Band¬ 
width  is  kept  low  at  around  100  Hz/pixel  to  keep  signal-to-noise  high. 
With  parallel  imaging,  whole  brain  coverage  for  the  lower  resolution 
of  0.5  x  1.0x2  mm3  takes  only  3  to  4  min.  An  example  of  a  processed 
MRA/SWI  data  set  at  3  T  is  shown  in  Fig.  6.  Until  recently,  only  one 
echo  was  used  with  SWI  scans,  but  there  are  now  a  number  of  papers 
extolling  the  virtues  of  multiple  echo  SWI  with  the  idea  that  the  first 
echo  can  be  used  for  MRA  (Barnes  and  Haacke,  2009b;  Deistung  et  al., 
2009;  Du  and  Jin,  2008;  Haacke  and  Reichenbach,  2011),  the  second 
for  SWI  and  the  combination  for  T2*.  Here  the  small  medullary  veins 
are  clearly  seen  in  the  MinIP  images  for  SWI.  With  sufficient  resolution, 
phase  can  be  measured  in  the  small  veins  and  is,  of  course,  proportional 
to  the  susceptibility,  which,  in  turn,  is  proportional  to  the  deoxyhemo¬ 
globin  content.  Phase  itself  has  been  used  to  measure  changes  in  flow 
indirectly  through  changes  in  oxygen  saturation  in  traumatic  brain  inju¬ 
ry  studies  (Shen  et  al.,  2007)  and  more  recently  in  stroke  (Li  et  al., 
2011b).  But  because  phase  is  non-local  and  depends  on  the  geometiy 
of  the  source  and  the  location  of  the  object  relative  to  the  main  field,  a 
direct  measure  of  the  venous  susceptibility  would  be  more  appropriate. 
Hence,  our  closing  section  will  be  about  the  use  of  susceptibility  map¬ 
ping  to  measure  oxygen  saturation  (Barnes  and  Haacke,  2009a;  Li  et 
al.,  2011a;  Reichenbach  et  al.,  1997). 

Susceptibility  weighted  imaging  and  mapping  (SWIM)  and  oxygen 
saturation  maps 

The  ability  to  map  oxygen  saturation  may  help  develop  the  potential 
of  high  resolution  SWI  as  a  means  to  detect  local  BOLD  changes,  especially 


a  b 


Fig.  6.  An  MRA  MaxIP  image  (a)  and  SWI  MinIP  image  (b)  from  the  same  data  set  with 
a  0.5  x  0.5  x  1  mm3  voxel  size.  The  original  data  was  acquired  using  a  single  echo  GRE 
sequence  with  TR/TE  =  30/20  ms.  Prior  to  the  SWI  MinIP  processing,  a  two-slice  sliding 
window  averaging  along  the  slice  direction  was  done  to  enhance  the  phase  effect. 
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at  high  fields.  The  presence  of  local  changes  in  susceptibility  causes  a  local 
change  in  magnetic  field  inside  the  object  and  a  non-local  change  outside 
the  object  of  interest.  The  phase  inside  the  object  is  a  constant  while  the 
change  outside  decays  as  1/r2  for  a  cylinder  and  1/r3  for  a  sphere.  These 
non-local  phases  are  generally  unwelcomed.  SWI  usually  emphasizes 
phase  of  one  sign  and  uses  this  to  enhance  the  magnitude  image 
(Reichenbach  et  al.,  1997;  Xu  and  Haacke,  2006). 

As  alluded  to  earlier,  the  phase  itself  is  of  great  interest.  It  has  been 
used  for  contrast  between  white  matter  and  gray  matter  (Haacke  et 
al„  1995,  1999;  Ogg  et  al.,  1999;  Shmueli  et  al.,  2009)  and  today  for 
contrast  between  a  variety  of  structures  in  the  brain  and  throughout 
the  body.  It  offers  the  potential  to  study  everything  from  vessel  wall 
to  even  imaging  nerves  (Haacke  et  al.,  1999).  The  problem  with 
phase  is  that  it  is  a  non-local  manifestation  of  the  magnetic  field. 
However,  this  very  non-local  behavior  contains  everything  we  need 
to  predict  the  inherent  local  susceptibility  source  (Deville  et  al., 
1979;  Haacke  et  al.,  2010;  Salomir  et  al.,  2003).  This  has  garnered 
great  interest  in  the  last  few  years  to  create  a  quantitative  susceptibil¬ 
ity  map  (QSM)  of  the  tissue  (de  Rochefort  et  al.,  2008;  Liu  et  al.,  2009; 
Marques  and  Bowtell,  2008;  Salomir  et  al.,  2003)  for  the  study  of 
brain  iron  and  oxygen  saturation.  The  latter  is  possible  because  the 
resulting  field  inside  the  veins  is  directly  related  to  the  hematocrit 
and  the  local  oxygen  saturation.  Therefore,  a  QSM  (or  as  we  refer  to 
it,  SWIM)  analysis  can  provide,  in  practice,  a  relative  oxygen  satura¬ 
tion  map  of  the  major  veins  in  the  body.  When  the  hematocrit  is 
known,  the  absolute  oxygen  saturation  can  be  found  (Barnes  and 
Haacke,  2009b;  Haacke  and  Reichenbach,  2011).  With  such  a  map, 
one  could,  in  principle,  open  a  new  means  to  study  fMRl  by  directly 
monitoring  changes  in  oxygen  saturation. 

As  a  means  to  demonstrate  a  systemic  BOLD  effect,  we  scanned  a 
subject  with  and  without  caffeine.  The  patient  ingested  a  200  mg 
NoDoze  pill  and  then  was  scanned  over  a  30  min  period.  The  SWI 
MinIP  images  are  shown  in  Figs.  7a  and  b,  and  the  SWIM  MaxIP  images 
are  shown  in  Figs.  7c  and  d,  both  before  and  after  caffeine  intake. 


Fig.  7.  Using  susceptibility  to  monitor  system  changes  in  oxygen  saturation.  The  SWI/ 
SWIM  results  before  (a,  c)  and  after  (b,  d)  200  mg  caffeine  intake  (one  NoDoz  pill) 
show  significant  systemic  susceptibility  changes  in  the  veins,  suggesting  a  correspond¬ 
ing  change  in  oxygenation  saturation.  In  the  future,  these  types  of  SWIM  images  may 
provide  a  source  of  monitoring  activation  in  high  resolution,  high  field  fMRl  studies. 


Although  the  SWI  projections  are  qualitatively  exquisite  images,  the 
QSM  or  SWIM  projections  are  quantitative  in  nature.  Despite  the  fact 
that  the  hematocrit  is  not  known,  a  measurement  of  the  venous  signal 
changes  normalized  to  the  pre-caffeine  data  provides  an  accurate  mea¬ 
sure  of  the  percent  change  in  oxygen  saturation  and  is  independent  of 
the  hematocrit.  It  is  very  difficult  to  measure  BOLD  changes  in  the  pa¬ 
renchyma  since  the  blood  volume  fraction  is  so  small  and  the  BOLD 
changes  are  small,  but  this  is  not  a  problem  in  the  major  veins 
(Sedlacik  et  al.,  2008).  Even  changes  in  the  medullary  veins  are  visible 
in  Figs.  7c  and  d.  Perhaps  this  new  approach  will  open  the  door  to 
using  the  veins  as  the  major  source  of  functional  activity  in  the  future 
rather  than  trying  to  discard  the  venous  information. 

Summary  and  conclusions 

With  the  advent  of  high  fields,  it  is  now  possible  to  image  the  vascu¬ 
lature  of  the  brain  with  high  resolution  rapidly.  Having  the  ability  do  to 
so  provide  a  means  to  understand  better  the  flow  and  function  of  the 
vasculature  system  and  to  model  the  vascular  system  for  any  individual. 
In  return,  modeling  the  system  can  predict  what  we  expect  to  see  from 
an  activated  to  resting  state  (Marques  and  Bowtell,  2008).  This  will  sure¬ 
ly  lead  to  a  better  understanding  of  why  different  individuals  respond 
differently  to  specific  stimuli,  and  certainly  to  separating  out  structural 
effects  from  real  functional  effects.  Along  these  lines,  we  would  recom¬ 
mend  that  all  fMRl  experiments  that  collect  T!  weighted  images  for  an¬ 
atomical  information  also  collect  a  high  resolution  SWI  scan  to  monitor 
the  venous  vasculature  in  the  activated  regions  of  interest. 

Further,  imaging  with  high  resolution  may  herald  the  day  when 
complicated  statistical  analysis  will  no  longer  be  needed  as  rapid  3D 
scanning  during  a  single  activated  and  resting  state  may  allow  for 
simple  subtractions  to  be  used  to  monitor  changes  in  the  vascular  sys¬ 
tem.  When  resolution  becomes  so  high  that  individual  local  veins 
draining  the  activated  region  can  be  discerned,  the  percent  change 
in  fMRl  experiments  can  jump  to  as  high  as  30%  (Cheng  and  Haacke, 
2001;  Haacke  et  al.,  1994;  Lai  et  al.,  1993). 
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Abstract 


Imaging  plays  a  key  role  in  the  diagnosis  and  longitudinal  follow  up  of  traumatic  brain  injury  (TBI). 
Among  injury  pathologies,  vascular  injury  is  associated  with  diffuse  axonal  injury  (DAI)  and  traumatic 


axonal  injury  (TAI).  The  vascular  network  is  ubiquitous  and  is  an  integral  part  of  the  tissue  structure. 

In  this  chapter,  we  focus  on  angiographic  and  venographic-related  imaging  methods  and  their  role  in 
assessing  mild,  moderate,  and  severe  TBI.  We  begin  with  an  introduction  to  susceptibility  weighted 
imaging  (SWI)  and  magnetic  resonance  angiography  (MRA)  and  then  provide  evidence  of  different  types 
of  vascular  damage.  Examples  of  TBI-induced  microbleeds  are  presented  along  with  the  concept  of  low- 
impact  medullary  vein  damage  (MVD).  This  MVD  has  been  seen  even  for  so-called  mild  TBI  cases. 
Vascular  damage  can  also  manifest  as  a  reduction  in  local  perfusion  even  when  no  clear  macroscopic 
vessel  damage  is  seen.  To  further  understand  the  role  of  vascular  abnormalities,  we  then  introduce  the 
different  perfusion  weighted  imaging  (PWI)  techniques  available  and  their  application  in  TBI.  The 
combination  of  SWI  and  PWI  should  make  it  possible  to  differentiate  the  role  of  local  thrombus  versus 
changes  in  oxygen  saturation  in  MVD.  for  example.  Since  MRA  and  SWI  are  able  to  provide  a  full 
description  of  the  brain’s  vasculature  in  3D,  we  briefly  discuss  the  presence  of  finite  element  modeling  in 
understanding  vascular  injury.  We  conclude  with  recommendations  related  to  the  use  of  perfusion  with 
MRA,  SWI,  and  oxygen  saturation  measurements  to  obtain  a  complete  picture  of  the  hemodynamics  of  the 
brain. 
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The  Presence  of  Venous  Damage  and  2 

Microbleeds  in  Traumatic  Brain  Injury  3 

and  the  Potential  Future  Role  of  Angiographic  4 

and  Perfusion  Magnetic  Resonance  Imaging  5 

E.M.  Haacke,  Waqar  Raza,  Wu  Bo,  and  Zhifeng  Kou  6 

4.1  Introduction  7 

Imaging  is  playing  a  more  and  more  important  role  in  the  diagnosis  and  longitudinal  8 
follow  up  of  traumatic  brain  injury  (TBI).  It  has  been  suggested  that  for  both  9 
anatomical  and  functional  purposes  magnetic  resonance  imaging  (MRI)  offers  one  10 
of  the  best  means  to  investigate  the  status  of  brain  function  and  hemodynamics  11 
(Duhaime  et  al.  2010;  Haacke  et  al.  2010a).  Among  injury  pathologies,  vascular  12 
injury  is  an  important  part  of  the  spectrum.  It  is  associated  with  diffuse  axonal  13 
injury  (DAI)  and  traumatic  axonal  injury  (TAI).  This  bleeding  has  become  a  sym-  14 
bolic  sign  of  brain  tissue  damage  from  DAI  or  TAI  in  diagnostic  radiology  (Grossman  1 5 
1996).  The  vascular  network  is  ubiquitous  and  is  present  at  all  levels  in  terms  of  size  16 
from  the  5  pm  capillaries,  to  the  50  pm  arterioles  and  venules,  to  the  several  100  pm  17 
large  vessels,  and  finally  those  vessels  on  the  order  of  a  millimeter.  Presumably,  18 
disruption  of  these  vessels  leads  to  a  disturbance  of  local  blood  perfusion.  19 

Vascular  damage  can  manifest  in  a  number  of  ways.  Recently,  Kou  et  al.  gave  a  20 
summary  of  the  different  patterns  of  vascular  bleeding  (Kou  et  al.  2008,  2010).  21 

Although  more  dated,  chapter  by  Grossman  et  al.  gives  a  classical  description  of  the  22 
different  types  of  hemorrhages  (Grossman  1996).  Anatomically,  at  a  macroscopic  23 
level,  both  hemorrhages  and  non-hemorrhagic  shearing  lesions  are  found  in  major  24 
white  matter  (WM)  tracts,  including  frontal  white  matter,  corpus  callosum,  internal  25 
capsule,  the  upper  part  of  the  brain  stem,  the  gray  matter  (GM)/(WM)  boundaries  26 
and  the  cerebellar  region.  Morphologically,  there  are  several  types  of  bleeding.  27 
Extra  axially,  there  are  epidural  and  subdural  hematomas,  which  could  be  caused  by  28 
direct  impact;  subarachnoid  hemorrhage,  which  can  be  seen  reasonably  well  visual-  29 
ized  with  computed  tomography  (CT)  and  now  with  MRI  as  well  (Wu  et  al.  2010)  30 
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31  and  periventricular  hemorrhage,  which  is  due  to  the  vascular  damage  along  the 

32  ventricle  walls  or  associated  with  damage  in  the  splenium  of  the  corpus  callosum. 

33  For  parenchymal  hemorrhage,  examples  include  (a)  petechial  microbleeds,  where 

34  blood  is  seen  to  build  up  at  isolated  points  after  TBI;  (b)  DAI,  where  multi-focal 

35  bleeding  is  seen  in  the  major  white  matter  tracts;  and  (c)  focal  contusions,  including 

36  coup  and  contra-coup  contusions,  which  are  often  seen  in  direct  impact  or  falls. 

37  If  the  focal  contusion  reaches  a  certain  size,  it  may  need  immediate  surgical  inter- 

38  vention  to  alleviate  the  rise  of  intracranial  pressure  (ICP). 

39  Of  course  any  impact  may  have  microvascular  repercussions  and,  as  such,  damage 

40  could  easily  occur  at  both  axonal  and  microvascular  levels  simultaneously  since 

41  axons  and  neurons  are  completely  intertwined  with  the  vasculature  (although  each 

42  has  its  own  unique  mechanical  properties).  From  a  biomechanical  point  of  view,  the 

43  neural  vasculature  provides  physical  and  structural  support  of  neuronal  and  glial 

44  cells.  From  a  physiological  point  of  view,  the  blood  flow  provides  nutritional  sup- 

45  port  for  neuronal  and  glial  cells.  By  the  time  of  axonal  injury,  it  is  hard  to  imagine 

46  that  there  is  pure  axonal  injury  without  any  vascular  damage  or  functional  distur- 

47  bance.  The  impairment  of  the  neurovascular  blood  supply  further  jeopardizes  the 

48  nutritional  support  and  the  physiological  hemostasis  environment  for  neuronal  cells 

49  to  survive  the  recovery  battle.  This  impairment  may  be  manifest  as  a  reduction  in 

50  perfusion  to  the  tissue.  An  immediate  perfusion  reduction  is  seen  when  using  the 

51  weight  drop  impact  model  in  rats  using  the  Marmarou  model  even  without  any  evi- 

52  dent  tissue  damage  (Shen  et  al.  2007).  Perfusion  usually  recovers  to  normal  levels 

53  within  72  h  in  rats.  Depending  on  the  injury  severity,  the  vascular  damage  could 

54  present  as  macroscopic  damage  or  bleeding  or  severe  reductions  in  deoxyhemoglo- 

55  bin  without  any  direct  evidence  of  tissue  damage.  Damage  to  the  vessel  could  lead 

56  to  an  impaired  neural  vascular  tone,  either  sympathetic  or  parasympathetic  vasocon- 

57  striction,  and  further  result  in  the  loss  of  control  of  blood  flow  or  to  reduced  perfu- 

58  sion  (Bouma  and  Muizelaar  1995).  Or  it  may  even  be  an  autoimmune  response 

59  through  which  the  brain  protects  itself  in  case  there  is  major  bleeding  after  impact. 

60  We  just  do  not  yet  understand  the  mechanisms  behind  these  immediate  responses  to 

61  mechanical  impact. 

62  Furthermore,  arteries  and  veins  are  different  in  many  aspects.  Anatomically, 

63  major  veins  are  located  on  the  surface  of  the  brain  with  branches  collecting  return- 

64  ing  blood  from  the  inside  of  the  brain,  while  major  arteries  are  located  inside  the 

65  brain  and  tend  to  branch  out.  From  Ommaya’s  Centripetal  Theory  or  Gennarelli’s 

66  Depth  of  Lesion  model,  with  the  increase  of  the  injury  severity,  the  injured  locations 

67  tend  to  penetrate  from  the  surface  deeper  into  the  brain.  That  means  the  brain  sur- 

68  face  is  always  the  most  vulnerable  location  of  injury,  which  makes  veins  exposed  to 

69  this  trauma  more  vulnerable.  Biomechanically,  arteries  are  smaller  but  much  stron- 

70  ger  than  veins  since  they  must  withstand  the  systolic  blood  pressure  during  the 

71  pumping  of  blood.  Weaker  venous  walls  makes  veins  more  vulnerable  than  arteries 

72  for  the  same  force  of  impact.  As  a  result,  the  veins  on  the  cortical  surface  tend  to  be 

73  the  most  vulnerable  part  of  the  vascular  structure.  Damage  here  may  also  lead  to 

74  deeper  damage  to  the  veins  draining  into  these  surface  veins.  Originally  developed 

75  as  a  venography  technique,  susceptibility  weighted  imaging  (SWI)  has  the  unique 
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capacity  to  visualize  both  venous  structures  and  hemorrhagic  bleeds.  This  allows 
researchers  to  visualize  the  spatial  relationship  between  bleeding  and  vascular  struc¬ 
tures.  Therefore,  it  is  possibly  to  identify  the  possible  bleeding  of  veins  and/or  the 
changes  in  local  oxygen  saturation  of  the  tissue  feeding  these  veins. 

Modern  imaging  technologies  now  make  it  possible  to  image  the  vasculature, 
perfusion,  and  oxygen  saturation  all  in  one  sitting  with  MRI  methods.  Although  the 
former  can  also  be  done  with  CT  angiography  (CTA)  using  a  contrast  agent,  CT  is 
based  on  the  use  of  ionizing  radiation  in  the  form  of  x-rays.  On  the  other  hand,  MR 
angiography  (MRA)  and  SWI  can  be  used  to  image  the  arteries  and  the  veins, 
respectively,  down  to  250  pm  without  a  contrast  agent.  Perfusion  imaging  can  also 
be  done  with  CT  with  the  help  of  a  contrast  agent.  MRI  also  requires  a  contrast 
agent  for  dynamic  contrast  enhanced  (DSC)  perfusion  weighted  imaging  (PWI). 
However  a  contrast  agent  is  not  required  for  arterial  spin  labeling  (ASL)  or  for  an 
average  cerebral  blood  flow  measurement  using  flow  quantification  in  the  neck 
(Haacke  2012).  Although  oxygen  saturation  can  be  measured  using  positron  emis¬ 
sion  tomography,  it  is  much  more  accessible  and  patient  friendly  to  use  when 
extracted  from  MRI  measures  such  as  T2(Lu  and  Ge  2008;  Ge  et  al.  2012),  T2*(Li 
et  al.  1998),  or  susceptibility  mapping  (Haacke  et  al.  2010b). 

In  this  chapter,  we  will  focus  on  these  angiographic-  and  venographic-related 
imaging  methods  and  their  role  in  assessing  mild,  moderate,  and  severe  TBI.  We 
will  begin  with  an  introduction  to  SWI  and  MRA  and  then  provide  evidence  of  dif¬ 
ferent  types  of  vascular  damage  at  different  levels  of  impact.  We  will  show  that  for 
mild  TBI,  even  when  other  damage  is  not  visible  with  most  clinical  techniques, 
there  is  potential  venous  damage.  We  will  then  introduce  the  different  perfusion 
techniques  available  and  their  application  in  TBI.  After  that,  we  will  discuss  the  role 
of  finite  element  modeling  in  understanding  the  vascular  injury,  and,  finally  we  will 
conclude  with  recommendations  related  to  the  use  of  perfusion  with  MRA,  SWI, 
and  oxygen  saturation  measurements.  Taken  together,  these  MRI  techniques  offer 
different  aspects  of  neural  vascular  injury  and  tend  to  complement  each  other  to 
depict  a  clear  panorama  of  this  phenomenon. 


4.2  Susceptibility  Weighted  Imaging  (SWI) 
and  Mapping  (SWIM) 

SWI  exploits  the  susceptibility  differences  between  tissues.  It  uses  phase  images  to 
detect  these  differences  and  subsequently  a  phase  mask  to  enhance  the  presence  of 
veins,  blood  products,  and  iron-containing  tissues.  Due  to  its  sensitivity  to  deoxyhe¬ 
moglobin  in  venous  blood  and  hemosiderin  in  hemorrhage,  SWI  is  commonly  used 
in  TBI  studies  (Haacke  et  al.  2010b;  Tong  et  al.  2008).  Compared  with  the  tradi¬ 
tional  invasive  manner  of  clinical  monitoring  of  cerebral  vascular  damage  and 
reductions  in  blood  flow,  this  method  offers  a  novel,  safe,  and  noninvasive  approach 
to  quantify  changes  in  oxygen  saturation  and  cerebral  blood  flow  (Shen  et  al.  2007; 
Fujima  et  al.  2011)  and  to  visualize  structural  changes  in  the  brain’s  vasculature 
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117  after  TBI.  These  three  pieces  of  information  can  be  used  to  improve  the  diagnosis 

118  of  the  state  of  health  of  brain  tissue  both  pre-  and  post-TBI.  The  detection  of  micro- 

1 1 9  hemorrhages  and  shearing  lesions  in  TAI/DAI  in  trauma  patients  is  often  difficult  as 

120  the  injuries  tend  to  be  relatively  small  and  can  be  easily  missed  by  low-resolution 

121  scans.  SWI  is  usually  run  at  relatively  high  resolution  (1  mm3  or  better)  and  is 

122  extremely  sensitive  to  bleeding  in  the  GMAVM  boundaries  making  it  possible  to  see 

123  very  small  lesions  and  increasing  the  ability  to  detect  more  subtle  injuries. 

124  The  phase  information  collected  in  SWI  represents  the  magnetic  field  distribu- 

125  tion  created  from  the  different  sources  such  as  veins,  microbleeds,  ferritin,  and 

126  hemosiderin.  It  is  possible  to  use  this  phase  information  to  find  the  local  source  of 

127  the  field  changes.  Using  a  rather  simple  dipole  approximation,  it  is  possible  to 

128  develop  an  inverse  reconstruction  approach  that  uses  the  phase  information  to  create 

129  a  source  or  susceptibility  map  (Haacke  et  al.  2010b;  Zaitsu  et  al.  201 1).  We  refer  to 

130  this  quantitative  susceptibility  mapping  as  SWIM.  The  problem  with  phase  is  that  a 

131  simple  spherical  dipole  source  will  create  a  non-local  phase  behavior  that  will  make 

132  it  very  difficult  to  do  a  good  job  with  SWI  by  itself.  By  creating  SWIM  data,  that 

133  phase  disappears  and  is  replaced  by  the  source  itself,  the  local  sphere.  This  new 

1 34  approach  should  make  it  possible  to  map  out  iron  content  in  microbleeds  and  to  map 

135  oxygen  saturation  in  the  major  veins  in  the  brain  (Haacke  et  al.  2010b,  2011). 

136  4.2.1  Visible  Vascular  Injury  in  TBI  and  the  Presence 

137  of  Damaged  or  Compromised  Medullary  Veins 

138  In  this  section,  we  will  review  what  can  be  seen  in  TBI  using  SWI  and  discuss  a  new 

139  finding,  the  presence  of  damaged  or  compromised  medullary  veins.  Whenever  there 

140  is  damage  to  the  vascular  system  that  leaves  behind  extravasated  blood,  there  is  a 

141  good  chance  that  it  will  be  seen  with  SWI.  Depending  on  the  scope  of  the  damage, 

142  one  may  see  local  damage  as  focal  lesions  or  one  may  see  widespread  damage  in  the 

143  form  of  DAI.  This  damage  may  or  may  not  produce  an  alteration  in  brain  function. 

144  Successfully  visualizing  the  damaged  veins  by  itself  is  an  important  biomarker  of 

145  TBI  but,  to  be  sure  if  there  is  or  is  not  a  perfusion  deficit,  one  also  needs  to  measure 

146  perfusion.  Even  if  there  is  a  perfusion  loss  in  the  putatively  damaged  tissue,  one 

147  needs  to  measure  oxygen  saturation  to  see  if  the  tissue  function  is  intact.  With  all 

148  four  components,  the  anatomy  of  the  vasculature,  the  presence  of  vascular  damage, 

149  the  state  of  perfusion,  and  the  level  of  oxygen  saturation,  it  should  be  possible  to 

150  obtain  a  relatively  complete  picture  of  the  hemodynamics  of  the  brain. 

151  Asa  demonstration  of  the  types  of  venous  damage  we  are  now  seeing  with  SWI, 

152  we  show  several  cases  of  medullary  vein  damage  (MVD)  from  frontal  and  side 

153  impacts.  Figure  4.1  shows  two  separate  cases  of  mild  TBI  with  clear  damage  to  the 

154  medullary  veins.  In  reviewing  more  than  100  TBI  cases,  including  mild,  moderate, 

155  and  severe,  we  have  seen  35  cases  with  MVD.  Of  the  100  cases,  15  of  these  were 

156  mild  and  of  these  15  only  3  showed  MVD.  However,  in  these  three  cases,  no  other 

157  technique  showed  the  presence  of  tissue  damage  (Figs.  4.1,  4.2,  and  4.3). 
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Fig.  4.1  Medullary  vein  damage  in  the  frontal  lobes  shown  with  SWI.  Resolution  is 
0.5  mm  x  0.5  mmx2  mm.  (a)  MVD  emanating  from  the  septal  vein  in  the  left  frontal  lobe  ( upper 
arrows).  This  image  also  shows  blood  in  the  ventricles  ( lower  arrows),  (b)  MVD  in  both  the  left 
and  right  frontal  lobes.  The  damaged  medullary  veins  are  depicted  by  arrows 


Fig.  4.2  Lateral  medullary  vein  damage  shown  with  SWI.  Resolution  is  0.5  mm  x  0.5  mm  x  2  mm. 
(a)  Lateral  MVD  in  the  frontal  lobes,  (b)  Lateral  MVD  in  the  left  temporal  lobe.  This  patient  also 
shows  damage  in  the  right  anterior  cerebral  artery  territory,  (c)  Lateral  MVD  in  the  left  frontal  lobe 
emanating  from  the  septal  vein.  The  damaged  medullary  veins  are  depicted  by  arrows 


4.3  MR  Angiography  isa 

Depending  on  the  impact  strength,  the  vasculature  can  be  severely  damaged  in  159 
TBI.  It  behooves  us  to  study  the  small  vessel  damage  and  eventually  compare  this  160 
to  any  perfusion  deficit.  The  SWI  data  above  already  suggest  that  veins  are  more  161 
sensitive  to  abnormal  stress/strain  changes  than  arteries  because  they  have  weaker  162 
vessel  walls.  However,  one  would  expect  that  would  also  be  true  for  smaller  arter-  163 
ies.  In  any  case,  to  really  do  a  complete  job  studying  TBI,  imaging  both  the  small  164 
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Fig.  4.3  Lateral  medullary  vein  damage  shown  with  both  SWI  and  SWIM.  Resolution  is 
0.5  minx 0.5  mmx2  mm.  The  same  subject  from  Fig.  4.2c  is  shown  here  with  (a)  the  usual  SWI 
minimum  intensity  projection  and  (b)  the  processed  phase  data  used  to  create  a  SWIM  image.  The 
SWIM  data  makes  it  possible  to  quantify  the  amount  of  iron  in  the  bleed  and  to  monitor  if  this  is 
resorbing  over  time,  remains  stable  or  continues  to  bleed.  The  damaged  medullary  veins  are 
depicted  by  long  arrows  and  the  microbleeds  by  short  arrows 


165  arteries  and  small  veins  would  seem  to  be  the  right  way  to  go.  Up  to  recently,  the 

166  major  focus  in  imaging  vessels  has  been  at  the  M1/M2/M3  levels  but  not  beyond. 

167  Now  the  capability  exists  to  see  much  smaller  arteries,  on  the  order  of  100-250  pm. 

168  The  methods  that  make  this  possible  are  the  usual  time-of-flight  (TOF)  MRA 

169  methods,  TOF  with  background  suppression  such  as  magnetization  transfer  con- 

170  trast  (MTC)  or  preferably,  at  high  field,  using  a  rephasing/dephasing  method,  and 

171  finally  the  use  of  a  contrast  agent  to  enhance  the  signal  from  the  vessels  far  above 

172  the  background  tissue.  In  principle,  the  latter  two  can  also  be  combined  for  poten- 

173  tial  further  improvement  in  vessel-tissue  contrast.  An  example  of  a  conventional 

174  MR  angiogram  and  the  same  individual  post-contrast  is  shown  in  Fig.  4.4a.  The 

175  data  are  collected  with  0.5  mm  isotropic  resolution  which  means  that  even  some 

176  2  50  pm-sized  vessels  may  be  visualized  if  their  signal  is  high  enough.  (Pushing 

177  this  to  250  pm  in  plane  resolution  makes  it  possible  to  see  even  100  pm-sized  arter- 

178  ies  but  takes  too  long  practically  at  the  moment.)  (Haacke  et  al.  201 1).  The  advan- 

179  tage  of  the  conventional  approach  is  that  there  are  no  veins  visible  in  the  images. 

180  When  a  contrast  agent  is  used,  both  arteries  and  veins  get  bright,  although  the 

181  larger  veins  tend  to  be  not  as  bright  at  3  T  as  they  have  a  much  lower  T2*.  Despite 

182  the  veins  causing  some  difficulty  in  visualizing  the  arteries,  with  high-resolution 

183  3D  imaging,  arteries  and  veins  can  be  separated  by  image  processing.  Figure  4.4b 

184  shows  that  there  are  indeed  smaller  arteries  now  visible  that  could  not  be  seen 

185  without  the  use  of  a  contrast  agent  (Fig.  4.4). 
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Fig.  4.4  MR  angiograms  with  0.5  mm  isotropic  resolution,  (a)  A  conventional  TOF  MRA.  (b)  A 
contrast-enhanced  TOF  MRA.  Clearly  many  veins  are  also  now  seen  in  the  contrast-enhanced 
image.  However,  there  are  also  many  small  arteries  now  visible  that  could  not  be  seen  before 


4.4  Clinical  Impact  of  Vascular  Damage  and  Disrupted 
Perfusion  to  the  Brain 

Damage  to  vessels  may  cause  a  change  in  the  brain’s  microcirculation  that  leads  to 
a  persistent  state  of  hypoperfusion  and  improper  delivery  of  vital  metabolites  to  the 
brain  parenchyma.  The  toxic  contents  may  also  leak  through  damaged  blood  brain 
barrier  (BBB)  into  the  already  vulnerable  neurons.  This  is  not  dissimilar  to  what  we 
see  with  stroke  patients,  and  it  would  behoove  us  to  collect  the  same  type  of  infor¬ 
mation  for  TBI  patients.  Early  impairment  of  cerebral  blood  flow  (CBF)  in  patients 
with  severe  injury  correlates  with  poor  brain  tissue  oxygenation  (Zwienenberg  and 
Muizelaar  2001).  Some  researchers  have  suggested  that  it  would  be  important  to 
normalize  cerebral  perfusion  and  blood  flow  after  TBI  for  better  outcome  (Gennarelli 
1987).  In  clinical  studies,  a  global  reduction  in  CBF  has  been  found  in  the  first  few 
hours  following  TBI  using  133Xe  (Zwienenberg  and  Muizelaar  2001;  Gennarelli 
1987;  Gennarelli  and  Graham  1998;  Kelly  and  Martin  1997).  This  is  similar  to  what 
is  found  in  animal  studies  as  well  (DeWitt  et  al.  1992).  Patients  who  have  a  global 
CBF  below  the  threshold  for  infarction  (approximately  18  mF/100  g/min)  (Bouma 
et  al.  1991)  have  a  worse  clinical  outcome.  However,  Zhao  et  al.  reported  a  thresh¬ 
old  for  ischemia  as  30  mL/100  g/min  in  his  study  (Zhao  and  Nowak  2006).  Kelly 
et  al.  showed  that  outcome  for  patients  with  CBF  less  than  this  were  poor  (Kelly  and 
Martin  1997).  CBF  and  brain  tissue  oxygenation  are  proving  to  be  crucial  parame¬ 
ters  to  monitor  in  intensive  care  units  along  with  others.  Moreover,  measurement  of 
absolute  regional  blood  flow  at  rest  allows  a  clear  estimation  of  the  magnitude  of 
blood  flow  changes  that  occur  with  drug  treatment  or  other  interventions.  There  are 
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209  several  perfusion-imaging  techniques  available  today  for  measuring  cerebral 

210  perfusion  in  brain  tissue  and  for  visualizing  the  healthy  and  damaged  vessels  in  the 

21 1  brain.  We  give  a  brief  overview  of  these  as  they  are  related  to  applications  in  TBI. 


212  4.4.1  Computer  Tomography 

213  Since  the  prohibition  of  the  use  of  Xenon  by  the  FDA,  contrast  enhanced  perfusion 

214  has  become  the  method  of  choice  in  CT.  Computed  tomography  perfusion  (CTP)  is 

215  initiated  by  injecting  iodinated  contrast  media  (iodixanol)  into  a  peripheral  vein. 

216  Only  a  few  studies  have  been  performed  on  TBI  using  this  technique  but  the  results 

217  show  that  there  were  significant  perfusion  changes.  One  study  reported  a  significant 

218  decrease  of  blood  flow  after  trauma  in  a  group  comparison  (Metting  et  al.  2010). 

219  However,  unlike  conventional  CT,  which  can  show  bleeding  or  mass  effect  on  indi- 

220  vidual  basis,  perfusion  changes  on  CTP  would  rely  on  quantitative  measurement  to 

221  show  group  differences.  Wintermark  et  al.  suggested  that  intermittent  perfusion-CT 

222  measurements  will  complement  the  information  acquired  by  continuous  cerebral 

223  perfusion  pressure  (CPP)  alone.  Perfusion-CT  gives  exclusive  information  about 

224  the  regional  heterogeneity  of  brain  perfusion.  This  information  will  be  very  useful 

225  in  differentiating  patients  with  preserved  auto-regulation  and  impaired  autoregula- 

226  tion  and  could  therefore  steer  interpretation  of  CPP  measurements  and  therapy 

227  (Wintermark  et  al.  2004). 


228  4.4.2  Positron  Emission  Computer  Tomography 

229  Positron  emission  tomography  (PET)  is  a  nuclear  medicine  test  in  which  tissue 

230  function  can  be  imaged.  In  A  PET  scan,  the  tracer  stays  in  the  blood  stream  rather 

231  than  being  absorbed  by  surrounding  tissues,  thereby  limiting  the  images  to  areas 

232  where  blood  flows.  Damaged  tissues  have  reduced  metabolic  activity;  therefore, 

233  gamma  radiation  from  these  areas  is  reduced  or  absent  (Alexander  1995).  Indeed, 

234  some  PET  data  suggest  discordance  between  blood  flow  and  metabolism  following 

235  TBI  (Obrist  et  al.  1984);  (Abu-Iudeh  et  al.  1998;  Mintun  et  al.  2001).  The  results  of 

236  the  Mintum  study  indicate  that  the  metabolic  state  of  the  traumatically  injured  brain 

237  should  be  defined  differentially  in  terms  of  glucose  and  oxygen  metabolism.  The 

238  use  of  FDG-PET  demonstrates  that  hyperglycolysis  occurs  both  locally  and  glob- 

239  ally  following  severe  head  injury  in  humans  cellular  pathophysiology  characteristic 

240  of  head  injury  (Bergsneider  et  al.  1997).  TBI  studies  in  animals  have  shown  a  well- 

241  defined  triphasic  pattern  of  the  cerebral  metabolic  rate  of  glucose  (CMRglc) 

242  (Bergsneider  et  al.  2001).  The  initial  brief  response  of  hyperglycolysis  is  followed 

243  by  a  second  relatively  prolonged  period  of  metabolic  depression.  The  second  period, 

244  characterized  by  metabolic  depression,  typically  lasts  several  hours  and  is  associ- 

245  ated  with  persistent  neurological  deficits  as  long  as  the  metabolic  depression  is 
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present.  The  third  period  is  characterized  by  a  recovery  behavioral  function  that 
parallels  that  of  the  recovery  of  metabolic  function.  No  such  pattern  was  seen  in 
human  data  (Bergsneider  et  al.  2001). 

Soustiel  et  al.  compared  the  respective  effects  of  established  measures  used  for 
management  of  TBI  patients  on  CBF  and  cerebral  metabolic  rates  of  oxygen 
(CMR02),  glucose  (CMRGlc),  and  lactate  (CMRLct)  and  found  out  that  moderate 
hyperventilation  could  aggravate  pre-existing  impairment  of  CBF  and  metabolism 
(Soustiel  et  al.  2006).  Therefore,  TBI  patients  should  be  carefully  monitored.  Some 
suggest  the  use  of  mannitol  for  ICP  control. 

It  is  now  well  established  that  PET  is  somewhat  superior  to  SPECT  in  terms  of 
spatial  resolution  and  permits  a  better  visualization  of  subtle  abnormalities.  With 
the  advent  and  implementation  of  novel  tracers  as  well  as  improvement  in  instru¬ 
mentation,  molecular  brain  imaging  using  PET  may  be  able  to  illustrate  TBI 
pathophysiology  and  point  to  novel  treatment  strategies.  However,  despite  these 
advantages  PET  requires  the  use  of  radiation. 


4.4.3  Single-Photon  Emission  Computed  Tomography  (SPECT) 

SPECT  is  a  nuclear  medicine  tomography  imaging  technique  that  uses  gamma  rays 
to  probe  the  tissue,  specifically  to  show  blood  flow  to  tissues  and  organs.  A  SPECT 
scan  integrates  two  technologies  to  view  the  body  namely,  CT  and  a  radioactive 
material  (tracer).  The  tracer  makes  it  possible  to  see  how  blood  flows  to  tissues  and 
organs.  Before  the  SPECT  scan,  a  gamma-emitting  radio-labeled  chemical  is 
injected  into  the  blood  stream  that  can  then  be  detected  by  the  scanner.  The  com¬ 
puter  collects  the  information  emitted  by  the  gamma  rays  and  translates  them  into 
two-dimensional  cross-sections.  These  cross-sections  can  be  merged  to  create  a  3D 
image  of  the  brain.  A  SPECT  scan  is  primarily  used  to  view  blood  flow  through  both 
arteries  and  veins  in  the  brain.  Tests  have  shown  that  it  might  be  more  sensitive  to 
brain  injury  than  either  MRI  or  CT  scanning  because  it  can  detect  reduced  blood 
flow  to  injured  sites  (Lewine  et  al.  2007;  Abdel-Dayem  et  al.  1998).  These 
qualifications  aside,  SPECT  offers  promise  with  regard  to  meeting  the  criteria  of 
sensitivity,  incremental  validity,  and  clinical  correlates.  SPECT  scans  are  cheaper 
and  more  readily  available  than  the  higher  resolution  PET  scans,  but  they  also  still 
involve  the  use  of  radiation  and  radioactive  agents.  Another  reason  for  SPECT  not 
gaining  so  much  popularity  is  that  it  requires  co-registration  with  MRI. 


4.4.4  Laser  Doppler  Flow 

Laser  Doppler  is  an  inexpensive,  noninvasive  method  of  measuring  the  continu¬ 
ous  circulation  of  blood  flow  on  a  microscopic  level.  The  beam  from  a  low-power 
laser  (usually  a  laser  diode)  penetrates  the  skin  sufficiently  to  be  scattered  with  a 
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283  Doppler  shift  by  the  red  blood  cells  and  returns  to  be  concentrated  on  a  detector. 

284  These  measurements  are  useful  to  monitor  the  effect  of  exercise,  drug  treatments, 

285  environmental,  or  physical  manipulations  on  targeted  micro-sized  vascular  areas. 

286  It  is  reported  by  Rafols  et  al.  using  Laser  Doppler  Flow  (LDF)  that  the  hypoper- 

287  fusion  following  TBI  is  a  result  of  the  loss  of  autoregulatory  response  of  reacting 

288  vessels  of  the  brain.  This  loss  of  autoregulation  causes  the  improper  delivery  of 

289  oxygen  and  metabolites  to  the  brain  tissue  and  results  in  secondary  injury  and  poor 

290  cognitive  outcome  (Rafols  et  al.  2007).  In  this  study,  he  used  the  Marmarou  model 

291  to  induce  TBI  to  avoid  focal  brain  contusions  and  tissue  necrosis.  This  study  revealed 

292  a  37%  reduction  in  cortical  blood  flow.  Ultrastructural  alterations  in  the  lumen  and 

293  wall  of  sensorimotor  cortical  microvessels,  including  endothelial  cell  distortion  and 

294  luminal  collapse,  were  seen  between  1  and  48  h  after  trauma.  Smooth  muscle  in  the 

295  wall  of  reacting  microvessels  showed  evidence  of  increased  contractility  that  coin- 

296  cided  temporally  with  the  decreased  perfusion  of  cortical  CBF. 

297  Thomale  et  al.  showed  that  the  acute  hypoperfusion  after  TBI  is  due  to  vasocon- 

298  striction  (Thomale  et  al.  2002).  This  study  also  found  hyperemia  occurring  between 

299  24  and  48  h  after  TBI  in  the  controlled  cortical  impact  model  using  LDF.  Muir  et  al. 

300  studied  the  continuous  pattern  of  CBF  following  TBI  in  a  lateral  percussion  fluid 

301  model  (Muir  et  al.  1992).  The  results  showed  that  seconds  after  injury,  blood  flow  in 

302  the  left  parietal  cortex  increased  about  200%  followed  falling  to  about  80%  of  control 

303  after  5  min,  and  to  56%  at  the  1  h  time  point.  Dewitt  et  al.  reported  that  L-arginine 

304  administered  after  TBI  prevents  post-traumatic  hypoperfusion  and  that  pre-treatment 

305  with  superoxide  dismutase  will  restore  CBF  after  a  brief  period  of  hypoperfusion 

306  (DeWitt  et  al.  1992).  Lam  et  al.  performed  a  study  to  monitor  autoregulation  using 

307  LDF  in  patients  with  head  injury.  It  was  concluded  from  the  study  that  maintenance  of 

308  CPP  is  important  in  preventing  cerebral  ischemia  (Lam  et  al.  1997).  Zweifel  et  al. 

309  showed  that  changes  in  CBF  detected  by  LDF  matched  well  with  flow  measurements 

310  performed  with  Transcranial  Doppler  Sonography  (Christian  et  al.  2010). 

31 1  The  clinical  application  of  LDF  is  limited  by  a  number  of  problems.  It  can  only 

312  access  a  small  volume  of  brain  tissue  and  thus  can  only  monitor  a  focal  area  and  LDF 

313  does  not  provide  absolute  values  for  blood  flow.  Further,  since  LDF  is  operator  depen- 

314  dent  especially  in  making  contact  on  the  skull  for  cases  of  subdural  hematoma  or  prom- 

31 5  inent  subarachnoid  space,  inaccurate  measurements  may  follow  (Lam  et  al.  1997). 


316  4.4.5  Transcranial  Doppler  Sonography 

31 7  Transcranial  Doppler  Ultrasound  (TCD)  is  a  noninvasive  method  for  the  measurement 

318  of  flow  in  the  major  vessels  in  the  head  such  as  the  middle  cerebral  artery  (MCA), 

319  petrosal  sinuses,  and  ophthalmic  vessels.  It  is  commonly  used  in  the  standard  care  of 

320  TBI  patients.  TCD  is  particularly  accurate  in  detection  of  episodes  of  hypoperfusion 

321  induced  by  low  CPP  (Brandi  et  al.  2010),  (Schmidt  et  al.  2000)  and  ICP  (Homburg 

322  et  al.  1993).  This  is  a  relatively  quick  and  inexpensive  test  and  is  growing  in  popularity 

323  in  the  USA.  The  equipment  is  portable,  making  it  possible  to  use  not  only  in  the 
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hospital  but  also  in  a  doctor’s  office,  nursing  home,  or  patient’s  home.  There  are  two 
methods  used  for  recording  the  data.  The  first  method  uses  “B-mode”  imaging,  which 
displays  a  two-dimensional  image  as  seen  by  the  ultrasound  probe.  Once  the  desired 
blood  vessel  is  found,  blood  flow  velocities  can  be  measured  with  a  pulsed  Doppler 
probe  to  plot  velocity  over  time.  Together,  these  make  a  duplex  test.  The  second 
method  of  recording  uses  only  the  second  probe  function,  relying  instead  on  the 
training  and  experience  of  the  clinician  in  finding  the  correct  vessels. 

A  number  of  observational  studies  have  suggested  that  vasospasm  develops  in 
some  patients  following  TBI  (Oertel  et  al.  2005).  Lee  et  al.  confirmed  that  hemody- 
namically  significant  vasospasm  with  low  CBF  was  a  significant  predictor  of  poor 
outcome  (Lee  et  al.  1997).  These  studies  have  almost  exclusively  used  TCD  to  iden¬ 
tify  vasospasm  criteria  similar  to  those  used  in  subarachnoid  hemorrhage.  The  limi¬ 
tation  of  transcranial  Doppler  is  penetrating  the  bone.  There  are  several  “windows” 
through  the  temporal  area  that  allow  the  visualization  of  the  petrosal  sinuses  and  the 
ophthalmic  vein.  TCD  can  also  be  used  to  measure  flow  in  the  middle,  anterior,  and 
posterior  cerebral  arteries. 

Cerebral  perfusion  should  be  estimated  as  soon  as  possible  after  severe  TBI  to 
deaccelerate  secondary  brain  injury.  TCD  is  a  simple  and  noninvasive  method  for 
cerebral  perfusion  assessment  and  is  particularly  accurate  to  detect  hypoperfusion. 
Catherine  et  al.  confirmed  that  early  TCD  makes  it  possible  to  identify  high-risk 
patients  with  impaired  cerebral  perfusion  and  poor  outcome.  Moreover,  this  study 
(Catherine  et  al.  2007)  has  also  discovered  that  TCD  goal-directed  therapy  has  the 
potential  to  improve  cerebral  perfusion  before  invasive  cerebral  monitoring  is  avail¬ 
able,  thereby  limiting  secondary  brain  injuries.  Although  study  suggested  the  need 
for  further  studies  to  define  optimal  treatment  TCD  threshold  values  for  the  initial 
management  of  TBI,  and  its  effect  on  outcome  (Catherine  et  al.  2007). 

4.4.6  Magnetic  Resonance  Imaging 

MRI  is  perhaps  the  most  used  imaging  modality  for  studying  the  anatomy  and  func¬ 
tion  of  the  brain  today.  The  findings  in  standard  MRI  studies  of  TBI  are  similar  to  CT, 
but  MRI  is  more  sensitive  than  CT  for  small  lesion  detection  such  as  TAI  and  DAI  as 
seen  with  SWI  (see  above).  MRI  also  offers  a  variety  of  different  measures  to  study 
TBI  such  as  Diffusion  Tensor  Imaging  (DTI),  Magnetic  Resonance  Spectroscopy 
Imaging  (MRSI),  and  functional  magnetic  resonance  imaging  (fMRI)  (Haacke  et  al. 
2011)  which  are  not  part  of  the  focus  of  this  chapter.  There  are  three  major  methods 
that  are  used  at  the  moment  to  study  either  local  or  global  perfusion  to  the  brain. 

4.4.6.1  DSC  Perfusion  Weighted  Imaging  (PWI) 

PWI  is  sensitive  to  microscopic  blood  flow.  The  method  is  based  on  monitoring  a 
non-diffusible  contrast  material  gadolinium  DTPA  (dimeglumine  pentaacetic  acid 
or  diethylenetriamine  pentaacetic  acid)  passing  through  brain  tissue.  The  signal 
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363  intensity  decreases  as  the  contrast  bolus  passes  through  the  tissue  and  returns  to 

364  normal  over  time.  Recording  this  time-varying  signal,  and  extracting  the  concen- 

365  tration  as  a  function  of  time,  it  is  possible  to  derive  the  following  measures  of 

366  blood  flow,  cerebral  blood  volume  (CBV);  CBF;  mean  transit  time  (MTT)  and  time 

367  to  peak  (TTP),  as  well  as  a  variety  of  other  hemodynamic  parameters.  To  date, 

368  there  are  few  studies  describing  the  use  of  MR  perfusion  in  TB1  evaluation  or  pre- 

369  dieting  outcomes  after  TBI.  One  small  study  (Garnett  et  al.  2001)  showed  in 

370  patients  with  focal  injuries  visible  in  T2  images  that  reductions  of  rCBV  were 

371  found  in  focal  lesions  and  significant  increases  of  rCBV  were  found  in  pericontu- 

372  sional  areas.  This  finding  is  also  consistent  with  the  concept  of  traumatic  penumbra 

373  (Menon  2003).  Furthermore,  a  subset  of  these  patients  also  showed  reduced  rCBV 

374  in  normal  appearing  white  matter,  and  these  patients  had  significantly  poorer  out- 

375  comes  (Garnett  et  al.  2001).  Secondary  damage  in  TBI  is  usually  a  result  of  altered 

376  hemodynamics  and  cerebrovascular  factors.  Hayward  et  al.  investigated  hemody- 

377  namic  coupling  between  CBF  and  CBV.  In  this  work,  surprising  results  were  found 

378  that  CBV  was  the  same  for  patients  and  controls  in  all  regions  considered. 

379  A  significant  decrease  in  CBF  was  found  at  48  h  after  TBI  but  the  CBV  remained 

380  unchanged  (Hayward  Nick  Mark  2011).  Our  own  work  in  this  direction  has  just 

381  begun.  However,  we  have  one  example  with  deep  MVD  and  PWI  that  appears  to 

382  show  several  regions  of  reduced  mean  transit  time  (Fig.  4.5).  If  this  pattern  is  found 

383  to  be  replicated  in  other  cases,  it  would  demonstrate  that  the  damaged  area  is  also 

384  associated  with  reduced  perfusion  (Fig.  4.5). 

385  4.4.6.2  Arterial  Spin  Labeling 

386  ASL  is  an  arterial  spin  tagging  technique  that  relies  on  the  detection  of  magnetically 

387  labeled  blood.  Once  the  magnetization  of  the  inflowing  arterial  blood  has  been 

388  modified  (generally  inverted)  upstream,  it  induces  a  small  MR  signal  change  down- 

389  stream  (a  few  percent  of  the  tissue  magnetization).  Meanwhile,  the  magnetization  of 

390  the  perfusion  tracer  (i.e.,  the  labeled  water)  is  rapidly  relaxing.  Return  of  the  longi- 

391  tudinal  magnetization  toward  its  equilibrium  values  takes  a  few  seconds.  Although 

392  ASL  techniques  have  not  entered  widespread  clinical  usage,  their  utility  has  been 

393  demonstrated  for  a  variety  of  acute  and  chronic  cerebrovascular  diseases  such  as 

394  stroke  and  epilepsy.  (Wintermark  et  al.  2005).  One  significant  advantage  of  ASL  is 

395  the  ability  to  perform  multiple  repeated  measurements  without  the  use  of  a  contrast 

396  agent.  This  might  be  necessary  before  and  after  treating  the  patient  with  a  cerebro- 

397  vascular  dilator  (acetazolamide)  or  vasoconstrictor  inhibitor  (clazosentan),  for 

398  example.  Kim  et  al.  used  ASL  to  measure  the  resting  state  CBF  in  28  moderate  to 

399  severe  TBI  patients  in  the  chronic  stage  and  found  that  (a)  there  was  a  global 

400  decrease  of  CBF  compared  to  controls  and  (b)  prominent  regional  hypoperfusion 

401  was  present  in  the  posterior  cingulate  cortices,  the  thalami,  and  multiple  locations  in 

402  the  frontal  cortices  (Kim  et  al.  2010).  However,  in  the  population  of  mTBI,  how 

403  cerebral  perfusion  is  being  affected  in  patients  and  their  association  with  patients’ 

404  neurocognitive  status  is  still  not  understood.  Given  the  milder  severity  in  mTBI 
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Fig.  4.5  An  example  of  TBI  case  with  frontal  medullary  vein  damage  and  microbleeds,  (a)  FLAIR 
image  showing  two  of  the  damaged  frontal  areas  (short  arrow )  and  enhanced  signal  in  the  sple- 
nium  of  the  corpus  callosum  ( long  arrow),  (b)  SWI  data  showing  multiple  bleeds  in  the  frontal, 
lateral,  and  the  splenium  of  the  corpus  callosum  (short  arrows),  (c)  T2  weighted  scan  showing 
similar  results  to  the  FLAIR  data,  (d)  Mean  transit  time  images  from  the  PWI  data  showing  pos¬ 
sible  delayed  arrival  times  in  frontal  regions  and  medial  gray  matter  regions  (black  arrows)  that 
may  represent  reduced  perfusion  to  these  tissues.  Note  the  dark  regions  correspond  to  the  bleeds 
and  because  of  the  long  echo  times  used  for  PWI,  the  signal  there  was  too  small  to  be  used  to 
determine  PWI  parameters  correctly  and  hence  are  set  to  zero  here 


patients,  perfusion  imaging  at  the  acute  stage  could  be  meaningful  both  for  clinical  405 
decision-making  and  patients’  outcome  prediction.  Ge  et  al.  also  employed  ASL  to  406 
study  21  mTBI  patients  at  the  chronic  stage  and  demonstrated  reduced  CBF  in  both  407 
sides  of  the  thalamus,  which  is  correlated  with  the  patients’  speed  of  information  408 
processing,  memory,  verbal,  and  executive  function  (Ge  et  al.  2009).  ASL  has  also  409 
been  used  in  conjunction  with  functional  brain  imaging  studies  and  with  animal  410 
studies,  but  these  directions  are  outside  the  scope  of  the  material  presented  in  this  41 1 
chapter.  However,  some  of  the  fMRI/ASL  research  does  suggest  that  there  are  412 
reductions  in  CBF  that  may  be  related  to  diffuse  pathology.  413 
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414  4.4.6.3  Blood  Flow  Measurements  in  the  Neck  and  the  Poor  Man’s  Bulk 

415  Perfusion  Measure 

416  Phase-contrast  MR  images  are  used  for  blood  flow  measurement  in  the  neck  and 

41 7  brain.  Flow  can  be  obtained  for  all  major  vessels  in  the  neck  thereby  allowing  a  mea- 

41 8  sure  of  the  cardiovascular  input  and  output  to  the  brain.  By  measuring  flow  at  C 1/C2, 

419  one  can  focus  on  just  the  internal  carotid  artery  and  vertebral  artery.  By  normalizing 

420  this  arterial  inflow  to  brain  volume,  it  is  possible  to  obtain  the  average  flow  in  units 

421  of  mL/100  g  tissue/min.  CSF  flow  can  also  be  measured  at  this  level.  Similarly,  by 

422  measuring  flow  at  C7/T1,  the  total  venous  outflow  can  be  determined.  Flow  is  actu- 

423  ally  measured  throughout  the  cardiac  cycle  (on  the  order  of  2CM-0  points  per  cardiac 

424  cycle).  A  number  of  flow  measures  can  be  then  determined  including,  flow  velocity 

425  (average,  peak  positive,  and  peak  negative),  volume  flow  rate  (average,  positive,  and 

426  negative),  and  total  flow  volume  per  cardiac  cycle.  The  following  arteries  and  veins 

427  of  left  and  right  sides  are  measured:  common  carotid  arteries,  internal/external  carotid 

428  arteries,  vertebral  arteries,  internal/external  jugular  veins,  vertebral  veins,  deep  cervi- 

429  cal  veins,  anterior  jugular  vein,  and  any  other  major  veins  that  are  apparent. 


430  4.5  Measuring  Oxygen  Saturation  Changes  with  SWIM 

431  Blood  oxygen  saturation  imaging  provides  a  means  to  monitor  changes  in  oxygen 

432  saturation  indirectly.  SWI  can  be  used  to  not  only  visualize  small  veins  and  microb- 

433  leeds  but  also  provide  phase  information  for  major  veins.  Such  changes  can  occur 

434  when  there  is  a  change  in  blood  flow.  When  the  local  field  changes  because  of 

435  deoxyhemoglobin  changes,  the  phase  changes  in  and  around  the  major  veins.  In 

436  Fig.  4.6  below,  we  show  the  effect  of  a  systemic  reduction  in  blood  flow  to  the  brain 

437  after  ingestion  of  200  mg  of  caffeine  (the  subject  ingested  a  200  mg  NoDoz  pill). 

438  Measuring  changes  in  phase  is  tantamount  to  measuring  changes  in  oxygen  satura- 

439  tion  (Shen  et  al.  2007).  Similarly,  reductions  in  perfusion  or  local  bleeding  can  lead 

440  to  the  apparent  enlargement  of  a  local  medullary  vein  (see  Figs.  4.1  and  4.2)  caused 

441  by  either  this  blood  oxygenation  level  dependent  (BOLD)  effect  or  by  local  bleed- 

442  ing  in  and  around  the  vessel. 

443  Applying  susceptibility  mapping  to  the  phase  yields  a  SWIM  image  (Fig.  4.7b) 

444  that  looks  similar  to  the  SWI  data  of  the  venous  system  (Fig.  4.7a).  Although  the 

445  latter  is  important  in  detecting  vascular  abnormalities,  it  is  not  quantitative.  With 

446  high-enough  resolution  and  improved  SWIM  algorithms,  it  is  expected  that  major 

447  veins  (roughly  2  mm  or  greater  in  diameter)  will  yield  reasonably  accurate  values 

448  for  the  oxygen  saturation  in  those  veins  (Figs.  4.6  and  4.7). 
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Fig.  4.6  SWI  projection  images  showing  the  same  individual  (a)  pre-caffeine  and  (b)  post-caffeine. 
The  patient  had  ingested  one  No-Doze  pill  containing  200  mg  of  caffeine.  Note  the  increased 
contrast  of  the  smaller  veins.  This  is  caused  by  the  increased  levels  of  deoxyhemoglobin  after 
caffeine.  Since  caffeine  is  a  vasoconstrictive  agent,  the  blood  flow  slows  to  the  brain.  However,  since 
the  cerebral  metabolic  rate  of  oxygen  utilization  (CMR02)  does  not  change,  the  amount  of  deoxy¬ 
hemoglobin  increases 


Fig.  4.7  SWI  versus  SWIM,  (a)  SWI  projection  image  from  Fig.  4.6b  showing  the  major  veins  in 
the  brain,  (b)  SWIM  image  (of  a  different  individual)  showing  a  quantitative  map  of  oxygen  satura¬ 
tion  in  the  brain 
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449  4.6  Role  of  Finite  element  Modeling  in  Understanding  TBI 

450  Finite  element  (FE)  models  of  brain  injury  are  essential  to  study  brain  injury 

451  mechanisms.  A  number  of  studies  have  been  conducted  to  investigate  the  contri- 

452  bution  of  the  vasculature  on  the  dynamic  response  of  the  brain  in  traumatic  injury 

453  (Zhang  et  al.  2002).  The  hypothesis  is  that  the  vascular  tree  contributes  to  the 

454  strength  of  the  neural  tissues  of  the  brain  when  it  is  subjected  to  applied  loads. 

455  Data  on  mechanical  properties  of  vascular  tissues  (Monson  et  al.  2000;  Arbogast 

456  and  Margulies  1997)  suggest  that  cerebral  vessels  are  much  stiffer  than  brain  tissue. 

457  Two-dimensional  FE  models  of  the  brain  with  and  without  the  simulation  of  major 

458  cerebral  arteries  were  evaluated  to  test  the  hypothesis  that  the  structural  stiffness 

459  of  the  human  brain  is  influenced  by  the  pressure  of  these  vessels.  Overall  results 

460  from  this  study  suggests  that  explicit  modeling  of  blood  vessels  increases  the  stiff- 

461  ness  of  the  brain  to  resist  both  rotational  and  linear  impulses  and  increases  the 

462  threshold  of  neural  injury.  Zhang  et  al.  assumed  in  her  study  that  the  blood  vessel 

463  walls  were  directly  connected  to  the  surrounding  brain  tissues  (Zhang  et  al.  2002). 

464  However,  there  is  no  data  that  describe  the  mechanical  interaction  between  the 

465  vessels  and  brain. 

466  Medullary  vein  damage  may  be  due  to  the  presence  of  internal  shear  forces 

467  (intracerebral  stresses).  How  this  occurs  likely  depends  on  the  stiffness  and  natu- 

468  ral  frequency  of  the  structure  impacted.  For  impacts  with  structures  of  relatively 

469  low  stiffness,  the  brain  tends  to  move  with  the  skull  because  of  the  predomi- 

470  nantly  low  frequencies  involved  in  these  impacts.  These  movements  result  in 

471  intracerebral  strains  (ICS)  that  can  produce  lesions  in  the  central  parts  of  the 

472  brain.  These  ICS  usually  occur  when  the  natural  frequency  of  the  impact  is  less 

473  than  about  250  Hz,  with  the  result  that  the  outer  parts  of  the  brain  move  with  the 

474  skull  causing  shear  strains  in  the  cerebral  tissue  as  the  inner  parts  of  the  brain  lag 

475  behind  the  outer  due  to  initial  effect.  A  theoretical  model  suggests  that  the  natu- 

476  ral  frequencies  of  the  head  and  the  structure  are  the  characteristics  with  most 

477  importance  in  determining  the  distribution  of  injury  in  the  brain  for  a  particular 

478  location  of  impact  on  the  head  (Willinger  et  al.  1994).  Clearly,  knowing  the  tissue 

479  properties  and  the  local  structures  of  the  vessels  and  how  they  are  embedded  in  the 

480  surrounding  tissue  is  necessary  for  proper  modeling  of  the  brain  injury  due  to  a 

481  given  loading  condition.  The  SWI  and  MRA  data  we  have  collected  makes  it 

482  possible  to  create  an  anatomically  detailed  FE  model  of  the  medullary  veins  down 

483  to  250  pm  in  size  (see  Fig.  4.8).  Armed  with  this  detailed  local  model  embedded  in 

484  the  whole  brain  model  and  estimates  for  the  tissue  properties,  it  may  be  possible  in  the 

485  future  to  predict  local  venous  damage  and  surrounding  neural  tissue  injury.  Some 

486  early  attempts  to  model  the  vasculature  from  MRA  data  were  presented  by  Parnaik 

487  et  al.  (2004)  (Fig.  4.8). 
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Fig.  4.8  Modeling  the  medullary  veins.  An  SWI  projection  image  is  used  to  find  medullary  veins 
and  then  they  are  extracted  in  3D.  (a)  An  example  image  from  SWI  data  showing  the  septal  vein 
and  the  confluence  of  several  medullary  veins,  (b)  A  projection  of  3D  model  of  medullary  veins 
branching  off  the  septal  vein.  The  resolution  in  the  model  is  62.5  pm  and  vessel  range  in  size  from 
1  mm  (the  septal  vein)  to  the  first  branchings  (0.5  mm)  to  the  second  branchings  (0.25  mm).  These 
veins  can  be  damaged  in  frontal  or  side  impact 

4.7  Conclusions  and  Future  Directions 

SWI  has  proven  to  be  one  of  the  methods  of  choice  along  with  DTI  and  T2  weighted 
FLAIR  (Fluid  Attenuation  Inversion  Recovery)  when  imaging  TBI.  SWI  detects 
abnormalities  caused  by  focal  lesions  and  DAI,  many  of  which  may  be  missed  in 
conventional  MRI  and  CT  scans.  Currently,  there  is  a  paucity  of  information  about 
perfusion  in  human  TBI  studies.  It  is  likely  that  perfusion  could  play  a  major  role  in 
understanding  patient  outcomes.  Perfusion  combined  with  oxygen  saturation  mea¬ 
surements  using  SWIM  may  also  open  the  door  for  a  better  understanding  of  the 
type  of  damage  done  in  TBI.  However,  there  is  still  one  missing  piece  to  the  puzzle 
that  is  a  critical  precursor  to  understanding  flow  changes  in  TBI  and  that  is  the  arte¬ 
rial  system.  Despite  our  focus  on  measuring  blood  products,  veins,  oxygen  satura¬ 
tion,  and  perfusion,  it  is  critical  to  gain  an  understanding  of  damage  to  the  arterial 
system  as  well.  Practically,  the  ideal  protocol  would  be  to  start  with  an  MR  angio¬ 
gram  (MRA)  and  venogram  (SWI)  at  the  beginning  of  the  scanning  for  TBI  patients. 
Recently,  we  have  implemented  a  double  echo  sequence  that  can  be  used  to  create 
both  MRA  and  SWI  data  for  such  an  analysis  (Haacke  et  al.  2012).  This  should  then 
be  followed  by  a  perfusion  scan  and  any  other  key  sequences  such  as  FLAIR  and 
DTI.  We  anticipate  that  it  may  then  be  possible  to  better  diagnose  many  more  mild 
TBI  patients  than  is  otherwise  possible  today. 
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2.1  Introduction 

Traumatic  brain  injury  (TBI)  affects  1.5  million  Americans  each  year,1,2  the 
majority  of  which  are  mild  TBI  (mTBI).3  Today,  more  people  who  might  have 
had  more  severe  TBI  or  died  in  motor  vehicle  accidents  suffer  milder  head 
injury,  thanks  to  the  advancements  in  motor  vehicle  safety.  mTBI  affects  over  1 
million  Americans  annually.  It  is  a  major  public  healthcare  burden,  yet  has  been 
overlooked  for  decades.1'2  Despite  its  name,  the  impact  of  mTBI  on  the  patients 
and  their  family  are  not  mild  at  all.3  mTBI  patients  typically  develop  a  variable 
constellation  of  physical,  cognitive,  and  emotional  symptoms,  collectively 
known  as  post-concussive  syndrome  (PCS),  that  significantly  impact  their 
quality  of  life.  The  direct  cost  of  mTBI  in  the  U.S.  is  approximately  $16.7  billion 
annually  and  this  does  not  include  the  indirect  costs  to  society  and  families.3  5 

Up  to  50%  of  patients  with  mTBI  have  persistent  neurocognitive  symptoms  at 
three  months  and  5%  to  15%  at  1  year.6  7  Meanwhile,  among  the  1.8  million 
troops  who  have  served  in  Iraq  and  Afghanistan,  it  is  estimated  that  at  least  20% 
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of  returning  troops  have  suffered  at  least  one,  and  some  multiple,  concussions. 
This  means  up  to  360000  veterans  may  have  brain  injuries  after  discharge.8 
However,  most  symptomatic  rnTBI  patients  have  normal  CT  and  conventional 
MRI  scans,910  even  with  the  greater  sensitivity  of  MRI  to  TB1. 1 1  In  addition  to 
neuroimaging,  clinical  indices  of  severity,  such  as  the  Glasgow  Coma  Scale 
(GCS)  and  duration  of  post-traumatic  amnesia  (PTA)  are  lacking  in  sensitivity  in 
mTBI  and  are  not  helpful  in  predicting  outcome.12  In  summary,  using  currently 
available  clinical  instruments,  it  is  difficult  to  determine  which  mTBI  patients  will 
have  prolonged  or  even  permanent  neurocognitive  symptoms. 

With  greater  appreciation  of  incidence  and  impact  of  “civilian,”  sports,  and 
military  mTBI,  there  is  an  increasing  need  to  identify  a  minimally  invasive  bio¬ 
marker  (or  biomarkers)  of  mTBI.  A  clinically  useful  biomarker  or  biomarkers 
should  ideally  be  very  sensitive  to  TBI  (including  mild  cases)  but  should  have 
reasonable  specificity.  That  TBI  pathology  is  not  singular  but  involves  multiple 
functional  and  structural  changes  affords  greater  opportunity  for  detection  and 
characterization  of  injury  with  sensitive  instruments.  Similarly,  biomarkers 
meeting  these  criteria  should  improve  clinical  trial  efficiency  since  patients  with 
similar  types  and  grades  of  pathology  could  be  grouped  together  in  order  to  reduce 
the  chances  of  a  negative  result  due  to  prohibitive  pathologic  heterogeneity. 

With  both  hardware  and  pulse  sequence  design  advances,  newer  MRI 
methods  have  demonstrated  the  ability  to  detect  and  localize  with  high  reso¬ 
lution  several  of  the  pathologic  and  pathophysiologic  consequences  of  mTBI. 
These  advanced  MR  technologies  include  susceptibility  weighted  imaging 
(SWI)  for  hemorrhage  detection,13  MR  spectroscopy  (MRS)  for  metabolite 
measurement,14  diffusion  weighted  and  diffusion  tensor  imaging  (DWI/DTI) 
for  edema  quantification1^  and  axonal  injury  detection,16  perfusion  weighted 
imaging  (PWI)  and  arterial  spin  labeling  (ASL)  to  measure  blood  flow  to  brain 
tissue,17  and  functional  MRI  (fMRI),  which  measures  changes  in  blood  oxygen 
level  locally  in  response  to  neuronal  activity.18  Having  a  number  of  imaging 
biomarkers,  all  of  which  are  obtained  in  a  single  scanning  session  (or  multiple 
for  longitudinal  study)  and  are  sensitive  to  different  consequence  of  traumatic 
injury  affords  great  advantages:  1)  enhanced  sensitivity;  2)  ability  to  study  inter¬ 
relationships  among  these  biomarkers  and  between  the  biomarkers  and  clinical/ 
neurocognitive  deficits;  3)  improved  clinical  management  resulting  from  more 
precise  characterization  of  injuries;  and  4)  enhanced  power  of  clinical  inter¬ 
ventional  studies  (see  above).  Furthermore,  the  information  comes  “cheaply” 
since  MRI  is  non-invasive  and  is  widely  available. 

2.2  Search  Criteria  and  Results 

The  purpose  of  this  chapter  is  to  review  the  current  status  and  clinical  impact  of 
MR  imaging  findings  in  mTBI.  BeLanger  et  al.19  provided  an  informative  review 
on  this  topic  in  2007.  Since  then  many  original  works  have  been  published. 

The  following  inclusion  criteria  were  used  to  identify  potential  papers:  a) 
must  be  cited  in  Medline  database;  b)  must  use  MRI  as  a  major  component  of 
the  study;  c)  only  clinical  (human)  studies  of  mTBI  patients;  d)  published  in  the 
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past  five  years  (i.e.  2006  and  later);  and  e)  published  in  peer-reviewed  English 
language  journals.  Readers  interested  in  moderate  to  severe  TBI  imaging  are 
encouraged  to  refer  to  another  recent  review  paper,  which  covers  SWI,  DTI, 
and  MRS.19 

The  authors  conducted  a  Medline  search  at  www.pubmed.com  in  September 
2011  by  using  the  search  phrases  “mild  traumatic  brain  injury”  and  “magnetic 
resonance  imaging.”  There  were  a  total  356  publications,  of  which  46  were 
reviews  and  310  were  original  studies.  After  excluding  articles  that  did  not  fit 
our  criteria,  we  found  a  total  62  articles;  11  reviews  and  51  original  studies.  Of 
original  studies,  32  used  structural  imaging,  DTI,  MRS,  and  perfusion  imaging 
and  19  articles  used  fMRI,  including  both  resting  state  fMRI  and  task-oriented 
fMRI.  This  review  focuses  on  clinical  MRI,  DTI,  MRS,  and  perfusion  imaging. 
The  following  was  categorized  according  to  the  technical  approaches  being 
used.  The  literature  is  listed  in  Table  2.1. 


2.3  Clinical  CT  and  MRI 

Although  CT  is  still  the  first  step  in  the  search  for  large  hemorrhages  that  may 
require  surgical  intervention,  it  is  insensitive  to  small  hemorrhages,  whether 
from  early  contusion  or  from  diffuse  axonal  injury.  Most  patients  with  mTBI 
have  normal  CT  and  conventional  MRI  findings.  Indeed,  they  usually  do  not 
get  clinical  MRI  scans  in  the  emergency  department,  which  sees  the  majority  of 
patients  with  mTBI.  Any  intracranial  abnormality  would  classify  the  patients 
as  complicated  mTBI,20  which  could  have  long-term  sequelae  similar  to  mod¬ 
erate  TBI.  A  clinical  radiological  sequence  protocol  typically  includes  Tl- 
weighted  imaging  for  structural  damage  and  brain  atrophy,  T2-weighted 
FLAIR  for  edema  and  acute  stage  subarachnoid  hemorrhage,  and  T2*  gradient 
recalled  echo  (GRE)  for  hemorrhages. 

It  would  be  of  clinical  value  to  know:  “How  many  patients  with  mTBI  have 
positive  MRI  findings?”  and  “How  these  findings  may  account  for  their  neu- 
rocognitive  or  functional  symptoms?”  Regarding  the  prevalence  of  positive 
MRI  findings  in  mTBI,  different  groups  reported  different  findings.  After 
reviewing  182  adolescent  patients,  Fay  et  al?x  reported  18%  had  positive  MRI 
findings.  Topal  et  air2  reported  that  of  40  patients  with  mTBI,  12.5%  had 
positive  MRI.  In  contrast,  Lee  and  colleagues23  reported  50%  with  positive  CT 
and  75%  with  positive  MRI  in  36  mTBI  patients  with  both  loss  of  consciousness 
and  post-traumatic  amnesia.  The  discrepancy  among  the  figures  could  be  due  to 
the  different  study  populations  and  different  patient  inclusion  criteria. 

Regarding  the  relationship  between  MRI  findings  and  patient  outcome, 
Muller  et  al.24  reported  that  GCS<14,  CT/MRI  positive,  and  S100B>0.14, 
together,  predicted  impaired  cognitive  performance  both  at  baseline  and  after  6 
months.  Kurca  et  al r5  reported  that  patients  with  mTBI  with  true  traumatic 
MRI  lesions  are  different  neuropsychologically  from  those  of  mTBI  with  non¬ 
specific  MR  lesions,  and  neurocognitive  and  symptom  signs  have  organic  bases 
that  can  be  detected  by  MRI.  Fay  et  al.2>  reported  that  ratings  of  PCS  in 
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Table  2.1  Summary  of  literature  search  result. 

First  author  Imaging  Number  of  Inclusion  Analytical 

and  year  approach  patients  criteria  Imaging  stage  approach  Major  findings 

Clinical  MRI 

Fay  201021  Clinical  MRI  182  children  LOC  or  Acute  Clinical  18%  with  positive  MRI  finding 
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Inglese  2006102  MRI  38  mTBI  Acute  and  chronic  Clinical  radi-  Average  number  of  Virchow  Robin 

ological  Space  was  significantly  higher  in 

diagnosis  mTBI  patients  than  in  controls. 

VRS  was  not  associated  with  neu- 
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Table  2.1  ( Continued ) 

First  author  Imaging  Number  of  Inclusion  Analytical 

and  year  approach  patients  criteria  Imaging  stage  approach  Major  findings 

Diffusion  imaging 

Grossman  Diffusion  22  mTBI  ACRM  Subacute  to  Diffusion  Kurtosis  Changed  DTI  and  DIG  measures  in 
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Table  2.1  ( Continued ) 

First  author  Imaging  Number  of  Inclusion  Analytical 

and  year  approach  patients  criteria  Imaging  stage  approach  Major  findings 

Niogi  20  0850  DTI,  T2*  34  mTBI  with  GCS  13-15,  Subacute  to  ROI  on  major  FA  decreased  in  several  WM  regions: 
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adolescents  were  moderated  by  both  cognitive  ability  and  injury  severity. 
Children  with  lower  cognitive  ability  with  complicated  mTBI  were  prone  to 
PCS.  Comparatively,  Lee  et  a/.23  reported  that  CT  and  MRI  findings  were  not 
associated  with  a  patient’s  neurocognitive  impairment.  Clearly,  the  seemingly 
conflicting  findings  on  the  predictive  role  of  structural  MRI  need  further 
investigation. 


2.3.1  Susceptibility  Weighted  Imaging  (SWI)  of  Hemorrhagic 
Lesions 

In  diagnostic  radiology,  intracranial  hemorrhage  has  been  sought  as  a  bio¬ 
marker  of  TB1.  Some  investigators  have  suggested  that  the  presence  of 
hemorrhage  in  DAI  is  predictive  of  poor  outcome  in  moderate  to  severe  TBI. 
SWI  was  developed  by  Haacke  et  air1  as  a  high  resolution  venography  method. 
It  has  been  used  to  evaluate  moderate  to  severe  TBI  patients  since  2003  by  Tong 
et  a/.,13  who  have  shown  that  SWI  is  3  to  6  times  more  sensitive  than  con¬ 
ventional  T2*  gradient  echo  imaging  (GRE)  for  detecting  suspected  DAI 
lesions  in  children.13’28  SWI  has  been  shown  to  detect  tiny  hemorrhages  that 
may  be  the  only  abnormal  finding  that  can  confirm  the  presence  of  brain  injury, 
and  change  management  of  the  patient.  In  addition,  lesion  number  and  volume 
identified  by  SWI  are  negatively  associated  with  patient  outcome28  and  neu¬ 
ropsychological  functions29  in  patients  with  moderate  to  severe  injury. 

After  brain  injury,  the  hemorrhagic  bleed  may  undergo  a  temporal  trans¬ 
formation  from  oxyhemoglobin  and  then  deoxyhemoglobin  in  the  acute  stage; 
intra-  and  then  extracellular  methemoglobin  in  the  subacute  stage;  and  finally, 
hemosiderin  in  the  chronic  stage.30  32  SWI  is  sensitive  to  deoxyhemoglobin  in 
the  acute  stage  and  extracellular  methemoglobin  in  the  subacute  stage  and 
hemosiderin  in  the  chronic  stage.30  Therefore,  positive  results  with  SWI  should 
provide  a  biomarker  for  hemorrhagic  brain  injury  at  any  stage.  However,  very 
few  data  have  been  reported  regarding  the  possible  role  of  SWI  in  an  improved 
detection  of  microhemorrhages  and  its  predictive  value  of  patient  outcome  in 
mTBI.  Unlike  DTI,  which  requires  complex  post-processing  and  comparison 
with  proper  controls,  SWI  hemorrhage  is  readily  available  for  radiological 
reading  immediately  after  MRI  scans.  This  makes  SWI  more  likely  to  have  a 
direct  impact  on  the  radiological  diagnosis.  Furthermore,  with  quantitative 
susceptibility  mapping,  the  hemorrhagic  lesion  load  could  be  more  acutely 
mapped  through  next  generation  SWI  analysis  and  potentially  improve  the 
prediction  of  a  patient’s  outcome.  Figure  2.1  demonstrates  an  exemplary  case 
to  show  the  SWI  detection  of  brain  hemorrhages  in  mTBI  in  the  acute  stage  in 
comparison  with  the  usual  clinical  T2*  GRE. 

In  summary,  there  is  still  a  need  to  evaluate  a  large  number  of  patients  with 
mTBI  by  using  structural  imaging,  particularly  SWI,  to  see  the  clinical  value  of 
intracranial  hemorrhage  detected  by  MRI.  This  should  be  done  along  with 
other  clinical  measures,  e.g.  GCS,  as  well  as  other  biological  markers,  such  as 
serum  biomarkers. 
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Figure  2.1  Comparison  of  T2*  with  SWI  on  the  detection  of  hemorrhagic  lesions. 

Image  order:  T2*  GRE  (left)  and  SWI  (right).  A  45-year-old  man  fell 
down  a  staircase  and  then  visited  ER  with  severe  headache.  Both  T2* 
GRE  and  SWI  detected  hemorrhages  on  the  right  side  of  inferior  temporal 
lobe.  However,  in  a  small  area  with  mixed  blood  and  edema  signal,  T2* 
GRE  only  detects  edema  (bright  signal)  and  SWI  detects  hemorrhage 
(dark  signal)  (see  arrows). 

2.4  Diffusion  Tensor  Imaging  (DTI)  of  Axonal  Injury 

Diffusion  imaging  sequences  are  sensitive  to  traumatic  axonal  injury  (TAI) 
secondary  to  stretch  and  shear  forces.  DTI  measures  the  bulk  motion  of  water 
molecular  diffusion  in  biological  tissues.  It  is  most  useful  when  tissues  are  ani¬ 
sotropic,  i.e.  when  diffusion  is  not  equivalent  in  all  directions,  such  as  in  skeletal 
muscle  or  axons  in  white  matter  of  the  central  nervous  system.  Histological 
correlates  have  validated  DTI’s  sensitivity  to  brain  injury  for  both  focal33  and 
diffuse  axonal  injury34  models.  The  apparent  diffusion  coefficient  (ADC)  and 
fractional  anisotropy  (FA)35'36  are  two  parameters  derived  from  DTI  that  have 
been  extensively  studied  in  TBI.  ADC  is  an  estimate  of  the  average  magnitude  of 
water  movement  in  a  voxel  (regardless  of  direction),  while  FA  is  an  index  of  the 
directional  non-uniformity,  or  anisotropy,  of  water  diffusion  within  a  voxel. 

FA  has  been  used  to  detect  alterations  in  directional  diffusion  resulting  from 
tissue  damage.  FA  in  white  matter  is  highest  when  fibers  are  long  (relative  to 
voxel  dimension)  and  oriented  uniformly  (collinear)  within  a  voxel  and  lowest 
when  fibers  are  not  collinear  (e.g.  “crossing  fibers”)  or  have  been  damaged. 
When  axons  are  injured,  as  in  acceleration/deceleration  injuries  (such  as  MVCs), 
diffusion  anisotropy  typically  decreases.  Loss  of  diffusion  anisotropy  is  the  result 
of  a  number  of  axonal  changes  after  injury  including:  1)  increased  permeability 
of  the  axonal  membrane;  2)  swelling  of  axons;  3)  decreased  diffusion  in  the  axial 
(long  axis)  direction;  4)  degeneration  and  loss  of  axons  in  the  chronic  stage.  In 
general,  any  pathological  alteration  of  white  matter  fibers  will  result  in  FA 
decrease,  since  one  or  more  of  these  axonal  changes  occurs  in  disorders  of  white 
matter.  Not  surprisingly,  most  clinical  studies  of  moderate  and  severe  TBI  have 
shown  FA  to  be  more  sensitive  than  ADC  to  traumatic  injury.  On  the  other 


1 

5 

10 

15 

20 

25 

30 

35 

40 

45 


Magnetic  Resonance  Imaging  Biomarkers  of  Mild  Traumatic  Brain  Injury 


29 


Figure  2.2  SWI  detects  subarachnoid  hemorrhage  (SAH)  in  an  acute  rnTBI  case. 

Image  order  from  left  to  right:  CT,  SWI  phase,  and  T2  FLAIR.  A  21-year- 
old  man  fell  from  a  15-foot  (4.5  m)  high  ladder.  Both  CT  and  SWI  phase 
detected  sub-arachnoid  hemorrhage  (see  arrows)  at  the  acute  stage,  while 
F2  FLAIR  failed  to  detect. 

hand,  ADC,  FA,  and  directionally  selective  difTusi vitics  (principal,  intermediate, 
and  minor  components  of  diffusion)  can  help  to  better  characterize  brain  injury 
pathologies.  Trace  and  mean  diffusivity  are  two  other  measurements  similar  to 
ADC  and  vary  similarly.  Changes  in  FA  in  association  with  ADC  changes  can 
differentiate  the  type  of  edema.  For  example,  in  the  acute  stage,  decreased  FA  in 
association  with  increased  ADC  suggests  vasogenic  edema,  while  increased  FA 
in  association  with  decreased  ADC  suggests  cytotoxic  edema.  Decreased  FA  in 
association  with  decreased  ADC  and  decreased  longitudinal  (parallel  to  the  long 
axis  of  the  axons)  water  diffusivity  suggests  axonal  transport  failure  as  occurs  in 
degenerative  neurological  diseases  such  as  in  ALS. 

Regarding  the  location  of  brain  lesions  detected  by  DTI,  Niogi  et  al ,37 
summarized  that  the  frontal  association  pathways,  including  anterior  corona 
radiata,  uncinate  fasciculus,  superior  longitudinal  fasciculus,  and  genu  of 
corpus  callosum  (CC),  were  the  mostly  frequently  injured  WM  structures  in  this 
cohort  of  mTBI  patients.  Single  subject  results  were  not  reported  but  it  can  be 
assumed  that  significant  inter-individual  variability  for  location  and  extent  of 
WM  injury  exists  due  to  varying  injury  mechanisms  (and  forces)  and  biological 
(neural  and  non-neural)  differences  across  patients. 

2.4.1  Imaging  at  Different  Pathological  Stages 

Interestingly,  despite  the  higher  incidence  of  milder  TBI  compared  with  more 
severe  TBI  in  Western  countries,  there  are  many  fewer  mTBI  imaging  studies 
reported  in  the  literature.  One  reason  is  the  fact  that  recruitment  of  patients  with 
mTBI  is  more  difficult,  since  they  are  typically  outpatients.  Another  reason  is  the 
common  conception  that  mTBI  is  a  transient  problem  from  which  virtually  all 
who  are  afflicted  will  recover  fully.  Therefore,  some  question  the  clinical 
importance  of  studying  mTBI.  Certainly,  insurance  reimbursement  for  MRI 
scanning  of  mTBI  is  rare  and  even  rarer  in  the  acute  setting,  so  that  “adding  on 
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research  sequences”  is  not  possible.  With  these  limitations  in  mind,  a  growing 
literature  on  DTI  in  TBI  has  begun  to  address  the  DTI  findings  at  different  stages. 

2. 4. 1.1  Acute  Stage 

There  are  conflicting  findings  for  FA  and  ADC  in  mTBI.  Chu  et  at.,38  Wu  et  a/.,39 
and  Wilde  et  at.40  from  Baylor  College  of  Medicine  scanned  10  to  12  adolescents 
with  mTBI  within  6  days  of  injury  and  reported  increased  FA,  reduced  ADC  and 
reduced  radial  ( short  axis )  diffusivity  in  WM  regions  and  left  thalamus.  Similarly, 
Bazarian  et  a/.41  studied  six  mTBI  patients  within  72  hours  and  reported 
increased  FA  in  the  posterior  CC  and  reduced  ADC  in  the  anterior  limb  of  the 
internal  capsule  (IC).  Similar  to  Chu  et  at.,38  and  Wilde  et  at.39  Mayer  et  at.42 
studied  22  patients  with  mTBI  within  12  days  of  injury  and  reported  FA  increase 
and  reduced  radial  diffusivity  in  the  CC  and  left  hemisphere  tracts. 

Inglese  et  at.,43  in  contradiction  to  Chu  et  at.38  Wilde  et  at. , 39,40  Bazarian  et  a/.,41 
and  Mayer  et  at.,42  found  reduced  FA  in  the  splenium  of  CC  and  posterior  limb  of 
IC  in  20  patients  with  mTBI  imaged  up  to  10  days  after  injury  (mean  =  4  days). 
Manually  drawn  regions  of  interest  were  used  to  assay  the  genu  and  splenium  of 
the  CC,  the  centrum  semiovale,  and  the  posterior  limb  IC  bilaterally.  In  the  same 
line  as  Inglese  et  at.,  Arfanakis  and  colleagues44  studied  a  handful  of  patients  with 
mTBI  at  the  acute  stage  and  reported  FA  decrease  in  major  WM  tracts. 

2. 4. 1.2  Subacute  Stage 

All  studies  reported  reduced  FA  and/or  increased  diffusivity ;  i.e.  ADC,  trace  or 
mean  diffusivity  (MD),  at  this  stage.  Messe  et  a/.45  studied  23  mTBI  at  the 
subacute  stage  and  found  significantly  increased  MD  in  mTBI  patients  with 
poor  outcome.45  The  authors  did  not  find  significant  changes  of  FA  values. 
Lipton  et  at.46  scanned  20  mTBI  patients  in  the  subacute  stage  and  demon¬ 
strated  reduced  FA  and  increased  MD  in  frontal  subcortical  WM.  Miles  et  at.41 
studied  17  acute  and  subacute  mTBI  patients  and  found  reduced  FA  and 
increased  MD.  All  of  these  studies  at  the  subacute  stage  reported  a  similar 
profile  of  DTI  measures  except  for  different  locations  of  injury. 

2.4.13  Chronic  Stage 

Three  studies  report  either  reduced  FA  or  increased  diffusivity  or  both.  Cubon 
et  a!.48  studied  10  collegiate  athletes  with  concussion  at  the  chronic  stage  and 
found  MD  increase  in  left  WM  tracts,  internal  capsule,  and  thalamic  acoustic 
radiations.  Lo  et  at.  studied  10  mTBI  patients  and  reported  reduced  FA  and 
increased  ADC  in  the  left  genu  of  the  CC  and  increased  FA  in  internal  capsule 
bilaterally  49  Niogi  et  at.50  studied  43  patients  with  chronic  mTBI  and  reported 
reduced  FA  in  a  large  number  of  WM  tracts. 

All  DTI  studies  of  moderate  to  severe  TBI  patients16'51-53  and  subacute/ 
chronic  mTBI  patients43’50  54  56  report  FA  decreases,  which  are  correlated  with 
clinical  or  neuropsychological  measures.  However,  there  are  seemingly 
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contradictory  findings  in  mTBI  in  the  acute  stage  (within  1  week  after  injury)  in 
the  literature:  Inglese43  and  Arfanakis44  both  reported  FA  decreases,  while 
Wilde40  and  Bazarian41  reported  FA  increases  and  decreased  radial  diffusivity. 
It  has  been  suggested  that  increased  FA  acutely  may  reflect  cytotoxic  edema,40 
which  would  shunt  extracellular  fluid  into  swollen  cells.  This  could  have  the 
effect  of  reducing  inter-axonal  free  water  and  therefore  increasing  anisotropy. 

2. 4. 1.4  Longitudinal  Studies 

Only  a  few  investigators  have  followed  FA  over  time  in  the  same  patients. 
Sidaros  et  al.51  studied  23  patients  with  severe  TBI  at  8  weeks  and  again  at  12 
months  and  found  that  partial  recovery  of  initially  depressed  FA  values  in  the 
internal  capsule  and  centrum  semiovale  predicted  a  favorable  outcome.  Kumar 
et  al.5S  studied  16  patients  with  moderate  to  severe  TBI  at  2  weeks  or  less,  6 
months,  and  2  years  and  found  persistently  reduced  FA  except  in  the  genu  of  the 
corpus  callosum,  where  there  was  partial  normalization  by  2  years.  Recently, 
two  studies  by  Mayer  et  al.42  and  Rutgers  et  al.55  reported  that  FA  may  partially 
normalize,  reflecting  recovery.  This  evidence  suggests  that  a  systematic  inves¬ 
tigation  of  a  large  number  of  patients  with  mTBI  at  acute,  subacute,  and  chronic 
stages  is  warranted  to  reveal  the  evolution  of  pathophysiology  in  mTBI. 

2.4.2  Correlation  Between  DTI-Derived  WM  Injury  Topography 
and  Neuropsychological  Deficits 

Patients  with  mTBI  often  develop  a  constellation  of  physical,  cognitive,  and 
emotional  symptoms  that  are  collectively  known  as  post-concussion  syndrome 
(PCS).  In  terms  of  neurocognitive  symptoms,  there  are  four  key  domains 
implicated  in  chronic  neuropsychological  impairment  after  mTBI.  These 
domains  include:  1)  higher-order  attention,  2)  executive  function,  3)  episodic 
memory,  and  4)  speed  of  information  processing.  To  date,  several  studies  have 
demonstrated  typical  mTBI  cognitive  profiles  and  association  with  DTI  find¬ 
ings.  Niogi  et  al.  summarized  the  topographic  and  neurocognitive  deficits.37 

Damage  to  the  frontal  WM  has  been  reported  to  be  associated  with  impaired 
executive  function.  Lipton  et  al.46  studied  20  acute  to  subacute  patients  and 
reported  that  reduced  FA  in  WM  of  dorsolateral  prefrontal  cortex  (DLPRC)  is 
correlated  with  worse  executive  function.  Frontal  WM  injury  is  also  associated 
with  attention  deficit.  Niogi  et  al.50  reported  that  reduced  FA  in  the  left  anterior 
corona  radiata  is  correlated  with  attention  control  in  chronic  mTBI  patients. 

Injury  at  the  temporal  WM  tracts  or  cingulum  bundle  may  cause  memory 
problems.  Niogi  et  al.  reported  that  reduced  FA  in  uncinate  fasciculus  corre¬ 
lated  with  memory  performance.59  Wu  et  al.29  reported  that  FA  measure  of  left 
cingulum  bundle  correlated  with  delayed  recall. 

Injury  of  the  callosal  fibers  has  been  reported  to  be  associated  with  PCS 
scores.  Wilde  et  al.40  studied  10  adolescent  patients  with  mTBI  in  the  acute 
stage  and  reported  that  increased  FA  and  decreased  ADC  and  MD  in  corpus 
callosum  is  correlated  with  patient  PCS  score.  Bazarian  et  al.41  studied  6 
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patients  with  mTBI  in  the  acute  stage  and  reported  a  lower  mean  trace  in  the 
left  anterior  IC  and  a  higher  FA  in  the  posterior  CC.  FA  values  correlated  with 
patient’s  72  hours  PCS  score  and  visual  motor  speed  and  impulse  control. 

The  overall  burden  or  extent  of  WM  injury  is  associated  with  both  speed  of 
information  processing  and  overall  functional  outcome.  Niogi  et  al.50  studied 
34  subacute  to  chronic  mTBI  patients  and  reported  that  FA  decreased  in 
several  WM  regions,  including  anterior  corona  radiata,  uncinate  fasciculus,  CC 
genu,  and  cingulum  bundle.  They  demonstrated  that  the  number  of  damaged 
WM  regions  is  correlated  with  patient’s  reaction  time.  Miles  et  al.47  studied  17 
mTBI  patients  at  the  acute  stage  and  followed  them  up  to  6  months  after  injury. 
They  reported  that,  at  the  acute  stage,  the  increased  mean  diffusivity  (MD)  in 
centra  semiovale,  the  genu  and  splenium  of  CC,  and  the  posterior  limb  of  IC 
tended  to  correlate  with  patient  response  speed  at  6  months  after  injury. 
Regarding  the  overall  outcome,  Messe  et  al 45  divided  patients  with  mTBI  into 
two  outcome  groups,  poor  outcome  versus  good  outcome.  Poor  outcome 
patients  showed  significantly  higher  mean  diffusivity  (MD)  values  than  both 
controls  and  good  outcome  patients  in  the  corpus  callosum,  the  right  anterior 
thalamic  radiations  and  the  superior  longitudinal  fasciculus,  the  inferior 
longitudinal  fasciculus  and  the  fronto-occipital  fasciculus  bilaterally. 

Interestingly,  injury  or  reduced  blood  supply  in  the  thalamus,  which  is  the 
relay  station  of  neuronal  pathways,  may  cause  a  constellation  of  symptoms  in 
speed  of  information  processing,  memory,  verbal,  and  executive  function. 
Grossman  et  al.  studied  22  patients  with  subacute  to  chronic  mTBI  by  using 
diffusion  kurtosis  imaging,  which  is  a  more  advanced  form  of  diffusion  analysis, 
and  demonstrated  that  overall  cognitive  impairment  is  associated  with  the 
diffusion  measurement  in  the  thalamus  and  internal  capsule.60  This  work  is 
along  the  same  line  of  a  perfusion  study  by  the  same  group,  which  demon¬ 
strated  that  reduced  blood  flow  in  the  thalamus  correlated  with  patient’s  overall 
neurocognitive  function.61 

In  summary,  significant  progress  has  been  made  by  researchers  in  recent 
years  regarding  the  prognostic  value  of  DTI  in  the  form  of  FA  for  neurocog¬ 
nitive  outcome  in  patients  with  mTBI.  There  is  still  a  need  to  evaluate  the 
neural  basis  of  a  patient’s  recovery  process,  which  requires  longitudinal  studies 
over  a  large  number  of  patients. 

2.4.3  Pros  and  Cons  of  Different  Analytical  Approaches 

One  challenge  of  DTI  is  post-processing  and  interpretation  of  the  images. 
There  are  several  approaches  available,  including  whole  brain  histogram  ana¬ 
lysis,  region  of  interest  analysis,  voxel-based  analysis,  and  tract-based  spatial 
statistics,  to  name  a  few.  Each  approach  has  pros  and  cons.  There  is  no  stan¬ 
dardized  procedure  for  DTI  analysis;  it  is  a  moving  target.  The  choice  of 
analytical  approach  depends  on  the  hypothesis  and  goal  of  the  study.  Table  2.2 
lists  the  pros  and  cons  of  each  approach. 

A  whole  brain  histogram  analysis  is  used  to  segment  out  the  whole  brain 
WM  and  plot  the  FA  frequency  distribution.16  Two  studies  have  demonstrated 
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Table  2.2  Comparison  of  different  DTI  analytical  approaches. 


Approach 

Pros 

Cons 

Whole  brain 
histogram 

Easy  to  implement 

Not  sensitive  to  subtle  changes  of 
mild  brain  injury;  lack  of  regional 
injury  information 

5 

Manual  ROI 

Easy  to  implement,  detect 
regional  injury 

Inter-rater  and  intra-rater  variability 
issue;  time  consuming 

Automated  or 

Objective,  consistent, 

Normalization  may  distort  DTI 

semi-auto- 
mated  ROI 

detect  regional  injury 

indices;  potential  misregistration 
issue;  small  fraction  of  voxels  may 
be  averaged  out 

10 

Voxel-based 

Data  driven,  automated, 

Warping  may  cause  DTI  indices’ 

analysis 

objective 

change;  misregistration  cause 
artificial  edging  effect 

Tract  based 

Data  driven,  automated, 

Only  look  at  the  center  line  of  WM 

spatial 

objective,  least  mis- 

tracts;  may  miss  the  target  of 

15 

statistics 

registration  issue. 

lesions  at  other  areas  of  WM 

DTI,  diffusion  tensor  imaging;  ROI,  region  of  interest;  WM,  white  matter. 


that  moderate  to  severe  TBI  patients  had  their  FA  histograms  shifted  toward 
zero,16'51  and  the  TBI  group  had  larger  variations  than  controls.  In  a  mTBI 
study,  Inglese  et  al.43  did  a  whole  brain  histogram  analysis  by  combining  both 
grey  and  white  matter,  and  did  not  find  a  group  difference  between  mTBI 
patients  and  controls.  Instead  of  combining  both  GM  and  WM,  Lipton  et  al 62 
segmented  out  WM  and  demonstrated  the  mTBI  patient  group’s  FA  histogram 
of  WM  only  shifted  to  lower  end.  However,  as  their  data  demonstrated,  it  is  not 
as  sensitive  as  that  of  the  voxel-based  analysis.  The  advantage  of  the  histogram 
analysis  is  that  it  is  relatively  easy  to  implement.  The  disadvantage  is  the  lack  of 
regional  injury  information  and  lack  of  sensitivity  to  subtle  changes  of  the  WM, 
particularly  for  mild  cases.  Unless  the  subject  has  widespread  multifocal  injury 
on  WM  tracts,  a  small  fraction  of  WM  injury  may  not  demonstrate  a  shift  in 
the  histogram  or  change  its  kurtosis. 

To  test  a  hypothesis  of  injury  at  specific  WM  tract(s)  in  association 
with  specific  domain  of  neurocognitive  or  neural  behavioral  symptoms, 
researchers  used  region  of  interest  (ROI)  analysis  by  either  manually 
drawing4143’44'47'49'50'59’63  or  segmenting  out  the  regions  in  a  prior  defined 
template  in  an  automated  or  semi-automated  way.42’64  This  could  also  provide 
needed  lesion  localization  information,  which  should  improve  biomechanical 
models  of  TBI.  Manually  drawn  ROI  is  relatively  easy  and  convenient  to  do. 
One  could  use  a  color-encoded  FA  diffusion  map  as  a  guide  to  define  the 
boundary  of  a  target  region.  It  works  on  the  native  space  of  the  subject’s  brain 
without  worrying  about  the  distortion  caused  by  the  warping  algorithm. 
However,  it  is  time  consuming  and  error  prone.  It  requires  the  rater’s  prior 
knowledge  of  which  specific  tracts  or  locations  to  define  the  region.  It  also  has  a 
potential  inter-rater  and  intra-rater  variability  problem. 

An  automated  ROI  analysis  approach  would  need  initially  to  transform  the 
images  into  a  standard  reference  space;  for  example.  International  Consortium 
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of  Brain  Mapping  (ICBM)  atlas,  by  using  SPM8,  FSL,  or  other  software,  in 
order  to  derive  statistical  maps  computed  from  comparing  “cases”  and  non- 
TBI  controls,  and  then  further  divide  the  WM  into  standardized  regions  that 
are  already  defined  on  the  atlas.  Then  each  regional  FA,  ADC,  or  Eigen  values 
could  be  computed.  In  this  category,  several  studies  reported  significant  FA 
differences  between  mTBI  patients  and  controls.42'64  This  method  overcomes 
the  subjective  issue  in  manually  drawn  ROIs  by  using  a  computer  algorithm  to 
warp  each  subject  to  a  standard  space.  The  result  is  objective  and  consistent. 
However,  the  use  of  a  computer  algorithm  for  image  warping  and  registration 
will  introduce  new  problems.  Due  to  individual  variations,  the  warping  of  one 
image  to  a  standard  space  could  change  the  DTI  indices’  values  from  the  native 
space.  Furthermore,  in  a  relatively  large  region  like  the  splenium  of  corpus 
callosum,  a  small  fraction  of  voxels  with  abnormal  FA  value  could  be  averaged 
out  by  the  surrounding  voxels  in  the  same  region,  rendering  this  method 
insensitive  to  small  lesions. 

Voxel-based  analysis  (VBA)  overcomes  the  sensitivity  issue  of  the  ROI 
approach  by  comparing  patient  and  control  DTI  images  voxel  by  voxel  to 
determine  the  abnormal  voxels  between  groups.  This  approach  is  fully  auto¬ 
mated,  data-driven,  and  reproducible.  Several  studies  reported  different  clusters 
of  WM  tracts  with  abnormal  FA  values  by  using  this  method.38'41,46'62’65’66 
However,  it  suffers  from  the  same  warping  and  misregistration  problems  as  that 
of  the  automated  ROI  approach.  Particularly  at  the  edge  of  WM  or  GM/WM 
junction  area,  any  misregistration  will  result  in  comparing  WM  of  one  group 
with  GM  of  another  group  and  result  in  abnormal  values  of  DTI  indices. 

To  overcome  the  misregistration  problem,  a  relatively  new  approach,  called 
tract-based  spatial  statistics  (TBSS),  has  been  introduced.67  First,  it  averages 
both  patient  and  control  brains  to  form  a  template  brain.  Second,  it  defines  a 
template  skeleton  of  the  WM  tracts  from  the  center  lines  of  WM  of  the  template 


Figure  2.3  A  mTBI  patient  scanned  on  the  day  of  injury  (5.5  hours  after  injury)  and  6 
weeks  later.  A  31 -year-old  male  fell  10  feet  (3.0  m)  off  a  ladder,  striking  the 
back  of  his  head  with  brief  loss  of  consciousness  and  confusion.  The  patient 
developed  persistent  mild  cognitive  symptoms  6  weeks  after  injury.  Note 
the  same  location  of  reduced  FA  in  left  corona  radiata.  Similar  finding  in 
splenium  of  CC  (not  shown).  Global  WM  FA  mean  was  within  normal. 
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Figure  2.4  A  demonstration  of  TBSS  analysis  to  detect  WM  injury  at  the  splenium  of 
corpus  callosum.  A  young  man  involved  in  a  motor  vehicle  accident  with 
negative  structural  MR  images.  TBSS  analysis  demonstrated  significantly 
lower  FA  values  at  the  splenium  of  corpus  callosum,  as  demonstrated  in 
the  highlighted  area.  As  shown  on  the  image,  it  focuses  on  the  skeleton 
(center  lines)  of  WM. 

brain.  Third,  all  of  the  subjects’  WM  skeletons  are  nonlinearly  registered  onto 
the  template  skeleton.  As  such,  the  WM  skeletons  of  all  subjects  will  be  aligned 
at  the  same  space.  Abnormal  DTI  indices  could  be  defined  along  the  skeletons 
by  using  a  group  comparison  on  a  voxel-based  analysis.  Two  groups  reported 
the  use  of  this  method  in  the  detection  of  mTBI45  48  in  association  with  patient 
outcome.  This  method  is  automated  and  reproducible  and  suitable  for  data- 
driven  analysis.  It  has  the  least  problem  of  misregistration.  However,  it  only 
looks  at  the  center  lines  of  WM  tracts.  In  thick  WM  bundles,  like  corpus  cal¬ 
losum,  it  might  miss  lesions  located  off  the  center  lines  of  WM. 


2.4.4  DTI  Caveat 

In  addition  to  the  post-processing  concerns  mentioned  above,  another  caveat  of 
DTI  is  that  it  is  experiment  dependent.  Different  magnet  field  strength  and 
image  acquisition  parameters,  among  others,  will  result  in  different  DTI  mea¬ 
sures.19  Furthermore,  to  date,  most  publications  are  performed  with  group 
comparisons  between  patients  and  controls  instead  of  an  individual-based 
analysis.  However,  radiological  diagnosis  is  case-based.  A  radiologist  or  a 
clinician  needs  to  make  a  decision  for  each  individual  subject  instead  of  group 
comparison.  The  low  sensitivity  of  DTI  to  individual  case  hinders  its  further 
implementation  as  a  routine  clinical  sequence.  A  relatively  standardized  pro¬ 
tocol  is  needed  to  test  a  large  number  of  controls  to  determine  the  normal 
variations  before  deciding  which  specific  case  has  brain  injury. 


1 

5 

10 

15 

20 

25 

30 

35 

40 

45 


Chapter  2 


36 

2.5  In  Vivo  1H  MR  Spectroscopic  Imaging 

Cellular  changes  to  neurons  and  glia  following  TBI  are  complex  and  dynamic. 
Proton  ( 1 H)  MRSI  has  the  advantage  of  measuring  brain  metabolism  in  vivo ,  and 
is  able  to  detect  various  biochemical  processes  of  brain  injury  such  as  the  loss 
(or  dysfunction)  of  neuronal  cells.  TBI  is  known  to  induce  changes  in  N- 
acetylaspartate  (NAA)  (a  neuronal  marker),  choline  (Cho)  (a  marker  for 
membrane  disruption,  synthesis,  or  repair),  and  lactate  (Lac)  (a  measure  of 
anaerobic  glycolysis).68  Many  studies  of  metabolic  changes  after  TBI  have 
proven  MRS  to  be  a  sensitive  tool  to  predict  neurologic,69  74  cognitive,75  79  and 
functional  outcome.80  Significant  decreases  of  NAA  and  increases  of  Cho  have 
been  observed  in  “normal  appearing”  WM  or  GM  following  moderate  to  severe 
TBI.74-76’81-83  Elevated  Cho  may  be  detected  in  WM  as  a  breakdown  product 
after  shearing  of  myelin84  and  the  presence  of  lactate  in  the  brain  suggests 
hypoxic/ischemic  injury.  Researchers  from  Loma  Linda  University  Medical 
Center  found  that  increased  frontal  white  matter  (FWM)  Cho/creatine  and 
decreased  FWM  NAA/Cho  at  acute  stage  are  useful  for  predicting  functional 
outcome  in  patients  with  moderate  to  severe  TBI.  In  addition,  they  also  found 
that  elevated  Glx  and  Cho  in  the  acute  stage  are  sensitive  indicators  of  injury  and 
predictors  of  poor  outcome  in  normal  appearing  brain  of  severe  TBI  patients.71 
They  also  found  that  the  decrease  of  NAA  at  subacute  stage  is  more  predictive  of 
patient  outcome. 

In  mTBI  studies,  there  are  relatively  less  data  reported  in  comparison  with 
moderate  to  severe  TBI  studies.  For  a  chronic  mTBI  study,  Kirov  et  a/.85 
scanned  20  mTBI  patients  and  17  age-  and  gender-matched  controls  and 
reported  that  the  mTBI  group  had  more  variations  in  NAA,  Cre,  and  Cho  than 
the  control  group.  Cohen  et  a/.86  investigated  20  patients  with  mTBI  from  9 
days  to  years  after  injury  and  found  that  whole  brain  NAA  decreased  by  12% 
along  with  significant  global  and  grey  matter  atrophy,  in  comparison  with  19 
controls.  This  further  demonstrates  NAA  as  a  biomarker  of  neuronal  injury  or 
cellular  loss.  Regarding  the  investigations  of  mTBI  in  the  subacute  or  semi¬ 
acute  stages,  a  few  studies  reported  an  NAA/Cr  decrease, 86-88  which  is  sug¬ 
gestive  of  neuronal  loss. 

A  recent  study  on  10  patients  with  mTBI  in  the  semi-acute  stage  showed 
creatine  increases  in  white  matter  and  Glx  (glutamate  and  glutamine)  decreases 
in  grey  matter.  The  creatine  (Cre)  increases  in  white  matter  were  correlated  with 
patient’s  executive  function.89  The  same  research  group  also  studied  a  relatively 
larger  group  of  patients  with  longitudinal  follow-up.  In  their  study  of  30  mTBI 
patients  at  the  first  visit  within  3  weeks  after  injury,  they  reported  Cre  and  Glx 
increases  in  WM  and  Glx  decreases  in  GM.90  Among  the  patients,  17  returned 
for  follow-up  at  the  chronic  stage;  their  GM  Glx  demonstrated  normalization 
and  their  WM  Cre  and  Glx  demonstrated  a  trend  of  normalization.  Of  parti¬ 
cular  note  is  that,  at  the  first  visit,  the  mTBI  group  did  not  have  any  difference 
from  the  controls  in  neuropsychological  assessment  except  for  more  somatic, 
cognitive,  and  emotional  symptoms.90  This  demonstrates  that  MRS  may  be 
sensitive  to  subtle  changes  of  the  brain,  which  may  not  yet  be  manifesting 
significant  neural  cognitive  deficits. 
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As  with  other  MRI  techniques,  MR  spectroscopy  has  also  undergone 
tremendous  advancements  in  the  past  two  decades,  evolving  from  single 
voxel  MR  spectroscopy  to  2D  spectroscopy  imaging  and  3D  whole  brain 
spectroscopy  imaging.  Single  voxel  spectroscopy  (SVS)  allows  acquisition  of 
a  single  spectrum  from  one  volume  element  (voxel)  typically  8  cm3  or  more, 
whereas  2D  or  3D  magnetic  resonance  spectroscopic  imaging  (2D-MRSI/3D- 
MRSI),  also  called  chemical  shift  imaging  (CSI),  allows  for  the  simultaneous 
acquisition  of  multiple  spectra  from  smaller  adjacent  voxels  through  multiple 
sections  of  the  brain.  MRSI  has  an  inherent  advantage  over  SVS  because  it  is 
better  able  to  evaluate  regional  distributions  of  neurochemical  alterations. 
Instead  of  hypothesizing  a  certain  region  or  certain  slice  of  the  brain  to  have 
abnormal  metabolite  signals;  3D  MRSI  allows  investigators  to  search  the 
whole  brain  to  identify  any  area  that  might  have  abnormal  metabolites.  This 
could  be  useful  in  mTBI  due  to  the  fact  that  choosing  an  arbitrary  region  to 
study  would  be  like  searching  for  a  needle  in  a  haystack.  Govind  et  a/.91  used 
a  3D-MRSI  sequence  to  scan  20  mild  and  9  moderate  TBI  patients  and  found 
a  widespread  decrease  of  NAA  and  NAA/Cre,  and  increases  of  Cho  and 
Cho/NAA,  within  all  lobes  of  the  TBI  subject  group,  and  with  the  largest 
differences  seen  in  WM.  Examination  of  the  association  between  all  of  the 
metabolite  measures  and  the  neuropsychological  test  scores  found  the 
strongest  negative  correlations  to  occur  in  the  frontal  lobe  and  for  Cho/ 
NAA. 

One  caveat  of  the  previous  MRS  work  on  TBI  is  that  Cre  has  been  used  as  an 
internal  reference  with  the  assumption  that  Cre  levels  will  remain  constant  after 
injury.9  However,  this  might  not  be  true  in  an  injured  brain,  which  could  have 
Cre  disturbance  after  TBI.  As  demonstrated  by  Gasparovic  et  a/.,89  90  the  Cre 
increase  in  mTBI  patients  at  the  semi-acute  stage  might  give  a  false  impression 
of  decreased  NAA/Cre  level.  It  is  therefore  important  to  quantitate  metabolite 
concentrations  carefully,  particularly  after  mTBI,  to  determine  possible  small 
changes  in  all  metabolites.  Comparing  metabolite  ratios  may  mask  important 
metabolite  changes. 

2.6  Perfusion  Imaging 

After  the  primary  brain  damage  occurs,  even  more  devastating  secondary  brain 
injuries  can  occur.92  This  usually  includes,  but  is  not  limited  to,  ischemic  and 
hypoxic  damage,  cerebral  swelling,  rise  of  intracranial  pressure,  hydrocephalus, 
and  infection93  in  moderate  to  severe  TBI  patients  at  the  acute  stage.  The  injury 
to  the  brain  can  result  in  disturbance  of  cerebral  blood  flow  (CBF),  which 
correlates  with  poor  brain  tissue  oxygenation,94  and  unfavorable  neurological 
outcome,95  which  is  implicated  in  rendering  the  brain  vulnerable  to  secondary 
damage.95  It  has  been  recognized  that  the  failure  of  brain  perfusion  is  the  most 
common  type  of  hypoxic-ischemic  brain  damage,  but  is  difficult  to  identify  even 
in  histological  studies.93  Cerebral  ischemia  (CBF<  18  mF/100  g/min)  is  present 
in  approximately  90%  of  patients  who  die  with  head  injury.92’96  97  Further¬ 
more,  after  TBI,  the  uncoupling  between  CBF  and  brain  metabolic  demand 


1 

5 

10 

15 

20 

25 

30 

35 

40 

45 


38 


Chapter  2 


occurs,  whereby  local  cerebral  glucose  demand  increases  and  local  CBF  is 
disturbed  due  to  damage  to  the  vascular  muscle  tone  and  arterioles.98  The  early 
hyper-glycolysis  is  thought  to  be  due  to  the  cellular  efforts  to  reestablish  ionic 
gradients.95  Different  aspects  of  ischemia  and/or  edema  after  injury  can  be 
quantified  using  perfusion  imaging,  including  perfusion  weighted  imaging 
(PWI)  and  arterial  spin  labeling  (ASL),  and  DW1/DT1.  Kim  et  al.17  used  ASL 
to  measure  the  resting  state  CBF  in  28  moderate  to  severe  TBI  at  chronic  stage 
and  found  a)  a  global  decrease  of  CBF  in  TBI  patients  than  that  in  controls, 
and  b)  prominent  regional  hypoperfusion  in  the  posterior  cingulate  cortices,  the 
thalami,  and  multiple  locations  in  the  frontal  cortices. 

However,  in  the  population  of  mTBI,  how  cerebral  perfusion  is  being 
affected  in  patients  and  their  association  with  patients’  neurocognitive  status  is 
still  less  understood.  Given  the  milder  severity  in  mTBI  patients,  perfusion 
imaging  at  the  acute  stage  could  be  meaningful  both  for  clinical  decision¬ 
making  and  patients’  outcome  prediction.  Ge  et  al.61  also  employed  ASL  to 
study  21  mTBI  patients  at  the  chronic  stage  and  demonstrated  reduced  CBF  in 
both  sides  of  the  thalamus,  which  is  correlated  with  the  patients’  speed  of 
information  processing,  memory,  verbal,  and  executive  function. 

In  summary,  our  knowledge  is  very  limited  in  the  understanding  of  cerebral 
blood  flow  and  its  impact  to  mTBI  patients.  Perfusion  imaging  of  mTBI, 
particularly  at  the  acute  stage,  might  shed  new  light  in  our  understanding  of 
this  disease. 


2.7  Delivering  Improved  Care  to  Patients  with  mTBI 

Ijg]  Despite  the  fact  that  the  ED  sees  the  majority  of  patients  with  mTBI,"  most 
stay  only  a  few  hours  and  then  are  discharged  home.  After  that,  most 
patients  fail  to  be  followed  up.  The  acute  stage,  within  24  hours  after  injury, 
is  the  critical  time  point  for  imaging  to  deliver  a  real  impact  to  medicine.100 
Either  improved  detection  or  outcome  prediction  will  greatly  help  emergency 
physicians  develop  a  better  referral  pattern  or  management  plan  for  the 
patient  or  notify  the  patient’s  family  members  to  arrange  necessary  resources 
in  advance.  However,  most  of  the  patients  with  mTBI  at  the  acute  stage  do 
not  get  an  MRI  scan  due  to  the  high  cost  for  medical  insurance  and  the  lack 
of  accessibility  to  an  MRI  magnet.  To  date,  very  few  studies  are  designed  to 
target  this  critical  time  point.  There  is  an  urgent  need  for  a  comprehensive 
use  of  advanced  MRI  techniques  to  evaluate  the  patients  from  the  acute 
setting  (within  24  hours)  to  the  chronic  stage  after  injury  to  identify  the 
association  between  injury  pathology  and  outcome.  In  the  future,  with  the 
availability  of  MR  magnets  in  the  ED  and  reduced  costs  for  an  MRI 
scan,  routinely  scanning  mTBI  patients  at  the  acute  stage  will  be  possible. 
However,  to  determine  who  will  get  MRI  scan  in  ED  will  be  another  issue, 
due  to  cost.  In  conjunction  with  other  biochemical  markers,  a  policy  on 
clinical  decision  making  would  be  expected  to  define  a  subgroup  of  patients 
for  MRI  scan. 
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MRI  has  demonstrated  superior  capabilities  over  CT  in  the  detection  of  subtle 
changes  of  the  brain  after  mTBI.  As  an  advanced  MRI  technique,  DTI  can 
detect  white  matter  abnormalities  that  are  unseen  in  structural  MRI  and  fur¬ 
ther  map  it  with  a  patient’s  specific  domain  of  neuropsychological  symptoms. 
Other  advanced  MRI  techniques  also  demonstrate  the  potential  to  detect  other 
types  of  brain  pathologies,  including  SWI  for  microhemorrhage  detection, 
MRS  for  measuring  abnormal  cerebral  metabolites,  and  perfusion  imaging  for 
detecting  cerebral  blood  flow  disturbances.  In  the  future,  these  advanced  MRI 
techniques  should  be  used  in  a  comprehensive  way  in  a  large  cohort  of  patients 
to  provide  a  panorama  view  of  brain  pathologies.  The  MRI  investigations  at 
the  acute  stage,  within  24  hours  after  injury,  will  most  likely  impact  emergency 
medicine,  which  is  at  the  forefront  of  mTBI  care.  Furthermore,  conjunctional 
use  of  imaging  with  other  biological  markers,  such  as  serum  biomarkers,  will 
likely  help  screen  patients  who  really  need  MRI  in  the  acute  setting. 
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